(1)噴火のダイナミクスの解明を目指した実験と理論研究
マグマ破砕過程を「粘弾性流体の破壊現象」と位置づけ,定量的モデル化に向けた粘弾性構成方程式の構築と数値計算手法の開発を進めた.単純なマクスウェル型の粘弾性を示す光弾性物質を用いた変形・破壊実験に着手し,加速を伴う3次元の変形場の中で,流動から破壊へと遷移する様子を,光弾性を利用した弾性歪の可視化を含めて観察した.また,気泡の膨張に伴う流体の破壊と流動挙動について,この粘弾性流体と降伏強度より小さい応力下で弾性を示すジェル状流体の比較を行った.気泡への気体供給速度を増加すると,粘弾性流体は脆性破壊を生じたが,ジェル状流体は流動速度が増加した.これらの流体のレオロジーと,水蒸気噴火の噴出物である火山泥のレオロジーを比較したところ,火山泥はジェル状流体に近い性質を持っていることが分かった.マグマにおいても,結晶を含む低粘性マグマはジェル状流体に近いレオロジーが報告されている.これらの結果は,これまで弾性の存在と歪み速度のみに基づいて判定されていた火山物質の脆性破砕条件を見直す必要のあることを示唆している.
これまで開発してきた気泡変形モデルと噴火の火道流モデルを組み合わせ,爆発的噴火噴出物の気泡変形度頻度分布と,噴火時の歪み速度プロファイルとの定量的関係づけを行った,この理論を,最も新しいカルデラ噴火であるタウポ火山の1800年前の噴火に適用し,プリニー式噴火時と大規模火砕流噴火時の噴出物気泡組織の違いが火道径の違いによって生じうることを示した.
(2)火山噴煙ダイナミクスのシミュレーション研究
爆発的火山噴火で見られる噴煙柱・火砕流の噴煙ダイナミクスと,火山灰輸送・堆積プロセスの解明を目指し,数値モデルの開発とそれを用いた大規模シミュレーション研究を進めている.火山灰は噴煙によって上空へと運ばれ,噴煙から離脱すると大気風によって広範囲に移流・拡散する.そのため,噴煙柱の高度・傘型噴煙の拡大速度・大気風による噴煙の変形度合い・火山灰粒子の噴煙からの離脱位置が,最終的な地表での降灰分布を決定する.この問題に取り組むため,噴出物と大気からなる2流体とトレーサー粒子をカップリングした数値モデルを用い,複数のスーパーコンピュータにおいて大規模シミュレーションを実施した.噴火強度が中程度と大きい場合について,無風状態の大気を仮定した場合のパラメータスタディを行った.その結果,火山灰の噴煙からの離脱高度が噴火強度によって定性的に変化する様子を捉えることができた.また,それぞれの噴火強度の火口条件について,大気の風速を変化させた場合のパラメータスタディを行った.計算の結果,火山灰は粒子サイズが小さくなるにつれ,噴煙柱の途中・噴煙柱上部・傘型噴煙と異なる地点から離脱することが分かった.噴火強度が大きい場合には細かな粒子は同心円状に広がる傘型噴煙によって運搬されるが,いったん傘型噴煙から離脱すると大気風の影響を強く受けて風下側に堆積した.これらの火山灰粒子挙動の解析は,噴火想定図やハザードマップ作成の基礎データにつながる.噴煙の流体計算結果に基づき,傘型噴煙や火砕流の拡大スケーリングについても解析を行った.
(3)大規模噴火に関する研究
南九州鬼界カルデラの活動履歴や7.3 ka鬼界アカホヤ噴火の推移を解明するための研究を進めている.とくに従来アカホヤ噴火の前駆的活動により形成された可能性が指摘されていた長浜溶岩(流紋岩質溶岩)やそれ以前の活動の実態を明らかにするために,2018年にボーリング掘削を実施した,それにより得られた試料の解析を進めた結果,長浜溶岩は深度11-190 m(水深130 mに相当)に存在し,その直下の深度190-230 mには貝殻を含む粗粒砂質層や複数枚のテフラ層を主体とした海成の地層が存在することがわかった.さらに下位(230 m以深)には斜長石斑晶に富む複数枚の安山岩質溶岩が存在する.長浜溶岩直下の砂層に含まれる複数の貝殻の14C年代測定を行ったところ,7000〜8300 calBPの年代値が得られた.これにより,長浜溶岩の活動が鬼界アカホヤ噴火に先行する活動であったことがはじめて地質学的・年代学的に明らかになった.長浜溶岩およびその下位のテフラ層(12-15 ka)とアカホヤ噴火の岩石学的関係,大規模噴火に先行する溶岩流活動の役割など,巨大噴火を起こしたマグマシステムとその進化について研究を進めている.
(1)火山の空振モニタリング手法の開発
火山噴火に伴う空振の波形や振幅を正確に計測するため,新しい空振計を開発している企業や工学系の研究者らと協力し,小型・低消費電力マイクロフォンやMEMSセンサー,高精度気圧計の比較試験および火山地域における長期評価試験を行い,必要な改良を進めた.
より効率のよい空振アレイ観測の方法として,従来のアレイ観測よりも一桁空間スケールの小さい,10メートルサイズの極小規模アレイの開発を行い,さらに,地上2~4m程度の高さに1要素加えることによって,方位角だけでなく仰角の分解能が向上させられることを示した.イタリアのストロンボリ火山において,極小規模空振アレイを用いて観測したデータの詳細な解析を行った.そして,アレイ解析と波動場数値計算を組み合わせることにより,これまでにない精度で空振源の高さに制約を与えることができた.また,活発な噴気活動を続ける霧島火山硫黄山周辺において,大小二つの噴気からの空振シグナルを計測・分離する観測実験を行った.3要素極小規模空振アレイと単独空振計による観測データの統合解析により,音源分解に成功し,両者の振幅比が15:1であると決められた.強度に差のある空振場で微弱な音源を把握することは容易ではなく,本研究の手法は,噴気地帯や複数の活動的火口を有する火山の監視や観測に役立つものと考えている.さらに,同じ極小規模空振アレイを用いて,霧島新燃岳の西側の山域で低周波音源探査を実施した.アレイの移動と設置を繰り返し,2時間以内に4地点で各15分程度の短時間観測を行い,噴気地帯から低周波音波が発生していること見出した.このような逐次移動観測は極小規模空振アレイの新たな応用法である.
(2)無人ヘリやドローンを活用した火口近傍観測システムの開発と応用
活動的な火山において,観測者を危険にさらすことなく火口周辺での様々な観測を実施することを目的として,2008年から無人ヘリを用いた火口近傍観測システムの開発を進めている.産業用無人ヘリを火山観測に利用するため,様々な火山での飛行実績を積むとともに,観測に必要な様々な周辺機器,静止画・動画撮影用の機器を搭載するための専用雲台,地震計やGPS観測装置をヘリから降下設置するウインチ,無人ヘリ設置用の地震計モジュール,GPSモジュールなどの開発を進めてきた.口之永良部島では2015年4月に火口近傍の4箇所に地震計を設置した.この地震計は2015年5月の噴火で失われたが2015年9月に再度5点を設置した.観測データから2015年5月29日の噴火に先行して火口近傍で地震が急増していたこと,単色地震も増加していたことがわかった.また,可視画像・熱映像・電磁気・ガス等の多項目データから,活動の大きな変化も捉えられた.火口に接近して得られたガスの分析により脱ガス時の見かけ平衡温度も推定された.2016年6月には,火口から1.5km内が警戒範囲となっている西之島において,気象庁と共同で無人ヘリ(船上より離発着および制御)により活動・噴出物の観察および岩石試料の採取を行った.また,2009年から2017年にかけて,桜島山頂付近に地震計およびGPS受信機を設置した.桜島山頂の地震計は2021年2月末時点も稼働を続けている.
無人ヘリによる空中磁気測量も精力的に行っている.2011年霧島新燃岳噴火後の山体の帯磁状態の変化を把握するため,2011年5月,11月,2013年11月,2014年10月,2015年11月,2017年11月,2018年11月の計7回,新燃岳およびその西側,およそ3㎞四方の領域において,繰り返し空中磁気測量を実施した.測線間隔および対地高度はおおよそ100mで一定として測定フライトを実施した.プログラムした航路に沿って正確に測定飛行できることは繰り返し測量にとって大きな利点である.解析の結果,新燃岳火口内の溶岩は平均として4.0 A/m帯磁したと想定すると観測された全磁力データをよく説明することが判り,火口に蓄積された溶岩が熱拡散過程で順調に冷却している様子を明確にとらえることに成功した.また,三宅島においては,今後の火山活動を把握するための基礎資料とするために無人ヘリを用いた詳細な空中磁気測量を2014年5月と2016年11月に実施し,2017年度に解析を進めた結果,山体北側で負,南側で正の変化を検出した.その後,2019年6月にも実施している.2018年1月 に噴火した草津白根山・本白根山においても無人ヘリによる空中磁気測量を実施し,過去有人機により得られたデータとの比較解析を進めている.
電動モーターを動力源とするいわゆる「ドローン」の性能が近年大幅に向上し,火山観測において活用できるレベルに達しつつある.火山センターではドローンを活用した火山観測も進めている.新燃岳においては,ドローンによる火口内への接近撮影を実施し,西之島においては船上から飛ばしたドローンによる画像撮影と試料採取を実施した.霧島山・硫黄山ではドローンによる繰り返し空中磁気測量の活用実験を開始し,2019年に複数回の測定を実施した.その結果,無人ヘリよりも低廉かつ機動的に観測を実施できることが確認できた.これを受けて,三宅島においてもドローンによる空中磁気測量を試みている.無人ヘリは広域をカバーする測量に適しているが,経費や機動性にやや問題がある.一方,電動ドローンの飛行性能は年々向上しており,観測対象によっては無人ヘリに置き換える観測手段となり得る.今後は,観測対象に応じて両者を使い分けることになろう.
(3)衛星技術を活用した火山活動の把握
ひまわり8号とJAXAのしきさい(GCOM-C/SGLI)の赤外画像を用いてアジア太平洋域の主要活火山のリアルタイム観測を行うと共に,これを基盤データとし各種高分解能画像・現地観測データ等を組合せ,噴火推移・噴火プロセス解明に関する研究を進めている.この一環として,高頻度観測が可能なひまわり8号の熱異常データを用いた噴出率推定方法の開発を行った.ひまわり8号の1.6-㎛画像での熱異常と噴出率の関係を検討し,両者の間に高い相関関係があり,この回帰式が,Y = 0.47 X(Y:噴出率 106 m3 day-1 ,X:輝度値 106 W m-2 sr-1 m-1 )と求められることを示した.この式を用いて,西之島2019-2020年噴火初期の噴出率を推定し,この時の噴出率は2013-2015年西之島噴火における平均噴出率より2-3倍高かったことを明らかにした.一方,しきさいのSGLI画像は分解能が250 mと比較的高く,溶岩流の拡大や火砕流の発生等,噴火状況の変化を高頻度で捉えることができる.このSGLI画像により,カムチャッカ半島に位置するシベルチ火山の2019年噴火の観測を行い,溶岩ドームの成長率が低い状態から急上昇する時,火砕流が発生していることを見出した.また,ジャワ島東端部にあるイジェン火山の火口湖の観測を行い,2019年5月中旬から6月にかけて湖水温が最高38℃まで上昇し,火山下のマグマあるいは熱水活動が活発化した可能性を示した.
火山センターでは,火山やその深部で進行する現象の素過程や基本原理を解き明かし火山噴火予知の基礎を築くことを目指して,火山や噴火に関連した諸現象の研究を行っている.その基本的な研究方針は地震研究所の2009年サイエンスプランで掲げられた「火山活動の統合的解明と噴火予測」および科学技術学術審議会による2019年1月に出された「災害の軽減に貢献するための地震火山観測研究計画(第2次)の推進について(建議)」に基づいている.
本センターは2004年度に「火山観測の将来構想」を作成し,その中で主たる観測対象とする火山を3つに分類し,a)観測網を強化し研究成果を上げるべき火山として,浅間山,伊豆大島の2火山を,b)研究成果が短期的には大きく望めないが,将来のために観測を継続・改良すべき火山として,三宅島・富士山・霧島山の3火山を,c)他機関が既に観測網を整備している等の理由で基本的には撤退する火山として草津白根山を挙げた.全国の火山噴火予知研究コミュニティーで了解を得つつ,この構想に基づいて順次更新・整備を進めてきた結果,浅間山・伊豆大島では多項目観測網の強化が進み,霧島山では広帯域地震観測網を火山体に集中することができた.富士山の観測は着実に継続しており,三宅島については次の噴火に備えた観測網の強化が進む.また草津白根山については撤退を完了した.2010年度以降は,観測所等の施設は観測開発基盤センターに移管されたが,同センターの火山担当教員との協力・共同の元に研究方針に沿った整備を進めている.
観測網の強化・整備は一段落したと言えるが,観測を担うことのできる人材が急速に減少しつつあり,現状の観測網をこれまでどおりに維持することは相当な困難を伴う.これは地震研に限らず全国的な傾向であり,国内の火山観測における大学の役割を見直す時期が来ていると考えられる.伊豆大島や三宅島等次の噴火が切迫する火山を念頭に置きつつ,新たな火山観測研究の将来構想を早急にまとめる必要がある.
本センターの観測研究の対象となっている主たる火山の近年の活動は以下のとおりである.浅間山では2004年の活発な活動以降に大きな活動は無く,2008年,2009年,2015,2019年に弱い噴火が発生したのみである.しかしながら,2019年8月に発生した小噴火は活動度が極めて低い状態で発生したものであり,気象庁の噴火警戒レベルは1であった.この不意打ちともいえる噴火は観測を継続する上で安全管理の問題に影響を及ぼしている.伊豆大島は顕著な活動は無いものの,マグマ蓄積を示す基線長の伸びは継続している.富士山では目立った活動は無い.霧島山・新燃岳では2011年1月26日に約300年ぶりの本格的な準プリニー式噴火が発生し,それ以降は活発な活動が継続している.新燃岳は2017年10月に再噴火し,翌2018年3月には山頂火口から溶岩が溢れて北西に流れ下った.爆発的な噴火活動は2018年6月まで続き,山頂火口を埋めた溶岩や西斜面の噴気孔からは今も噴気が立ち上る.霧島山・硫黄山では2018年4月に小規模な水蒸気噴火が発生し,噴気活動は消長を繰り返しつつ2021年中も続いている.口永良部島では2015年5月に全島避難となる噴火が発生し,その後2017年,2018,2019年にも爆発的噴火が散発的に発生している.草津白根山・本白根山では2018年1月に小規模な水蒸気噴火が発生した.伊豆・小笠原諸島の西之島で2013年11月から始まった噴火は,周辺の浅海を溶岩で埋め立て新しい火山島を作り出し,約2年の活動後一旦終息したが,2017年4月と2018年8月に再度活発化し流出した溶岩により西方と南方への拡大が進んだ.2019年12月には再度溶岩流出が始まり,2020年8月まで続いた活動により旧島部分は完全に埋まり,面積が大幅に増加した.
本センターでは主たる観測対象火山以外についても様々な観測研究を行っており,新たな火山活動があれば国内外を問わず観測あるいは調査研究を実施している.また,実験・理論,シミュレーション,地質学,物質科学的手法等に基づく火山の基礎研究も実施している.以下,火山毎および研究手法毎に,最近の主たる研究成果を紹介する.
(1)地震・電磁気・測地観測網(海半球観測ネットワーク)の展開・維持
(1-1)海洋島地震観測網
ジャヤプラ(インドネシア),パラパト(インドネシア),デジャン(韓国),ポナペ(ミクロネシア),マジュロ(ミクロネシア),犬山(日本),石垣(日本),パラオ(パラオ),バギオ(フィリッピン),父島(日本),カメンスコエ(ロシア),サパ(ベトナム),ハイフォン(ベトナム),ビン(ベトナム)の9ヵ国14定常観測点における観測を, 海洋研究開発機構と共同で継続した.このうちマジュロ(ミクロネシア),父島(日本),カメンスコエ(ロシア)を除く11観測点からはリアルタイムで地震波形データを収集した.
(1-2)海洋島電磁気観測網
ポナペ(ミクロネシア連邦),アテーレ(トンガ王国),モンテンルパ(フィリピン),カンチャナブリ(タイ),ワンカイヨ(ペルー),南鳥島の各観測点における地磁気3成分と全磁力の観測を海洋研究開発機構と共同で継続した.絶対観測値を用いて2017年以降の地磁気三成分確定値の検討を開始した.また,2019 年までの観測値の公開準備を行った.通常はポナペおよびアテーレ観測点において年に1度の頻度で絶対観測と関連する観測を行うが,COVID-19 による各国の渡航制限により,この観測が実施できなかった.
(1-3)海底ケーブルネットワークによる電位差観測
フィリピン-グアム,二宮沖-グアム(TPC-1),グアム-沖縄(TPC-2),上海沖-苓北(上海ケーブル)の海底ケーブルについて電位差観測を継続し,これらの電位差に含まれる長期変動成分の解析を継続して行った.電位差成分の永年変動(時間1階微分)と,短期主磁場変動の地磁気ジャークや海流変動との関連の調査を継続した.また,電位差変動から地下電気伝導度構造の推定を目的として,海洋潮汐による電磁誘導数値モデリング手法の開発も継続して行なうとともに,海洋潮汐起源信号の強度と環境変化に伴う海水電気伝導度変化の関係に関する研究を開始した.
(2)海半球観測網を補完する長期アレイ観測
(2-1)海底地震観測
海底観測網直下の構造を浅部から深部まで決定する「広帯域海底地震探査」の手法開発を継続して行った.周期3–30秒においては地震波干渉法を,周期30–100秒においては遠地地震のアレイ解析手法をもちいることで,地震波異方性も含めた深さ10–150 kmの構造の定量的な議論が,浅部の構造を仮定せずに行うことが可能となった.既存の海底観測アレイにこの手法を適用する中で,使用するデータの選別の重要性が明らかになった.また,海洋底を伝播する Love波は複数のモードがほぼ同時に到着するためモード分離が困難であったため、位相速度測定が困難であった.観測波形を基本モードと1次高次モードの合成波として扱う新手法を開発し,Love波の基本モード位相速度の測定を可能とした.
本センターが実施した海底地震観測の記録は,「ふつうの海洋マントル計画」までの記録がOHPデータセンターより公開済みである
(2-2)海底電磁気観測
三陸沖日本海溝では,太平洋プレートの沈み込みに伴う変遷と地震発生との関連を電磁気学的手法と熱学的手法で解明することを目的とした研究を,2007年よりJAMSTECと共同で進めた.またこの海域での観測は,2009年度以降は,「地殻流体」計画の一環として継続している.2012年度までに海溝軸を横切る複数の測線上の合計31観測点でデータを取得し,2次元構造解析を進めている.なお,本研究で2010年に設置したOBEMは,2011年3月11日の東北沖地震に伴って生じた大津波によって誘導された磁場変動を記録しており,巨大振幅津波の波源域推定に貢献した(Ichihara et al., 2013, Earth Planet. Sci. Lett.).更に2013年4月から8月にかけて,新潟・秋田県沖日本海でも6台のOBEMを用いた観測を行った.同時に周辺の島で観測したデータ,過去に秋田県沖日本海で取得したデータを加えて3次元解析が進行中である.これらの観測データを統合的に解析し,最終的には日本海溝から日本海にかけての島弧断面の電気伝導度構造を明らかにすることを目指している.
(2-3)陸上電磁気観測
1998年以来,中国地震局地質研究所の協力を得て中国東北部吉林省中部および遼寧省西部・中部においてネットワークMT観測を行ってきた.そのデータの解析から,吉林省内の4地点においてマントル遷移層の深さで電気伝導度が他地域に比べて有意に高くなる傾向が認められた.ただ,その解析において,深部構造を決定する鍵となる数日以上の長周期データは,長春の磁場観測に基づく鉛直磁場-水平磁場変換関数のみを用いていたという問題点があった.このため,上記の変換関数の空間的な分布特性を調べるために,2007年より,中国全域にわたる既存磁場データのコンパイルと解析を始め,周期数日から100日程度の超長周期の変換関数推定を試み,誤差の小さな質の良い応答関数を推定した.その応答関数に基づき,1次元層構造を仮定した構造推定を試みた.その結果,中国東北部の広域にわたってマントル遷移層が高い電気伝導度をもつことが明らかとなった.しかし一方で,特に低磁気緯度地域で,下部マントルに至るまで異常に低電気伝導度となる結果が得られた.このため,昨年度に引き続き,3次元順計算によって,大規模な低比抵抗域に電流が集中することが磁場ー磁場変換関数にどのような影響を与えるかを見積もった(地震予知研究センターと共同).
(3)海半球ネットワークデータの編集・公開
Boulder Real Time Technologies社のAntelopeというソフトウェアを用い,オーストラリア地質調査所,台湾中央研究院地球化学研究所,及びIRISとリアルタイムデータ交換を継続した.
各種機動観測データの公開を継続した.定常観測点データに関しては,海洋研究開発機構と共同で,広帯域地震データ,GPSデータ,電磁気データの公開を継続した.
(1)次世代の海底地震・測地観測システムの開発
本所において共に海域地震観測を行う観測開発基盤センターと共同し,海底地震観測の高度化として複数次元での観測帯域拡大を進めている.現在,広帯域地震観測での機器の高機能化,機動的海底観測での測地学的帯域への拡大,および水深6000 mを越える超深海域での地震観測の実現,の3項目を具体的課題として機器開発を実行中である.
広帯域海底地震計(BBOBS)の平均的ノイズレベルを評価すると,長周期側での水平動のノイズレベルが陸上観測点での統計的上限に対して数倍以上高い.この対策として,低背なセンサー部をデータ記録部から独立させ海底面に突入させて自己埋設する構造の新型広帯域海底地震計(BBOBS-NX)を,ROV等の潜水艇による支援(設置・回収時)を要する運用方式で実用化した.2010年以降での複数の観測結果から,陸上観測点並のノイズレベルが確保できることを確認した.更に,このBBOBS-NXと同等の観測がROVを使用せず,自律動作により可能となる次世代機(NX-2G)の開発研究を科研費基盤研究(A)の補助を受け2015年から進めた.2016年10月にNX-2G試験機での実海域試験を実施,2017年4月に福島県沖日本海溝陸側斜面にて,既設置のBBOBS近傍にNX-2G試験機を設置,長期試験観測を開始し,2018年10月に回収,基本的な自律動作の機能が想定通りであることを検証した.更なる改良によるデータの質向上を確認するため、再試験を2023年度以降に実施する予定である.
また,BBOBS-NXを基に,機動的に広帯域地震・傾斜同時観測を行うBBOBST-NXの開発・実用化を進めると共に,海底での条件次第ではこれまでのBBOBSでも傾斜変動が計測可能であることも,複数地点での試験的観測データにより分かってきた.使用している広帯域地震センサーの長期間での安定性には問題は無さそうで,観測対象次第では有用と考えられる.2020年10月に,房総半島南東沖に2015年7月に設置したBBOBST-NXを5年ぶりに回収し[図3.7.4],2年間の地震・傾斜連続データを得ることに成功した.水温データも4年間分を取得した.これらのデータの解析処理を進めたところ,過去の傾斜観測時には不明瞭であったBBOBST-NXのセンサー部での温度依存性が明確になり,今後はセンサー内部で精密な温度記録を得て,より高精度な傾斜変動データの取得を狙う.なお,上記のNX-2Gでも傾斜観測は可能であり,機動的で高密度な海底地震・地殻変動観測アレイの実現性が見えてきた.
(2)最先端の海底電場観測装置(EFOS)の開発
電磁気探査の到達可能深度は,測定する電磁場変動の周期によって制御される(表皮効果).OBEM観測データのインバージョンによる最大探査深度は,周期1日以上で電場のS/Nが悪くなるために上部マントルの数百kmに限定される.新しい長基線電場観測装置(EFOS)は,長いケーブル(EFOS-6は6km,EFOS-2は2km)を海底に展張して良質な長周期電場データを取得する目的で開発された.上記「ふつうの海洋マントル計画」では,海域Aに合計3台のEFOS-2と1台のEFOS-6を設置し,2014年9月に3台のEFOS-2を,2015年9月に1台のEFOS-6を回収した.観測点NM16に設置したEFOS-2[図3.7.5]とOBEMの電場データのノイズスペクトルを比較すると,105秒よりも長い周期でEFOS-2のノイズが約1桁低いことが示された.このデータを用いて遷移層の電気伝導度を求め,地震波のレシーバ関数解析結果と統合して,遷移層に存在しうる水の量の上限を推定することができた.
今後進めるべき方向の一つは,EFOSによる観測を世界中の様々な海域で実施して,遷移層の水のグローバルな分布を明らかにすることである.しかし現状のEFOSは,設置および回収に無人探査機(ROV)を必要とし,このことがEFOS観測のグローバル展開を困難にする要因となっている.現在のEFOSは耐圧容器にガラス球を用いているために深海有人探査機での取り扱いができない.2017年度はこの点を改善して有人探査機でも扱えるよう,耐圧容器を金属製に変更した.2018年度および2019年度には,科研費基盤研究(B)により有人探査機による展張・回収システムを検討し,作成した.このシステムを有人潜水調査船「しんかい6500」により設置する観測航海がこれまでに2度採択された(2019年8月・小笠原海盆,2021年6 月・伊豆諸島青ヶ島東方沖)が,いずれも海況の不良により,機器設置を行うことができなかった.次回の設置機会に向けて,さらなる開発と観測準備を継続している.
深海でのEFOSの設置・回収作業が可能な有人/無人探査機は世界中を見ても,極めて数が限られる.一方,マニピュレータがないため複雑な作業はできないが,深海底でケーブルを展張する機能はある各種曳航体が使用可能な研究船は,多くの国で保有している.これらの曳航体を用いた設置・回収が可能になれば,EFOSによる観測の機会が格段に増えることが期待される.我々は,深海曳航体(ディープトウ)によって設置/回収できるよう,EFOSの全面的設計変更を行い,このシステムについても本格的な開発を進めている.
(a)ミュオグラフィ検出器 - 並列ミュオグラフィの強化
2006年に地震研究所が火山内部を世界に先駆けて描き出して以来,ミュオグラフィは急速に世界に広まりつつある.ミュオグラフィとは,宇宙線に含まれる高エネルギー素粒子・ミューオンの強い透過力を利用して,キロメートルを超えるサイズの巨大物体内部を透視し,その内部の密度構造を可視化する技術である.これまで第2世代システムのノイズ低減能力を強化することで2013年に薩摩硫黄島で発生した噴火において,マグマの昇降をとらえることに成功しているが,薩摩硫黄島は小規模火山として位置付けられるため,ミュオグラフィを桜島のような中規模火山に適用しようとすると,より厚い岩盤を通り抜けることができる極めて低強度のミューオンを一定時間内にできるだけ多く記録する必要がある.そのために2014年に設置された桜島ミュオグラフィ観測所(SMO)を観測装置の並列化により継続的に強化してきた.
2015年から2017年にかけて学術交流協定,知的財産協定など種々の協定を締結してきたハンガリー科学アカデミーウィグナー物理学研究センターとの協働により,2017年には軽量高解像度ミュオグラフィ観測システム(Multi-wire-proportional-chamber-based Muography Observation System; MMOS)を開発した.これは軽量でありながらも第2世代システム以上の高いノイズ低減能力と従来技術を一桁以上凌駕する解像力を実現した.ただ,有感面積が不十分であったため,2018~2019年にかけて口径を順次拡大し,現在では5.9m2となっている.2019年度はこれをさらに拡大し,2020年に入るまでに総有感面積は9m2に到達した.また,2019年度には並列化に起因する故障率を低減する目的で複数台の観測装置すべての通信系統を無線化することで通信故障率が軽減されたが,2020年度は電気系統においても,安定運用を妨げる要因があることが明らかとなり,その対策を講じている.
一方,並列化の段階で得られたデータについても解析・解釈が進んだ.2017年終わりから2018年初めにかけて桜島における噴火が昭和火口から南岳火口へと推移したが,それに合わせて観測された昭和火口底直下における直径200m程度の密度上昇現象について考察を行い,それがプラグ様の物体であることが分かった.2020年度も引き続き後継を拡大することで時間分解能を上げ,時系列画像を取得していった結果,南岳火口下にプラグ形成を示唆する高密度構造物の成長が見られた.このプラグは南岳火口の活発化に伴って形成されつつあるものであることが想定されるが,今後更に時間分解能を上げた解析によって,切迫性評価にどう活用できるか引き続き火山学の各分野の研究者とさらに連携して検討していく.
さらに2021年度には,活発化した南岳火口付近におけるミュオグラフィを活用したテフラ(火山性噴出降下物)の増減モニタを開発した.テフラが堆積することで火山表面における土壌の安定性が弱まり,重力に駆動される土石流などが誘引されやすくなる.桜島において図3.8.1 A,Bに示す5つの領域をモニターした結果が図3.8.1 Cである.南岳火口付近におけるテフラの増減が最も顕著なことがわかるが,それ以外の部分でも同様の変動が見られ,幅広い領域において南岳火口噴火に伴うテフラの影響があることが考えられる.この結果は,ミュオグラフィを用いることで比較的短時間でテフラの堆積・浸食プロセスをモニタリングできることを示しており,火山表面のモデリングや火山の警戒レベルのアセスメントへの貢献が期待される.
(b)ミュオグラフィの海への展開
東京大学国際ミュオグラフィ連携研究機構(主管部局地震研究所)は,同大学生産技術研究所,大学院新領域創成科学研究科,および九州大学,関西大学,シェフィールド大学,英国科学技術施設会議ボルビー地下実験施設,ウィグナー物理学研究センター,日本電気株式会社と共同で,世界初となる海底ミュオグラフィセンサーアレイの一部を東京湾アクアライン海底トンネル内部の100 mにわたって設置し,東京湾における天文潮位のリアルタイム測定に成功した.
東京湾の海水を貫通し,海底下の東京湾アクアライン海底トンネルにまで到達した素粒子ミュオンは,センサーモジュールにて検知され,TS-HKMSDDの中央に位置するデータ収集センターにて記録される(図3.8.2).この記録されるミュオン数の時間変化を測定することにより,TS-HKMSDD上部に位置する海水の動きや海底岩盤内部の変化をイメージングすることが可能となる.
(c)ボアホール設置型ラジオグラフィー
宇宙線ミューオンは上空からのみ飛来する.したがって,断層破砕帯や地滑り面等の地下構造を透視するためには,測定対象を見上げるように,ミューオン検出器を地下深く掘削坑(ボアホール)等に埋設することが必要となる.しかし,ボアホールのような狭隘な空間では,センサーの有効面積を大きくとることが困難であり,ミューオン・フラックスは限られた量しか得られないので,それを有効に活用する観測技術の開発が不可欠となる.
2014年度までに,跡津川断層(岐阜県飛騨市の山中)近傍に掘削された最大深度350mのボアホールを利用して,深度100mまでのミューオン・フラックスデータを取得した.2020年度は検出器の較正と解析にとりくんだ.これまでは検出器の感度分布を完全に再現することができなかったが,原因がチャンネル間のクロストークにあることを突き止め,クロストークの効果を含めたシミュレーションとデータ解析ツールを開発し,問題を解決した.解析の結果,シミュレーションは観測データを良く再現し,得られた断層の姿勢(走向・傾斜角・深さ)は,過去の地質調査結果と良く一致した(走向:北から時計回り72.5±0.4°・傾斜角:斜面北向き85.1±0.4°・断層とボアホールとの交点:-57.5±3.3 m).また,観測結果から予想される断層露頭の位置と,過去の地質調査で見つかっている露頭の位置も一致した.加えて,断層破砕帯とみられる低密度領域の幅は,140±40 mという結果が得られた(図3.8.3)一般の断層と比べると非常に幅が広いが,過去の電磁気探査の結果とは調和的な結果となった.
跡津川断層の観測の結果,我々が開発したボアホール型ミューオン検出器は断層の透視を行う上で十分な性能を持ち、断層の姿勢及び破砕帯の幅の100mスケールでのリモートセンシングという,これまでの手法では不可能であった観測が可能となることが分かった.今年度は解析に加え,より小型な検出器の開発にも取り組み,電子回路の製造,シンチレータ検出器の製造を行った.今後は新たに開発した検出器を用いて、観測事例を増やし,防災研究への応用を行っていく.
分散型音響センシング(Distributed Acoustic Sensing,以下DAS)技術を用いて,四国中央部において超稠密な地震観測を2021年11月中旬から開始した.既設の光ファイバの末端にDAS装置を接続し,コヒーレントな光パルスを光ファイバに連続的に伝送し後方散乱信号を測定することで,ファイバ軸方向の動的ひずみデータを取得した.観測には,徳島県三好市池田町を起点に,国道32号線と国道192号線に沿う2本の光ファイバを使用した.たとえば,観測期間中の2021年12月30日にインドネシアで発生したMw7.3の地震による良好な波形データが取得された.さらに,測線近傍で発生した微小地震による双曲線形状の波面が光ファイバに沿ってシームレス且つ明瞭に記録されている.
太陽系の岩石惑星の中でも、地球は、海と陸、活発な地震・火山活動、プレート運動と大陸移動、地球磁場を有し、生命を宿す「にぎやかな惑星」である。なぜ兄弟惑星である火星や金星と異なりこれほど活発で多様性に富むのか、その仕組み・鍵の一つは水にあると考えている。物質科学系研究部門・岩森研究室では、これらのユニークな地球の営み(=地球ダイナミクス)について、特に水と固体地球の相互作用に注目しながら、温泉や火山の調査、室内分析、データ統計解析、数値シミュレーションなど、さまざまなフィールドや手法を組み合わせて研究している。2020-2021年には、昨年度の内容(1.と2.)の発展に加え、新たに3.を実施した。
- 島弧火山岩および地下水組成解析(地球化学解析および統計解析)に基づき、沈み込んだプレートから物質が供給され,マグマや深部流体が地表に達するまでのプロセスを、特にカムチャッカ半島、箱根火山、阿蘇火山などでの岩石・流体試料とそれらの組成解析・統計解析に基づき、定量的にとらえた。
地球表層を覆うリソスフェアが,地球内部のアセノスフェアと熱的・物質的にどのように相互作用するかは,プレート運動や地球内部物質循環,および地球の熱的進化を規定する重要なプロセスである。アフリカ・カメルーンでは,大陸リソスフェアがプルームあるいはマントル対流の上昇流と相互作用し,火山(カメルーン火山列)を生み出していると考えられている。カメルーン火山列のマグマの組成およびその時空間から、リソスフェア-アセノスフェアの相互作用とその時間変化を明らかにした。
地殻-上部マントルでの液体(水溶液、超臨界流体、マグマ)の存在量や分布形態をとらえるため、観測される地震波速度と電気伝導度を同時解析する手法を開発した。ある与えられた、しかし多様な温度・圧力、岩質、液体の種類・組成と分布量・形態に対して、地震波速度と電気伝導度を予測するモデルを構築した。またこのモデルを用い、逆問題として、すなわち観測される地震波速度と電気伝導度に基づいて岩質や液相の量・分布形態が推定可能であることを示した。