日本屈斜路（Kussharo）火山口下方岩板（Slab）来源岩浆水（Magmatic Water）的氢氧同位素证据

《Journal of Volcanology and Geothermal Research》：Isotopic evidence for slab-derived magmatic water beneath the Kussharo Caldera, Japan

　　

众所周知，随俯冲海洋板块（岩板，Slab）携带的水在岛弧岩浆形成中起关键作用，但直接证明岩浆水源自岩板的证据仍有限。为此，研究人员对日本北海道屈斜路（Kussharo）火山口内及周边的温泉水、冷泉水、河水、喷气孔蒸气及热水液（Hydrothermal Water）进行了同位素和溶解离子分析。该研究区太平洋板块于千岛海沟俯冲至北美板块之下，十余座第四纪火山沿与海沟平行的火山前锋（Volcanic Front）分布。结果显示，位于火山口内活熔岩穹丘附近的川汤（Kawayu）温泉并非单纯源于大气降水（Meteoric Water）；δ2H、δ18O、Cl?及SO42?浓度间呈强线性相关性，指示其为大气降水与岩浆水的混合产物，且岩浆端元的同位素组成与该处模拟的太平洋岩板来源水（Pacific Slab-derived Water）数值吻合，为岩浆水源自岩板提供了直接证据。相比之下，喷气孔蒸气与热水液的同位素组成因受汽—液分离（Vapor–Liquid Separation）伴随的同位素分馏（Isotopic Fractionation）影响而发生部分偏移。本研究强调，尽管围岩同位素组成可能因同位素交换（Isotopic Exchange）显著改变，大量岩板来源水在上升过程中仍可保留其原始同位素特征。

《Journal of Volcanology and Geothermal Research》刊载论文解读——日本Kussharo火山口下方岩板来源岩浆水的同位素证据

一、研究背景与意义

岛弧岩浆成因的经典模型认为，俯冲海洋板块（Slab）在深部经脱水反应（Dehydration Reaction）释放富水流体进入地幔楔（Mantle Wedge），降低地幔固相线从而诱发部分熔融形成岛弧岩浆。以往研究多依靠熔岩中B、Sr、Nd、Pb等同位素示踪判定岩板组分加入，但直接约束岩浆水中"水本身"来源于俯冲岩板（Slab-derived Water）的证据极为匮乏。此外，前人（如Giggenbach, 1992）提出的安山质岩浆水（Andesitic Magma Water）典型值（δ²H = ?20±10‰，δ¹⁸O = +8±2‰）未考虑两点局限：一是热液系统中汽—液分离造成的同位素分馏会改造原始信号；二是岩板水同位素组成随俯冲深度、温度及渐进脱水发生空间变化，非常量。Adachi and Yamanaka (2024)已建立考虑氢氧同位素交换与渐进脱水的岩板水同位素演化数值模型，并在日本非火山前弧区成功识别菲律宾海板块水与太平洋板块水，但该模型能否应用于火山弧区地表出露流体尚未验证。本研究选取北海道东侧屈斜路（Kussharo）火山口——太平洋板块于此经千岛海沟俯冲、火山前锋发育、具活跃熔岩穹丘（Atosanupuri）与温泉群——通过地表流体氢氧稳定同位素与地球化学分析，结合数值模拟，检验岩浆水端元同位素组成是否与当地125 km深处太平洋岩板脱水释放水模拟值吻合，以直接证实岩浆水之岩板来源，并探讨其在上升过程中抗同位素改造的机制。

二、主要关键技术方法

研究人员于2017–2024年采集屈斜路火山口川汤温泉（H1–H3，n=11次）、周边12处温泉、2处冷泉、8条河水，以及Atosanupuri活火山穹丘喷气孔的热水液（Fw）与冷凝蒸气（Fv）。氢氧稳定同位素（δ²H、δ¹⁸O）用Picarro L2130-i腔衰荡光谱仪测定；主要阴离子（Cl^?、SO₄^2?等）用离子色谱测定。采用两种途径确定大气降水—岩浆水混合线（Mixing Line, MXL）：①基于川汤温泉多年重复同位素数据的时空追踪回归；②引入流域尺度降水同位素景观模型（Isoscape Model）估算平均补给大气水同位素值，二者联合回归以降低外推误差。岩浆水端元（Magmatic End-member）取MXL与海源岩石圈水曲线（Ocean-origin Lithospheric Water Curve, OLWC）交点求解。对喷气孔样品作Rayleigh分馏校正（平均池温91.6℃）以还原原始流体同位素。用Adachi and Yamanaka (2024)岩板水同位素演化模型做蒙特卡洛模拟（1000次），输入太平洋板块俯冲温—深关系，输出125 km深度（屈斜路下方岩板估算深度）处岩板水δ²H、δ¹⁸O并与混合分析重建的岩浆水端元对比。

三、研究结果

3.1. 温泉、冷泉与河水样品（Hot/cold springs and river samples）

河水和冷泉水δ²H–δ¹⁸O投影落于当地大气水线（Local Meteoric Water Line, LMWL: δ²H = 7.52δ¹⁸O + 5.75）上；而川汤温泉（H1–H3）明显偏右（高δ¹⁸O侧）。川汤温泉δ²H、δ¹⁸O与Cl^?（r²=0.972–0.998）、SO₄^2?（r²=0.972–0.998）浓度呈强线性正相关，排除水—岩氧同位素交换（不伴阴离子协同变化）及蒸发富集主导可能，支持为低δ低离子大气降水与高δ高离子岩浆水两元混合。据此得混合线MXL：δ²H = 2.63δ¹⁸O ? 46.26（r²=0.988, n=12）。湖畔H4温泉无此偏移，基本为大气降水来源。距火山弧外缘知床硫黄山（Mt. Shiretoko-Iozan）下H8温泉同样右偏且有高Cl^?，暗示也存在岩浆水混入；H11高Cl^?低δ偏移归为演化化石海水而非岩浆水；其余外围温泉均属大气降水来源。

3.2. 喷气孔样品（Fumarolic samples）

喷气孔热水液与蒸气δ值较川汤温泉更偏离LMWL，其中残余热水液偏右上、蒸气偏左下，二者差值可用91.6℃平衡同位素分馏（汽—液分离）解释，表明原始流体本位于MXL上，经开放系统沸腾蒸发导致轻同位素优先进入气相而产生分馏偏移，故喷气孔样品不宜直接用于岩浆水端元反演。

3.3. 端元混合分析结果（End-member mixing analysis）

MXL与LMWL交点得大气水端元δ²H = ?74.23‰，δ¹⁸O = ?10.64‰，与流域isoscapes估算均值吻合。MXL与OLWC（δ¹⁸O_fin=11±1‰，δ²H_off=0±10‰）交点得岩浆水端元同位素组成δ²H_ma= ?23.93‰，δ¹⁸O_ma= +8.49‰，落入经典安山质岩浆水范围。川汤温泉中岩浆水占比13.4%–16.7%，Atosanupuri喷气孔热水液换算占比31.4%–43.4%。由δ–Cl^?/SO₄^2?混合关系得岩浆水端元Cl^?= 9234(±74) mg L^?1，SO₄^2?= 17525(±173) mg L^?1，具中等Cl、极高分SO₄特征，区别于前弧非火山热泉（高Cl低SO₄）。

3.4. 与数值模拟结果对比（Comparison with numerical simulation）

岩板水同位素演化模型对太平洋板块俯冲过程（升温+渐进脱水）做1000次蒙特卡洛模拟，125 km深度（屈斜路下方岩板深度）中位数预测值δ²H = ?23.61‰，δ¹⁸O = +8.02‰，不确定区间与混合分析重建岩浆水端元（δ²H = ?23.93±~0.6‰，δ¹⁸O = +8.49±~0.3‰）在误差范围内重合，表明岩浆水同位素特征可完全由该处太平洋岩板脱水释放水解释。

四、讨论与结论翻译

研究人员进一步将Kussharo所得岩浆水端元代入北海道西侧Usu火山Showa-shinzan熔岩穹丘历史喷气孔数据（Mizutani, 1978；同岩板深度近似），发现原解释为"岩浆水+¹⁸O偏移大气水"可重新解释为仅混合+汽—液分馏效应，无需引入氧漂移，支持本文重建值的普适性。尽管岩板来源水上涌经地幔楔和岩浆时理论上可发生同位素再平衡（与地幔橄榄岩平衡δ²H≈?60±20‰），但实际未观测到此偏移，很可能是因为大量岩板来源水具高累积水/岩比（Cumulative Water/Rock Ratio >0.05对δ²H，>1对δ¹⁸O，700℃开放体系），流体改写了围岩同位素而非被围岩改写，从而保留了原始同位素信号。按川汤温泉流量估算，屈斜路活动期岩板来源水释放总量下限约94.6 km³。

结论（翻译）：

位于屈斜路火山口Atosanupuri熔岩穹丘附近的川汤温泉含岩浆组分，其特征为显著偏高δ¹⁸O（+8.49‰）、中等偏高δ²H（?23.93‰）、极高SO₄^2?（17525 mg L^?1）及中等偏高Cl^?（9234 mg L^?1）。重建之岩浆水同位素组成与当地125 km深处太平洋岩板脱水释放水模拟值高度吻合；此一致性亦见于北海道西侧Showa-shinzan熔岩穹丘喷气孔历史数据。结果表明岩浆水确系岩板来源，且岩板来源水同位素组成在上涌过程中不易受与地幔橄榄岩或硅酸盐熔体相互作用的次生改造；累积岩板来源水排放量相对岩浆及水—岩交互围岩体积较大，得以保存原始信号。虽结论基于北海道北部，其适用性尚需在其他板块构造背景与岛弧中检验。