火山性水銀ガスの熱力学‑動力学的スペシエーション解析

概要

火山から放出される水銀は、単なるフラックスだけでは環境挙動を予測できません。大気中への輸送・沈着は、水銀がどの化学形態（元素状、ハロゲン化物、硫化物など）で存在するかに強く依存します。本稿では、主要ガス系（H‑O‑C‑S‑Cl‑Br）を熱力学的背景として扱い、その上で微量水銀を各リザーバーへ分配する手法を提示しています。

背景と課題

従来、火山性水銀は「モル数」や「キログラム／時間」といったフラックスで報告されてきましたが、同じモル数でも元素状水銀とハロゲン化物・硫化物では環境影響が大きく異なります。そのため、フラックスだけでは化学的な完全性が欠け、輸送・沈着・保存過程との橋渡しが不十分でした。

熱力学的アプローチ

本研究では、主要ガス系の圧力・温度・組成を基に熱平衡図を構築し、そこから水銀がどのリザーバーに優先的に入るかを決定します。ハロゲン供給量や硫黄活性度が変化すると、ハロゲン化物が支配的になる領域と硫化物が支配的になる領域が明確に分かれます。

動力学的クエンチ評価

実際の火山プラズマは急速に冷却・希釈され、反応速度が平衡到達より遅くなることが多いです。この「クエンチ」を、平衡曲線からの逸脱度合い（対数活量比空間での遅れ）として定量化し、リザーバー切替の指標としました。

エトナ型ベンチマーク

簡易数値モデルを用いてエトナ火山のケースを再現。高温のマグマガスでは元素状水銀が優勢ですが、冷却・酸化・HCl濃度上昇に伴いハロゲン化物が増加する傾向が確認されました。これにより、熱力学的予測と実測データの整合性が示されています。

ワークフローと結論

1. ガス組成と温度・圧力から熱平衡境界を算出。
2. 冷却・混合速度を考慮し、クエンチ指数を導入。
3. 得られた境界とクエンチ指標を組み合わせ、火山ごとのスペシエーション図を生成。
この手順は汎用的であり、個別の火山条件に合わせて数値入力を変更すれば、再現可能な平衡・クエンチ図が得られます。今後は、より詳細な熱化学データベースや実測観測と統合し、モデル精度の向上が期待されます。

補足資料

付録には、使用した熱化学定数表、数式の導出過程、および追加診断結果が掲載されています。

元記事: https://www.frontiersin.org/journals/earth-science/articles/10.3389/feart.2026.1902534/full