A simulation study of electromagnetic wave reflection and refraction to monitor subsurface electrical resistivity time variation

• Project/Area Number: 17K05624
• Research Category: Grant-in-Aid for Scientific Research (C)
• Allocation Type: Multi-year Fund
• Section: 一般
• Research Field: Solid earth and planetary physics
• Research Institution: Tohoku University
• Principal Investigator: Ichiki Masahiro 東北大学, 理学研究科, 助教 (80359182)
• Project Period (FY): 2017-04-01 – 2020-03-31
• Project Status: Completed (Fiscal Year 2019)
• Budget Amount: ¥3,510,000 (Direct Cost: ¥2,700,000、Indirect Cost: ¥810,000); Fiscal Year 2019: ¥2,340,000 (Direct Cost: ¥1,800,000、Indirect Cost: ¥540,000); Fiscal Year 2018: ¥650,000 (Direct Cost: ¥500,000、Indirect Cost: ¥150,000); Fiscal Year 2017: ¥520,000 (Direct Cost: ¥400,000、Indirect Cost: ¥120,000)
• Keywords: 地震発生予測・火山噴火予測 / 観測手法開発 / 火山噴火予測 / 地震発生予測 / 火山噴火予測・地震発生予測

Outline of Final Research Achievements

Electrical resistivity time variation of the earth's crust reflects transport of crustal fluids, which more or less involves phreatic eruption or seismicity. To develop a technique to monitor the resistivity variation, this project aims to study loop-loop transient electromagnetic method. The objective is (1) to develop a finite element time domain simulation code to calculate horizontal transient electromotive force in the 3D resistivity distribution, and (2) to calculate expected transient electromotive force in the 3D resistivity distribution model of Zao volcano, NE Japan by use of the code.
We took a finite element time domain approach to develop the code, and found still insufficient numerical precision of loop source formulation to be composed line dipole segments. In uniform half-space resistivity, we can expect 0.987 mV at 1 ms, 0.987 us at 10 ms on 100×100 m^2 in 100 Ohm-m, and 98.7 mV at 1 ms, 98.7 uV at 10 ms in 10 Ohm-m, induced by 100×100 m^2 loop of 10 A.

Academic Significance and Societal Importance of the Research Achievements

2014年の御嶽山の噴火や2016年の阿蘇山の噴火によって、地震計や傾斜計、GNSSによる力学的変動の観測では水蒸気噴火の短期的噴火時間予測が難しいケースが再認識された。水蒸気噴火の予測には火山体下の熱水流体の動きを捕捉すると考えられる電気比抵抗の変化の監視で精度を向上させることが期待されており、本研究では時間領域電磁法と呼ばれるうちLoop-Loop法と呼ばれる手法を向上させることで、地下の電気比抵抗が変化した場所を特定できるだけでなく、時間、空間的に分解能を改善させ、火山噴火予測、さらには流体駆動型の内陸地震などの予測を多項目で実装できることが期待できる。

Research Products

• [Presentation] A 3-D crust and uppermost mantle electrical conductivity model of subduction zone beneath NE Japan (2019)
  • Author(s): Ichiki, M., T. Kaida, Y. Ogawa
  • Organizer: 日本地球惑星科学連合
• [Presentation] A 3-D crust and uppermost mantle electrical conductivity model of subduction zone beneath NE Japan (2019)
  • Author(s): Ichiki, M., T. Kaida, Y. Ogawa
  • Organizer: The Second International Symposium Crustal Dynamics
  • Int'l Joint Research
• [Presentation] 吾妻山・安達太良山周辺の地殻電気伝導度構造 (2018)
  • Author(s): 市來雅啓, 海田俊輝, 中山貴史, 山本希, 福島洋, 三浦哲
  • Organizer: 日本地球惑星科学連合2018年大会
• [Presentation] EM ACROSS experiment at Kusatsu-Shirane volcano (2018)
  • Author(s): 曾國軒・小川康雄・國友孝洋・深井雅斗・Ted Bertrand・木下雄介・南拓人
  • Organizer: 平成29年度 Conductivity Anomaly 研究会
• [Presentation] Audio-frequency magnetotelluric, and total magnetic intensity observations in 2014-2016, at Zao volcano, NE Japan (2017)
  • Author(s): Ichiki, M., T. Kaida, T. Demachi, S. Hirahara, T. Nakayama, M. Yamamoto, S. Miura, T. Moriyama, M. Akutagawa, T. Kobayashi, S. Matsuura, S. Omori, K. Ono, S. Seki, W. Kanda, K. Seki, Y. Ogawa, M. Ushioda, M. Uyeshima
  • Organizer: AGU 2017 Fall Meeting
  • Int'l Joint Research