プチスポット火山群を通した水循環・変質作用による海洋プレート改変の数値モデリング

2022 年度 実施状況報告書

• 研究課題/領域番号: 22K03765
• 研究機関: 国立研究開発法人海洋研究開発機構
• 研究代表者: 川田 佳史 国立研究開発法人海洋研究開発機構, 海洋機能利用部門(海底資源センター), 臨時研究補助員 (50402558)
• 研究期間 (年度): 2022-04-01 – 2025-03-31
• キーワード: 海底熱水循環 / 太平洋プレート / 地殻熱流量 / 数値計算法 / 有限要素法 / 有限差分法 / 陽解法

研究実績の概要

太平洋プレートのような古いプレートが沈み込む直前には、プチスポットと呼ばれる海底火山が新たに顔を出すことがある。太平洋プレート上のプチスポット近傍では高い地殻熱流量が観測されており、海底を通した水循環が起こってると考えられている。本研究課題は、このプチスポット海底火山の中を流れる熱水循環が、海底地殻をどのように改変するかを調べる研究である。本研究により、例えば沈み込むプレートの温度構造がどのように変わるかを議論することができ、その結果として、沈み込み帯の巨大地震に温度構造が与える影響を見積もることができるようになる。
プチスポットの中を流れる熱水循環を計算するための数値計算法についての検討を行った。本研究では実装が簡単な加速化陽解法を適用することを計画している。陽解法で重要な点は、許される時間ステップを知ることである。そこで、まず通常の流体についての加速化陽解法について、時間ステップが決まる仕組みを明らかにした。本研究では山体の変形に関連する要素である。
次に、浸透流の熱対流計算を陽解法を用いて解くための手法を開発した。本研究が扱うのは山体の中の水循環であるから、計算メッシュは不規則にならざるを得ない。そこで、有限要素法を用いることが必須である。そこで加速化された陽解法の有限要素法への実装を行った。しかし陽解法を不規則メッシュの系にそのまま適用すると実は効率が悪いことがわかっている。一方このような系を差分化された式（計算機が扱う形式）のレベルで見ると、計算メッシュが均一で浸透率が不規則な場合に似ている。そこで、浸透率が不規則な場合の計算法を計算メッシュが不規則な場合に適用した。その結果、この方法が不規則メッシュの系を効率よく計算できることを確かめることができた。以上より、本研究課題が扱う予定の数値計算方法を実装することができた。

現在までの達成度 (区分)

3: やや遅れている

理由

ここまで行ってきた数値計算法の検討により、今後の数値計算を効率よく行うためのいくつかの方法を見い出すことが出来た。具体的には、加速化陽解法の時間ステップが決まるメカニズムを明らかにすること、有限要素法に加速化陽解法を組み込むこと、不規則メッシュ系について陽解法を効率的に適用すること、等が実現できた。以上により、本研究課題の遂行に必須である、プチスポットのような複雑な形状をもつ山体を流れる熱水循環の計算を効率的に行える方法論が完成した。同様に、山体の変形についての計算についても効率よい計算が行うことができる。他方、当初計画していた、山体の大きさなどをパラメータをとした組織的な計算はまだ完全には行えていない状況である。以上から、研究はやや遅れていると判断した。

今後の研究の推進方策

今後は、プチスポットの中を流れる熱水循環について、山体の大きさや透水性、複数の山体を考慮する場合は密集度をパラメータとして、組織的な数値計算を行うことを計画している。計算機環境の構築は済んでいるため、計算の実行は十分可能であると考えている。ここまでの計算ができれば、研究開始2年度としてはほぼ当初の計画通りということになる。
また、熱水変質により空隙が詰まる過程に対する数値計算を開始したい。ただし、具体的な反応を組み込む前に、研究代表者が過去の研究において行ったように、概念的な反応を用いることで空隙が詰まる過程を表現することを検討している。熱水変質の詳細（反応速度などのパ熱力学的パラメータ）については今後文献調査を行う予定であり、これを用いた計算を行うのは、最終年度になると考えている。