日本地球惑星科学連合2022年大会

講演情報

[J] 口頭発表

セッション記号 S (固体地球科学) » S-SS 地震学

[S-SS06] 地震波伝播:理論と応用

2022年5月23日(月) 15:30 〜 17:00 301A (幕張メッセ国際会議場)

コンビーナ:澤崎 郁(防災科学技術研究所)、コンビーナ:西田 究(東京大学地震研究所)、新部 貴夫(石油資源開発株式会社)、コンビーナ:岡本 京祐(産業技術総合研究所)、座長:江本 賢太郎(東北大学大学院理学研究科)、西田 究(東京大学地震研究所)

15:30 〜 15:45

[SSS06-06] 火山性地震のエンベロープ幅を用いた逆問題解析(1):空間重み関数を用いた散乱・非弾性減衰の空間分布推定

*熊谷 博之1浜本 未希1 (1.名古屋大学大学院環境学研究科)

The envelope widths of natural and active seismic sources have been used to estimate one-dimensional (1D) layered seismic scattering and attenuation structures at volcanoes, which indicate that a surficial, highly heterogenous and attenuative layer up to around 1 km thickness commonly exists. The envelope width (p) is defined by the ratio of the cumulative to peak amplitudes for S waves in high-frequency (5−10 Hz) envelope seismograms. In this study, we developed an inversion method to estimate spatial distributions of the scattering mean free path (l0) and the quality factor for medium attenuation (Qi) using the residuals between observed p values and those calculated with 1D l0 and Qi models. One of the difficulties in performing such an inversion is the computational cost to calculate the partial derivatives of the envelope widths with respective to l0 and Qi perturbations in individual grids for all combinations of source-station pairs to construct the data kernel. To overcome this difficulty, we proposed an efficient approach to calculate the partial derivatives using the 3D space-weighting functions (SWFs) for the 1D models. In our inversion, we used the quality factor for scattering attenuation (Qs) instead of l0, and focused on the first surficial layer because the envelope width is insensitive to the second layer. We first calculated the energy distributions in the individual 3D grids for the 1D models by Monte Carlo envelope simulations for individual source locations. The resultant 3D energy distribution for each source location is regarded as the SWF. Then, we calculated the partial derivatives with respective to Qs-1 and Qi-1 perturbations in the entire first layer for all combinations of source-station pairs. These partial derivatives were weighted by the SWF to obtain the partial derivatives in the individual grids for each source-station pair, which were used to construct the data kernel. We solved the damped least-squares problem for the envelope width residuals with the hyper-parameters for Qs-1 and Qi-1 corrections in the individual grids. Our tests of this inversion method with synthetic envelope waveforms indicated that iterative inversions are required to suppress artifacts tended to appear beneath station locations. Our inversion results for the observed envelope width data at Taal (Philippines), Nevado del Ruiz, and Galeras (Colombia) indicate that there exist smaller l0 and larger Qi areas near the summit craters at all these volcanoes, which may correspond to heterogeneous, unconsolidated sediments with cracks filled with steam in volcanic conduits. Smaller Qi areas found beneath the eastern flank at Taal matched with an S-wave attenuation region, where the existence of actively degassing magma was inferred by a previous study. These results demonstrate that our inversion method using the envelope widths is useful to characterize spatial distributions of scattering and attenuation structures at volcanoes.