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

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

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