Pub Date : 2025-12-20DOI: 10.1016/j.chemgeo.2025.123217
Shaoji Yang , Pan Wang , Haijin Xu , Fan Feng
Subduction zone carbon cycling governs the carbon budget balance between Earth's interior and surficial systems and constrains long-term climatic-environmental evolution. While extensive research has focused on oceanic subduction zones, CO2 release mechanisms in continental subduction zones remain poorly understood. The Sulu orogen, which experienced cold subduction followed by hot exhumation and extensive anatexis, provides an ideal natural laboratory for investigating carbon release processes in continental subduction zones. Petrology, mineral chemistry, and titanite geochronology were conducted to investigate a melt-metacarbonate interaction profile, where the migmatite – contact zone (calc-silicate rock) - marble system clearly records CO2 release and migration. Petrographic and mineral chemical analysis reveal diverse metasomatic reactions between silicate melts and marble. In the innermost marble, limited melt infiltration led to olivine and phlogopite formation with CO2 release. At the contact zone, intense metasomatic reactions produced diopside, tremolite, scapolite, and titanite, releasing substantial CO2. The metasomatism also induced carbon migration through calcite-bearing melt/fluid, including calcite-bearing silicate melts near the migmatite, quartz + calcite ± tremolite veins and calcite veins in the contact zone, and carbon-bearing serpentine veins in the marble. Titanite U-Pb age (206 ± 0.7 Ma) and Zr-in-titanite thermometer (T = 792–862 °C) suggest the CO2 release occurred at the exhumation stage. These findings demonstrate that metamorphic decarbonation occurs not only in hot subduction zones but also during hot exhumation of cold subducted slabs, providing new insights into carbon cycling in continental subduction zones and helping address carbon budget imbalance in global carbon cycling models.
{"title":"Melt‑carbonate metasomatism in the Sulu orogen: Implications for CO2 release in the continental subduction zone","authors":"Shaoji Yang , Pan Wang , Haijin Xu , Fan Feng","doi":"10.1016/j.chemgeo.2025.123217","DOIUrl":"10.1016/j.chemgeo.2025.123217","url":null,"abstract":"<div><div>Subduction zone carbon cycling governs the carbon budget balance between Earth's interior and surficial systems and constrains long-term climatic-environmental evolution. While extensive research has focused on oceanic subduction zones, CO<sub>2</sub> release mechanisms in continental subduction zones remain poorly understood. The Sulu orogen, which experienced cold subduction followed by hot exhumation and extensive anatexis, provides an ideal natural laboratory for investigating carbon release processes in continental subduction zones. Petrology, mineral chemistry, and titanite geochronology were conducted to investigate a melt-metacarbonate interaction profile, where the migmatite – contact zone (calc-silicate rock) - marble system clearly records CO<sub>2</sub> release and migration. Petrographic and mineral chemical analysis reveal diverse metasomatic reactions between silicate melts and marble. In the innermost marble, limited melt infiltration led to olivine and phlogopite formation with CO<sub>2</sub> release. At the contact zone, intense metasomatic reactions produced diopside, tremolite, scapolite, and titanite, releasing substantial CO<sub>2</sub>. The metasomatism also induced carbon migration through calcite-bearing melt/fluid, including calcite-bearing silicate melts near the migmatite, quartz + calcite ± tremolite veins and calcite veins in the contact zone, and carbon-bearing serpentine veins in the marble. Titanite U-Pb age (206 ± 0.7 Ma) and Zr-in-titanite thermometer (<em>T</em> = 792–862 °C) suggest the CO<sub>2</sub> release occurred at the exhumation stage. These findings demonstrate that metamorphic decarbonation occurs not only in hot subduction zones but also during hot exhumation of cold subducted slabs, providing new insights into carbon cycling in continental subduction zones and helping address carbon budget imbalance in global carbon cycling models.</div></div>","PeriodicalId":9847,"journal":{"name":"Chemical Geology","volume":"702 ","pages":"Article 123217"},"PeriodicalIF":3.6,"publicationDate":"2025-12-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145796051","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 : 2025-12-18DOI: 10.1016/j.chemgeo.2025.123210
V.B. Rivera Banuchi , J. Gong , E.N. Mansbach , T. Bosak , J.A. Hurowitz
Aqueous iron UV photooxidation, a proposed pathway for early Martian iron deposits produced in surface environments, produces hydrogen gas, a potentially important greenhouse gas for early Mars. This experimental study focused on iron photooxidation in the presence of dissolved inorganic carbon (DIC) at a range of concentrations that represent partial pressures of CO2 (0.1–1 bar), consistent with estimated atmospheric pressures and compositions on early Mars. We report in-situ sampling of gas headspace to understand the extent to which DIC impacts the production of hydrogen gas that may be produced from aqueous iron photooxidation. The experiments demonstrate that ferrous iron photooxidation does not occur at any of the tested DIC concentrations (3.4–35 mM). Ferrous iron carbonate minerals were precipitated at pH conditions above 6, while maintaining relatively high dissolved iron concentrations at siderite saturation. We observed hydrogen and methane gas production at varying initial conditions, but no clear relationships to either iron loss to oxidation and/or precipitation or DIC concentrations could be inferred. Instead, we suggest that hydrogen is a byproduct of methane photolysis. The methane itself was sourced from an unknown background organic carbon compound from residual organic impurities in the deionized water source utilized in the experiments and is unrelated to iron photochemistry. These observations indicate that iron photooxidation is an unlikely mechanism to have produced hydrogen and precipitated iron oxide minerals at near-neutral pH in the presence of 0.1–1 bar CO2 on early Mars.
{"title":"Photochemistry of Fe (II) and carbonate-bearing waters and the influence on Greenhouse Gas production in Early Mars","authors":"V.B. Rivera Banuchi , J. Gong , E.N. Mansbach , T. Bosak , J.A. Hurowitz","doi":"10.1016/j.chemgeo.2025.123210","DOIUrl":"10.1016/j.chemgeo.2025.123210","url":null,"abstract":"<div><div>Aqueous iron UV photooxidation, a proposed pathway for early Martian iron deposits produced in surface environments, produces hydrogen gas, a potentially important greenhouse gas for early Mars. This experimental study focused on iron photooxidation in the presence of dissolved inorganic carbon (DIC) at a range of concentrations that represent partial pressures of CO<sub>2</sub> (0.1–1 bar), consistent with estimated atmospheric pressures and compositions on early Mars. We report in-situ sampling of gas headspace to understand the extent to which DIC impacts the production of hydrogen gas that may be produced from aqueous iron photooxidation. The experiments demonstrate that ferrous iron photooxidation does not occur at any of the tested DIC concentrations (3.4–35 mM). Ferrous iron carbonate minerals were precipitated at pH conditions above 6, while maintaining relatively high dissolved iron concentrations at siderite saturation. We observed hydrogen and methane gas production at varying initial conditions, but no clear relationships to either iron loss to oxidation and/or precipitation or DIC concentrations could be inferred. Instead, we suggest that hydrogen is a byproduct of methane photolysis. The methane itself was sourced from an unknown background organic carbon compound from residual organic impurities in the deionized water source utilized in the experiments and is unrelated to iron photochemistry. These observations indicate that iron photooxidation is an unlikely mechanism to have produced hydrogen and precipitated iron oxide minerals at near-neutral pH in the presence of 0.1–1 bar CO<sub>2</sub> on early Mars.</div></div>","PeriodicalId":9847,"journal":{"name":"Chemical Geology","volume":"703 ","pages":"Article 123210"},"PeriodicalIF":3.6,"publicationDate":"2025-12-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145786145","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 : 2025-12-18DOI: 10.1016/j.chemgeo.2025.123209
Qasid Ahmad , Hongjie Zhang , Sune G. Nielsen
Precise and accurate stable barium isotope (δ138/134Ba relative to the NIST3104a standard) measurements are crucial for geochemical studies, as natural variations can be small. Double spike Ba isotope analysis by MC-ICP-MS is typically the preferred method for generating reliable Ba isotope data as it effectively corrects for instrumental mass fractionation and isotope fractionation during chemical separation, while also providing precise concentrations via isotope dilution. However, measurement of 50 ng or less of Ba using high sensitivity interface cones can be challenging due to the effects of isobaric interferences from trace Xe in the Ar plasma. Here we demonstrate that interferences from Ba-hydrides (BaH+) also have significant effects on double spike corrected δ138/134Ba-values and, therefore, must be assessed and corrected.
We present a multi-isotope deconvolution strategy utilizing a 135Ba–136Ba double spike and three independent inversions using different un-spiked isotope pairs (134Ba–137Ba, 134Ba–138Ba, and 137Ba–138Ba) to solve the double spike equation. Each inversion independently calculates δ138/134Ba that is differently affected by both background Xe and hydride interferences during the measurement. We show that the Xe correction is smallest for the 137Ba–138Ba inversion, which allows accurate Ba isotope data to be obtained for this inversion even when <5 ng of Ba is consumed during triplicate sample measurements. However, the 137Ba–138Ba inversion also carries the largest effect from hydride formation, which we correct for using a combined 134Ba–138Ba/137Ba–138Ba inversion method. This approach produces long-term external reproducibility for rock reference materials of ±0.025 ‰ (2SD) when consuming less than 50 ng of Ba, which represents an improvement over previously reported Ba double-spike MC-ICP-MS methods at similar Ba consumption. We conclude that the double spike multi-isotope deconvolution approach can accurately assess and correct for elemental and hydride interferences for isotope systems with 5 or more stable isotopes suggesting that isotope systems such as Zr, Mo, Cd, Sn and Te might benefit from this strategy.
{"title":"A multi-isotope deconvolution approach to correct for hydride and xenon interferences during double spike barium isotope analyses with MC-ICP-MS","authors":"Qasid Ahmad , Hongjie Zhang , Sune G. Nielsen","doi":"10.1016/j.chemgeo.2025.123209","DOIUrl":"10.1016/j.chemgeo.2025.123209","url":null,"abstract":"<div><div>Precise and accurate stable barium isotope (δ<sup>138/134</sup>Ba relative to the NIST3104a standard) measurements are crucial for geochemical studies, as natural variations can be small. Double spike Ba isotope analysis by MC-ICP-MS is typically the preferred method for generating reliable Ba isotope data as it effectively corrects for instrumental mass fractionation and isotope fractionation during chemical separation, while also providing precise concentrations via isotope dilution. However, measurement of 50 ng or less of Ba using high sensitivity interface cones can be challenging due to the effects of isobaric interferences from trace Xe in the Ar plasma. Here we demonstrate that interferences from Ba-hydrides (BaH<sup>+</sup>) also have significant effects on double spike corrected δ<sup>138/134</sup>Ba-values and, therefore, must be assessed and corrected.</div><div>We present a multi-isotope deconvolution strategy utilizing a <sup>135</sup>Ba–<sup>136</sup>Ba double spike and three independent inversions using different un-spiked isotope pairs (<sup>134</sup>Ba–<sup>137</sup>Ba, <sup>134</sup>Ba–<sup>138</sup>Ba, and <sup>137</sup>Ba–<sup>138</sup>Ba) to solve the double spike equation. Each inversion independently calculates δ<sup>138/134</sup>Ba that is differently affected by both background Xe and hydride interferences during the measurement. We show that the Xe correction is smallest for the <sup>137</sup>Ba–<sup>138</sup>Ba inversion, which allows accurate Ba isotope data to be obtained for this inversion even when <5 ng of Ba is consumed during triplicate sample measurements. However, the <sup>137</sup>Ba–<sup>138</sup>Ba inversion also carries the largest effect from hydride formation, which we correct for using a combined <sup>134</sup>Ba–<sup>138</sup>Ba/<sup>137</sup>Ba–<sup>138</sup>Ba inversion method. This approach produces long-term external reproducibility for rock reference materials of ±0.025 ‰ (2SD) when consuming less than 50 ng of Ba, which represents an improvement over previously reported Ba double-spike MC-ICP-MS methods at similar Ba consumption. We conclude that the double spike multi-isotope deconvolution approach can accurately assess and correct for elemental and hydride interferences for isotope systems with 5 or more stable isotopes suggesting that isotope systems such as Zr, Mo, Cd, Sn and Te might benefit from this strategy.</div></div>","PeriodicalId":9847,"journal":{"name":"Chemical Geology","volume":"703 ","pages":"Article 123209"},"PeriodicalIF":3.6,"publicationDate":"2025-12-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145785270","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 : 2025-12-18DOI: 10.1016/j.chemgeo.2025.123208
Arka Chakraborty , Sourav Ganguly , Arild Folkvord , Steven E. Campana , Prosenjit Ghosh , Ramananda Chakrabarti
We report stable (δ88/86Sr) and radiogenic (87Sr/86Sr) Sr isotopic compositions of natural and laboratory-cultured aragonitic fish otoliths. The natural samples are from five different species collected from the Atlantic Ocean that grew over temperature ranges of 2–25 °C while the cultured samples grew in water temperatures of 4–16 °C. The otoliths exhibit significant variability in δ88/86Sr values (−0.095 ‰ to 0.102 ‰) and incorporate the lighter 86Sr isotope relative to ambient seawater in both natural and experimental settings. The δ88/86Sr values in natural otoliths negatively correlate with ambient seawater temperature. They show positive and negative correlations with previously reported Δ47 and δ44/40Ca values, respectively, measured on the same powdered samples. A multi-proxy approach to determine seawater temperature combining δ88/86Sr, δ44/40Ca, and Δ47 values and Sr/Ca ratio is proposed, which has an uncertainty of ±0.75 °C and is better than uncertainties when these proxies are used in isolation. Strong correlations are also observed between Sr elemental partitioning (KDSr/Ca) and Δ88/86Sr (δ88/86Srotolith - δ88/86Srseawater) values in natural otoliths. Although the lighter isotopes of both Sr and Ca are preferentially incorporated in the otolith carbonates, a negative correlation is observed between δ88/86Sr and δ44/40Ca values. Contrasting correlations are also observed between KDSr/Ca and Δ88/86Sr values as well as KDSr/Ca and Δ44/40Ca (δ44/40Caotolith - δ44/40Caseawater) values, suggesting different carbonate crystallization pathways for Sr and Ca isotopes.
The δ88/86Sr values of cultured otoliths overlap with natural otolith samples but their KDSr/Ca values are lower than the natural samples. While no relationship is observed between δ88/86Sr values, water temperatures, and KDSr/Ca values, growth rate shows a positive correlation with water temperature in cultured otoliths. The Δ88/86Sr* values (Δ88/86Sr normalized to growth rate) of the cultured otoliths show positive and negative correlations with KDSr/Ca* (KDSr/Ca corrected for growth rate) values and temperature, respectively, consistent with observations from natural otoliths. The δ88/86Sr values do not show any correlation with 87Sr/86Sr ratios. Variability in 87Sr/86Sr ratios in the natural otoliths reflect their provenance from coastal oceans, which receive varying amounts of continent-derived freshwater.
{"title":"δ88/86Sr and 87Sr/86Sr variability in fish otoliths: Implications for seawater temperature, growth-rate, and provenance","authors":"Arka Chakraborty , Sourav Ganguly , Arild Folkvord , Steven E. Campana , Prosenjit Ghosh , Ramananda Chakrabarti","doi":"10.1016/j.chemgeo.2025.123208","DOIUrl":"10.1016/j.chemgeo.2025.123208","url":null,"abstract":"<div><div>We report stable (δ<sup>88/86</sup>Sr) and radiogenic (<sup>87</sup>Sr/<sup>86</sup>Sr) Sr isotopic compositions of natural and laboratory-cultured aragonitic fish otoliths. The natural samples are from five different species collected from the Atlantic Ocean that grew over temperature ranges of 2–25 °C while the cultured samples grew in water temperatures of 4–16 °C. The otoliths exhibit significant variability in δ<sup>88/86</sup>Sr values (−0.095 ‰ to 0.102 ‰) and incorporate the lighter <sup>86</sup>Sr isotope relative to ambient seawater in both natural and experimental settings. The δ<sup>88/86</sup>Sr values in natural otoliths negatively correlate with ambient seawater temperature. They show positive and negative correlations with previously reported Δ<sub>47</sub> and δ<sup>44/40</sup>Ca values, respectively, measured on the same powdered samples. A multi-proxy approach to determine seawater temperature combining δ<sup>88/86</sup>Sr, δ<sup>44/40</sup>Ca, and Δ<sub>47</sub> values and Sr/Ca ratio is proposed, which has an uncertainty of ±0.75 °C and is better than uncertainties when these proxies are used in isolation. Strong correlations are also observed between Sr elemental partitioning (K<sub>D</sub><sup>Sr/Ca</sup>) and Δ<sup>88/86</sup>Sr (δ<sup>88/86</sup>Sr<sub>otolith</sub> - δ<sup>88/86</sup>Sr<sub>seawater</sub>) values in natural otoliths. Although the lighter isotopes of both Sr and Ca are preferentially incorporated in the otolith carbonates, a negative correlation is observed between δ<sup>88/86</sup>Sr and δ<sup>44/40</sup>Ca values. Contrasting correlations are also observed between K<sub>D</sub><sup>Sr/Ca</sup> and Δ<sup>88/86</sup>Sr values as well as K<sub>D</sub><sup>Sr/Ca</sup> and Δ<sup>44/40</sup>Ca (δ<sup>44/40</sup>Ca<sub>otolith</sub> - δ<sup>44/40</sup>Ca<sub>seawater</sub>) values, suggesting different carbonate crystallization pathways for Sr and Ca isotopes.</div><div>The δ<sup>88/86</sup>Sr values of cultured otoliths overlap with natural otolith samples but their K<sub>D</sub><sup>Sr/Ca</sup> values are lower than the natural samples. While no relationship is observed between δ<sup>88/86</sup>Sr values, water temperatures, and K<sub>D</sub><sup>Sr/Ca</sup> values, growth rate shows a positive correlation with water temperature in cultured otoliths. The Δ<sup>88/86</sup>Sr* values (Δ<sup>88/86</sup>Sr normalized to growth rate) of the cultured otoliths show positive and negative correlations with K<sub>D</sub><sup>Sr/Ca</sup>* (K<sub>D</sub><sup>Sr/Ca</sup> corrected for growth rate) values and temperature, respectively, consistent with observations from natural otoliths. The δ<sup>88/86</sup>Sr values do not show any correlation with <sup>87</sup>Sr/<sup>86</sup>Sr ratios. Variability in <sup>87</sup>Sr/<sup>86</sup>Sr ratios in the natural otoliths reflect their provenance from coastal oceans, which receive varying amounts of continent-derived freshwater.</div></div>","PeriodicalId":9847,"journal":{"name":"Chemical Geology","volume":"702 ","pages":"Article 123208"},"PeriodicalIF":3.6,"publicationDate":"2025-12-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145785757","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 : 2025-12-18DOI: 10.1016/j.chemgeo.2025.123207
Kyle M. Henderson , James R. Clark , Anthony E. Williams-Jones , Ryan Sharpe , Mostafa Fayek
Hyper-enriched black shales (HEBS) are thin stratiform and stratabound sulfide-rich (∼70 vol%) units that have percentage level concentrations of metals such as Ni, Mo and Zn. In the Peel River area of Yukon, Canada, HEBS are highly enriched in Ni (∼3 wt%), Zn (∼0.2 wt%), and Mo (0.2 wt%), hosted mainly by millerite (NiS), sphalerite (ZnS) and pyrite (FeS2). The processes that control millerite-sphalerite precipitation in HEBS are poorly constrained. Here, we report the first in-situ sulfur isotope measurements of sulfide minerals in the Peel River HEBS. Early diagenetic pyrite (δ34S -31.9 ± 9.0 ‰) formed in response to open-system bacterial sulfate reduction (BSR), whereas later sphalerite (16.5 ± 5.1 ‰) and millerite (+10.0 ± 1.4 ‰) crystallized via thermochemical sulfate reduction. We have developed a three-stage model for HEBS formation. Stage 1 involves high primary productivity, a redox-stratified water column, a low sedimentation rate, and an active Fe-oxide shuttle, which promoted BSR-driven pyrite precipitation. Organic matter scavenged Ni and Zn from the water column and pore-water, producing metal-rich kerogen. Stages 2 and 3 occur during late diagenesis, with the infiltration of an acidic, oxidized, metal-poor, sulfate-bearing brine that interacted with the metal-bearing kerogen, triggering thermochemical sulfate reduction and millerite-sphalerite precipitation. Subsequent fracturing and increased fluid-flow during peak hydrocarbon generation caused acid-driven replacement of sphalerite by a second generation of millerite. This model accounts for HEBS formation in both Canada and China, demonstrating that, although late-stage fluids did not transport the metals, brines were essential to their genesis.
{"title":"From bacterial to thermochemical sulfate reduction: Sulfur isotope constraints on the genesis of hyper-enriched black shales","authors":"Kyle M. Henderson , James R. Clark , Anthony E. Williams-Jones , Ryan Sharpe , Mostafa Fayek","doi":"10.1016/j.chemgeo.2025.123207","DOIUrl":"10.1016/j.chemgeo.2025.123207","url":null,"abstract":"<div><div>Hyper-enriched black shales (HEBS) are thin stratiform and stratabound sulfide-rich (∼70 vol%) units that have percentage level concentrations of metals such as Ni, Mo and Zn. In the Peel River area of Yukon, Canada, HEBS are highly enriched in Ni (∼3 wt%), Zn (∼0.2 wt%), and Mo (0.2 wt%), hosted mainly by millerite (NiS), sphalerite (ZnS) and pyrite (FeS<sub>2</sub>). The processes that control millerite-sphalerite precipitation in HEBS are poorly constrained. Here, we report the first in-situ sulfur isotope measurements of sulfide minerals in the Peel River HEBS. Early diagenetic pyrite (δ<sup>34</sup>S -31.9 ± 9.0 ‰) formed in response to open-system bacterial sulfate reduction (BSR), whereas later sphalerite (16.5 ± 5.1 ‰) and millerite (+10.0 ± 1.4 ‰) crystallized via thermochemical sulfate reduction. We have developed a three-stage model for HEBS formation. Stage 1 involves high primary productivity, a redox-stratified water column, a low sedimentation rate, and an active Fe-oxide shuttle, which promoted BSR-driven pyrite precipitation. Organic matter scavenged Ni and Zn from the water column and pore-water, producing metal-rich kerogen. Stages 2 and 3 occur during late diagenesis, with the infiltration of an acidic, oxidized, metal-poor, sulfate-bearing brine that interacted with the metal-bearing kerogen, triggering thermochemical sulfate reduction and millerite-sphalerite precipitation. Subsequent fracturing and increased fluid-flow during peak hydrocarbon generation caused acid-driven replacement of sphalerite by a second generation of millerite. This model accounts for HEBS formation in both Canada and China, demonstrating that, although late-stage fluids did not transport the metals, brines were essential to their genesis.</div></div>","PeriodicalId":9847,"journal":{"name":"Chemical Geology","volume":"702 ","pages":"Article 123207"},"PeriodicalIF":3.6,"publicationDate":"2025-12-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145785268","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 : 2025-12-17DOI: 10.1016/j.chemgeo.2025.123194
Georg Spiekermann , Maximilian Schulze
Ab initio molecular dynamics (AIMD) simulations of melt and glass of MgSiO3 and of other compositions at high pressure consistently observe the predominance of fivefold-coordinated silicon ( ) over and in the pressure range of the mantle transition zone. However, the common approach to compute the abundance of different coordination states is only based on a static ‘sequential-3D’ structural analysis. The absence of a time scale in the analysis has led to concern about the validity of its results. To validate AIMD results, here we investigate the coordination dynamics beyond the static limit, by carrying out a dynamical coordination species analysis in ‘4D’ on the basis of the short-term part of the intermittent coordination autocorrelation function described in an earlier study. We found that the predominance of SiO5 species at pressures of the transition zone is persistent even in a dynamical analysis. However, towards lower pressures and most pronounced at ambient pressure, the majority of all SiO5 detected by the common static sequential-3D analysis are so short-lived that they should be considered ‘spurious’ rather than chemical species.
{"title":"The predominance of SiO5 species in MgSiO3 melt at transition zone pressures from a dynamics point of view","authors":"Georg Spiekermann , Maximilian Schulze","doi":"10.1016/j.chemgeo.2025.123194","DOIUrl":"10.1016/j.chemgeo.2025.123194","url":null,"abstract":"<div><div>Ab initio molecular dynamics (AIMD) simulations of melt and glass of MgSiO<sub>3</sub> and of other compositions at high pressure consistently observe the predominance of fivefold-coordinated silicon ( <figure><img></figure> ) over <figure><img></figure> and <figure><img></figure> in the pressure range of the mantle transition zone. However, the common approach to compute the abundance of different coordination states is only based on a static ‘sequential-3D’ structural analysis. The absence of a time scale in the analysis has led to concern about the validity of its results. To validate AIMD results, here we investigate the coordination dynamics beyond the static limit, by carrying out a dynamical coordination species analysis in ‘4D’ on the basis of the short-term part of the intermittent coordination autocorrelation function described in an earlier study. We found that the predominance of SiO<sub>5</sub> species at pressures of the transition zone is persistent even in a dynamical analysis. However, towards lower pressures and most pronounced at ambient pressure, the majority of all SiO<sub>5</sub> detected by the common static sequential-3D analysis are so short-lived that they should be considered ‘spurious’ rather than chemical species.</div></div>","PeriodicalId":9847,"journal":{"name":"Chemical Geology","volume":"702 ","pages":"Article 123194"},"PeriodicalIF":3.6,"publicationDate":"2025-12-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145785274","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 : 2025-12-16DOI: 10.1016/j.chemgeo.2025.123206
Jing Gao , Xueyin Yuan , Xinyang Li , Chaowen Xu , Yue Wang , Wanghua Wu
Geophysical data has revealed anomalous water enrichment in the mantle transition zone (MTZ, 410–660 km depths), which is widely attributed to crustal water released from stagnant subducted slabs. However, the majority of hydrous minerals dehydrate completely at depths <∼300 km, leaving the mechanism of water transport to the deeper mantle unresolved. The discovery of ‘epidote-phases’ in superdeep diamonds from the MTZ offers a potential solution, yet critical uncertainties persist regarding its stability under subduction-zone thermal gradients and MTZ conditions. This study experimentally demonstrates that the hydroxyl groups in epidote strengthen along cold-to-hot subduction geotherms (5–25 °C/km) across ∼250 km. In the deeper reduced mantle (buffered by FeO-Fe), epidote undergoes complete dehydration at ∼250 km via decomposition into grossular, andradite and wustite. Yet, kinetic barriers and the limited extent of reduced domains within slabs probably inhibit this process. Under oxidized conditions (Fe2O3-Fe3O4 buffer), epidote destabilizes at thermally perturbed slab interfaces, and is more likely to decompose via prolonged mantle-derived heating in the MTZ. This breakdown partially traps water in the newly formed grossular garnet. These findings redefine epidote as a dynamic agent in the hydration of the MTZ, offering a novel mineralogical perspective on the water-rich characteristics of this enigmatic region.
{"title":"Epidote emerges as a pivotal player in the hydration dynamics of the mantle transition zone","authors":"Jing Gao , Xueyin Yuan , Xinyang Li , Chaowen Xu , Yue Wang , Wanghua Wu","doi":"10.1016/j.chemgeo.2025.123206","DOIUrl":"10.1016/j.chemgeo.2025.123206","url":null,"abstract":"<div><div>Geophysical data has revealed anomalous water enrichment in the mantle transition zone (MTZ, 410–660 km depths), which is widely attributed to crustal water released from stagnant subducted slabs. However, the majority of hydrous minerals dehydrate completely at depths <∼300 km, leaving the mechanism of water transport to the deeper mantle unresolved. The discovery of ‘epidote-phases’ in superdeep diamonds from the MTZ offers a potential solution, yet critical uncertainties persist regarding its stability under subduction-zone thermal gradients and MTZ conditions. This study experimentally demonstrates that the hydroxyl groups in epidote strengthen along cold-to-hot subduction geotherms (5–25 °C/km) across ∼250 km. In the deeper reduced mantle (buffered by FeO-Fe), epidote undergoes complete dehydration at ∼250 km via decomposition into grossular, andradite and wustite. Yet, kinetic barriers and the limited extent of reduced domains within slabs probably inhibit this process. Under oxidized conditions (Fe<sub>2</sub>O<sub>3</sub>-Fe<sub>3</sub>O<sub>4</sub> buffer), epidote destabilizes at thermally perturbed slab interfaces, and is more likely to decompose via prolonged mantle-derived heating in the MTZ. This breakdown partially traps water in the newly formed grossular garnet. These findings redefine epidote as a dynamic agent in the hydration of the MTZ, offering a novel mineralogical perspective on the water-rich characteristics of this enigmatic region.</div></div>","PeriodicalId":9847,"journal":{"name":"Chemical Geology","volume":"702 ","pages":"Article 123206"},"PeriodicalIF":3.6,"publicationDate":"2025-12-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145766047","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 : 2025-12-13DOI: 10.1016/j.chemgeo.2025.123205
Ruirui Meng , Zuozhen Han , Wen Nie , Chao Han , Yu Han , Maurice E. Tucker
Previous studies have demonstrated that bacteria can induce proto-dolomite precipitation at ambient temperatures. Yet, the mineralization of the bacteria during this process has rarely been reported. Here we describe mineralized bacteria formed during proto-dolomite precipitation mediated by the moderate halophilic bacteria Bacillus licheniformis Y1. The XRD data and SEM images show that proto-dolomite precipitated mainly as aggregates in media with Mg/Ca ratios of 6, 8 and 10 cultured at temperatures of 30 °C and 40 °C, but this was not observed in solutions cultured at 20 °C. Solution chemistry analysis showed that all solutions maintained a higher CO32− concentration and correspondingly higher saturation indices for proto-dolomite at 40 °C. The HRTEM results illustrate that the crystal structures and Mg/Ca atomic ratio of the proto-dolomite more closely match those of dolomite, and some superstructure reflections appear in the HRTEM and SAED images. These features indicate that a few crystal units of proto-dolomite may have been formed with some degree of cation ordering. The SED/EDS and HRTEM images also demonstrate that the bacteria are mineralized and form proto-dolomite crusts. The HTREM analyses of bacterial ultrathin slices show that the proto-dolomite crystals formed within cells, on the cell surface and within EPS. We assume that during biomineralization, Ca2+ and Mg2+ ions are trapped by the negatively charged EPS and transferred into the cells through diffusion, leading to the extracellular and intracellular precipitation of proto-dolomite; the bacteria themselves are mineralized during the constant formation of proto-dolomite crystals within and upon the surface of the cells. This study provides insights into our understanding of the formation of biotic proto-dolomite, and contributes to the origin of “microbial dolomite”.
{"title":"Mineralization of bacterial cells during proto-dolomite precipitation mediated by the halophilic bacteria Bacillus licheniformis","authors":"Ruirui Meng , Zuozhen Han , Wen Nie , Chao Han , Yu Han , Maurice E. Tucker","doi":"10.1016/j.chemgeo.2025.123205","DOIUrl":"10.1016/j.chemgeo.2025.123205","url":null,"abstract":"<div><div>Previous studies have demonstrated that bacteria can induce proto-dolomite precipitation at ambient temperatures. Yet, the mineralization of the bacteria during this process has rarely been reported. Here we describe mineralized bacteria formed during proto-dolomite precipitation mediated by the moderate halophilic bacteria <em>Bacillus licheniformis</em> Y1. The XRD data and SEM images show that proto-dolomite precipitated mainly as aggregates in media with Mg/Ca ratios of 6, 8 and 10 cultured at temperatures of 30 °C and 40 °C, but this was not observed in solutions cultured at 20 °C. Solution chemistry analysis showed that all solutions maintained a higher CO<sub>3</sub><sup>2−</sup> concentration and correspondingly higher saturation indices for proto-dolomite at 40 °C. The HRTEM results illustrate that the crystal structures and Mg/Ca atomic ratio of the proto-dolomite more closely match those of dolomite, and some superstructure reflections appear in the HRTEM and SAED images. These features indicate that a few crystal units of proto-dolomite may have been formed with some degree of cation ordering. The SED/EDS and HRTEM images also demonstrate that the bacteria are mineralized and form proto-dolomite crusts. The HTREM analyses of bacterial ultrathin slices show that the proto-dolomite crystals formed within cells, on the cell surface and within EPS. We assume that during biomineralization, Ca<sup>2+</sup> and Mg<sup>2+</sup> ions are trapped by the negatively charged EPS and transferred into the cells through diffusion, leading to the extracellular and intracellular precipitation of proto-dolomite; the bacteria themselves are mineralized during the constant formation of proto-dolomite crystals within and upon the surface of the cells. This study provides insights into our understanding of the formation of biotic proto-dolomite, and contributes to the origin of “microbial dolomite”.</div></div>","PeriodicalId":9847,"journal":{"name":"Chemical Geology","volume":"701 ","pages":"Article 123205"},"PeriodicalIF":3.6,"publicationDate":"2025-12-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145753515","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 : 2025-12-13DOI: 10.1016/j.chemgeo.2025.123201
Nadine P. Cooper , D. Graham Pearson , Andreas Auer , Tristan Bongartz , Yannick Bussweiler , Mario Fischer-Gödde , Max Hellers , Geoffrey Howarth , Andrew J. Locock , Ambre Luguet , Carsten Münker , Thomas Stachel , Ashlea N. Wainwright , Sarah Woodland , James M. Scott
The c. 86 Ma Jagersfontein kimberlite in the southern portion of the Kalahari craton in South Africa erupted a suite of refractory peridotite xenoliths with regionally anomalous rhenium+osmium (ReOs) isotope systematics and highly siderophile element (HSE) abundances. A portion of this suite comprises peridotites with 2.7 to 3.0 Ga Re- depletion 187Os/188Os model ages, refractory major elements and HSE abundances that represent some of the most depleted and least metasomatised known components of Archean cratonic nuclei. Other peridotites in the Jagersfontein suite show enrichment in Os, Ir, Pd, Pt and, especially, Ru, as well as having more radiogenic 187Os/188Os. Ruthenium abundances reach up to ∼22 times “Primitive Upper Mantle” levels and correspond to samples containing microscopic laurite plus pentlandite. The more radiogenic 187Os/188Os in HSE-enriched samples, along with the modern mantle like Ru isotope composition (ε1000.00 ± 0.12, ε102–0.03 ± 0.16), indicate that the prominent HSE enrichment is not related to an excess of unequilibrated Late Veneer material in the mantle source. Instead, the HSE characteristics and Os isotope ratios of these peridotites are inferred to have been caused by metasomatic addition of small volumes of HSE-rich sulfide/alloy to the refractory Archean peridotites. Although there is some variation in the metasomatic characteristics, the calculated 187Os/188Os of the contaminant causing the strongly Ru-enriched peridotites is ∼0.116, which is equivalent to the convecting mantle composition in the early Mesoproterozoic. Unlike the erosion and Proterozoic replacement of older mantle interpreted to explain post-Archean 187Os/188Os mantle beneath the central Kaapvaal nucleus, the Archean mantle lithosphere beneath Jagersfontein therefore appears to have been modified by the passage of Proterozoic melts. The 2400 km2 1.9 Ga layered mafic Trompsburg Complex, less than 10 km southeast of Jagersfontein, indicates that the mantle lithosphere in the area has experienced substantive transfer of large volumes of post-Archean melt, and we suggest that the passage of those melts may have caused the mantle HSE enrichment.
南非Kalahari克拉通南部的Ma Jagersfontein金伯利岩c. 86喷发出一套具有区域异常的铼锇(ReOs)同位素系统和高亲铁元素(HSE)丰度的难熔橄榄岩包体。该套的一部分包括2.7 ~ 3.0 Ga Re- depletion 187Os/188Os模式年龄的橄榄岩、难熔主元素和HSE丰度,代表了太古宙克拉通核中一些已知最贫化和最不交代的成分。Jagersfontein组的其他橄榄岩均富集Os、Ir、Pd、Pt,尤其是Ru,且具有较多的187Os/188Os放射性成因。钌丰度高达“原始上地幔”水平的22倍,对应于含有显微月辉石和镍黄铁矿的样品。富HSE样品的187Os/188Os具有较高的放射性成因,以及现代地幔中的Ru同位素组成(ε1000.00±0.12,ε102-0.03±0.16)表明,HSE的富集与地幔源中过量的不平衡晚期饰面物质无关。相反,推断这些橄榄岩的HSE特征和Os同位素比值是由于在难熔太古宙橄榄岩中交代添加了少量富HSE的硫化物/合金所致。虽然交代特征有一定的变化,但计算得到的强富钌橄榄岩污染物的187Os/188Os值为~ 0.116,相当于中元古代早期对流地幔成分。与解释Kaapvaal中央核下1870 / 1880年代后太古代地幔的侵蚀和旧地幔的元古代替换不同,Jagersfontein下的太古代地幔岩石圈似乎被元古代熔体的通过所改变。位于Jagersfontein东南不到10 km处的2400 km2 1.9 Ga层状基性特朗普斯堡杂岩表明,该地区地幔岩石圈经历了大量太古宙后熔融体的实质性转移,这些熔融体的通过可能导致了地幔HSE的富集。
{"title":"Proterozoic highly siderophile element enrichment of the Kaapvaal cratonic Archean mantle","authors":"Nadine P. Cooper , D. Graham Pearson , Andreas Auer , Tristan Bongartz , Yannick Bussweiler , Mario Fischer-Gödde , Max Hellers , Geoffrey Howarth , Andrew J. Locock , Ambre Luguet , Carsten Münker , Thomas Stachel , Ashlea N. Wainwright , Sarah Woodland , James M. Scott","doi":"10.1016/j.chemgeo.2025.123201","DOIUrl":"10.1016/j.chemgeo.2025.123201","url":null,"abstract":"<div><div>The c. 86 Ma Jagersfontein kimberlite in the southern portion of the Kalahari craton in South Africa erupted a suite of refractory peridotite xenoliths with regionally anomalous rhenium+osmium (Re<img>Os) isotope systematics and highly siderophile element (HSE) abundances. A portion of this suite comprises peridotites with 2.7 to 3.0 Ga Re- depletion <sup>187</sup>Os/<sup>188</sup>Os model ages, refractory major elements and HSE abundances that represent some of the most depleted and least metasomatised known components of Archean cratonic nuclei. Other peridotites in the Jagersfontein suite show enrichment in Os, Ir, Pd, Pt and, especially, Ru, as well as having more radiogenic <sup>187</sup>Os/<sup>188</sup>Os. Ruthenium abundances reach up to ∼22 times “Primitive Upper Mantle” levels and correspond to samples containing microscopic laurite plus pentlandite. The more radiogenic <sup>187</sup>Os/<sup>188</sup>Os in HSE-enriched samples, along with the modern mantle like Ru isotope composition (ε<sup>100</sup>0.00 ± 0.12, ε<sup>102</sup>–0.03 ± 0.16), indicate that the prominent HSE enrichment is not related to an excess of unequilibrated Late Veneer material in the mantle source. Instead, the HSE characteristics and Os isotope ratios of these peridotites are inferred to have been caused by metasomatic addition of small volumes of HSE-rich sulfide/alloy to the refractory Archean peridotites. Although there is some variation in the metasomatic characteristics, the calculated <sup>187</sup>Os/<sup>188</sup>Os of the contaminant causing the strongly Ru-enriched peridotites is ∼0.116, which is equivalent to the convecting mantle composition in the early Mesoproterozoic. Unlike the erosion and Proterozoic replacement of older mantle interpreted to explain post-Archean <sup>187</sup>Os/<sup>188</sup>Os mantle beneath the central Kaapvaal nucleus, the Archean mantle lithosphere beneath Jagersfontein therefore appears to have been modified by the passage of Proterozoic melts. The 2400 km<sup>2</sup> 1.9 Ga layered mafic Trompsburg Complex, less than 10 km southeast of Jagersfontein, indicates that the mantle lithosphere in the area has experienced substantive transfer of large volumes of post-Archean melt, and we suggest that the passage of those melts may have caused the mantle HSE enrichment.</div></div>","PeriodicalId":9847,"journal":{"name":"Chemical Geology","volume":"702 ","pages":"Article 123201"},"PeriodicalIF":3.6,"publicationDate":"2025-12-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145731849","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 : 2025-12-12DOI: 10.1016/j.chemgeo.2025.123204
J. Ahmadi , E. Widom , M.-N. Guilbaud , C. Siebe , D.C. Kuentz , R. Sanchez , M. Medina-Jaen , J.J. Zamorano-Orozco , S.M. Straub
Xitle Volcano, an ∼1700-year-old scoria cone located at the southern edge of Mexico City and within the central Trans-Mexican Volcanic Belt (TMVB), provides a unique opportunity to investigate the causes of geochemical diversity in monogenetic eruptions. We present the first integrated major-trace element and Sr–Nd–Pb–Os isotopic dataset which, coupled with new stratigraphic mapping and petrologic observations of Xitle's eruptive products, provides constraints on the roles of mantle heterogeneity, subduction-related metasomatism, and crustal assimilation in Xitle magma genesis. The lavas and tephras range in composition from basalt to basaltic andesite and basaltic andesite/basaltic trachyandesite. All samples have arc-like trace element patterns with enrichments in large ion lithophile elements and depletions in high field strength elements, with early-stage magmas exhibiting stronger arc-like signatures, while later magmas trend toward intraplate affinities. Isotopic compositions span 87Sr/86Sr = 0.70335–0.70398, 143Nd/144Nd = 0.512790–0.512953, 206Pb/204Pb = 18.667–18.770, 207Pb/204Pb 15.570–15.612, and 208Pb/204Pb = 38.355–38.554. 187Os/188Os ratios range from 0.1368 to 0.1968, with Os concentrations between 10.7 and 168.4 pg/g. These 187Os/188Os values are more radiogenic than the primitive upper mantle (∼0.1296) and fall within the range observed for other arc magmas in the TMVB. The samples collectively define three distinct compositional groups, but the compositional variations do not evolve monotonically with time; rather abrupt compositional and isotopic shifts occur over the course of the eruption. Neither closed system fractional crystallization nor assimilation-fractional crystallization of a single parental magma can reproduce the observed geochemical range. Instead, we propose an open-system scenario involving high-temperature fractional crystallization periodically interrupted by recharge events from compositionally and isotopically distinct magmas, supported by MELTS and Magma Chamber Simulator (MCS) modeling. Isotope systematics argue against significant crustal assimilation and point to a heterogeneous, subduction-metasomatized mantle wedge as the primary source of the magmas from Xitle. The mantle source beneath Xitle was variably fluxed by two distinct slab-derived components, both dominated by terrigenous sediment with variable inputs of altered oceanic crust (AOC)-derived fluid/melts, in relative proportions that evolved over time. Our results demonstrate that Xitle's chemical complexity reflects serial addition of subduction fluids, partial melting of a compositionally evolving mantle source, and repeated magma recharge, emphasizing that even small monogenetic systems can record prolonged and multi-stage magmatic histories of mantle derived melts.
{"title":"Deciphering chemical complexity in monogenetic eruptions: Xitle Volcano, Mexico","authors":"J. Ahmadi , E. Widom , M.-N. Guilbaud , C. Siebe , D.C. Kuentz , R. Sanchez , M. Medina-Jaen , J.J. Zamorano-Orozco , S.M. Straub","doi":"10.1016/j.chemgeo.2025.123204","DOIUrl":"10.1016/j.chemgeo.2025.123204","url":null,"abstract":"<div><div>Xitle Volcano, an ∼1700-year-old scoria cone located at the southern edge of Mexico City and within the central Trans-Mexican Volcanic Belt (TMVB), provides a unique opportunity to investigate the causes of geochemical diversity in monogenetic eruptions. We present the first integrated major-trace element and Sr–Nd–Pb–Os isotopic dataset which, coupled with new stratigraphic mapping and petrologic observations of Xitle's eruptive products, provides constraints on the roles of mantle heterogeneity, subduction-related metasomatism, and crustal assimilation in Xitle magma genesis. The lavas and tephras range in composition from basalt to basaltic andesite and basaltic andesite/basaltic trachyandesite. All samples have arc-like trace element patterns with enrichments in large ion lithophile elements and depletions in high field strength elements, with early-stage magmas exhibiting stronger arc-like signatures, while later magmas trend toward intraplate affinities. Isotopic compositions span <sup>87</sup>Sr/<sup>86</sup>Sr = 0.70335–0.70398, <sup>143</sup>Nd/<sup>144</sup>Nd = 0.512790–0.512953, <sup>206</sup>Pb/<sup>204</sup>Pb = 18.667–18.770, <sup>207</sup>Pb/<sup>204</sup>Pb 15.570–15.612, and <sup>208</sup>Pb/<sup>204</sup>Pb = 38.355–38.554. <sup>187</sup>Os/<sup>188</sup>Os ratios range from 0.1368 to 0.1968, with Os concentrations between 10.7 and 168.4 pg/g. These <sup>187</sup>Os/<sup>188</sup>Os values are more radiogenic than the primitive upper mantle (∼0.1296) and fall within the range observed for other arc magmas in the TMVB. The samples collectively define three distinct compositional groups, but the compositional variations do not evolve monotonically with time; rather abrupt compositional and isotopic shifts occur over the course of the eruption. Neither closed system fractional crystallization nor assimilation-fractional crystallization of a single parental magma can reproduce the observed geochemical range. Instead, we propose an open-system scenario involving high-temperature fractional crystallization periodically interrupted by recharge events from compositionally and isotopically distinct magmas, supported by MELTS and Magma Chamber Simulator (MCS) modeling. Isotope systematics argue against significant crustal assimilation and point to a heterogeneous, subduction-metasomatized mantle wedge as the primary source of the magmas from Xitle. The mantle source beneath Xitle was variably fluxed by two distinct slab-derived components, both dominated by terrigenous sediment with variable inputs of altered oceanic crust (AOC)-derived fluid/melts, in relative proportions that evolved over time. Our results demonstrate that Xitle's chemical complexity reflects serial addition of subduction fluids, partial melting of a compositionally evolving mantle source, and repeated magma recharge, emphasizing that even small monogenetic systems can record prolonged and multi-stage magmatic histories of mantle derived melts.</div></div>","PeriodicalId":9847,"journal":{"name":"Chemical Geology","volume":"703 ","pages":"Article 123204"},"PeriodicalIF":3.6,"publicationDate":"2025-12-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145731851","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}