SSS27-01

会場:103

時間:5 月 24 日 10:00-10:15

地震発生の予測可能性があるかないか？ 〜前震活動の統計的性質に着目した検討
Statistical investigation of spatio-temporal densities of foreshocks to understand earthquake predictability
鴨川 仁 1∗ ; 田中 利佳 1 ; 織原 義明 1 ; 橋本 哲 1
KAMOGAWA, Masashi1∗ ; TANAKA, Rika1 ; ORIHARA, Yoshiaki1 ; HASHIMOTO, Satoshi111

東京学芸大学教育学部物理学科
Dpt. of Phys., Tokyo Gakugei Univ.

The relation between the size of earthquake preparation zone and the magnitude of forthcoming earthquake is different between
nucleation and domino-like cascade models. The former model indicates that the magnitude is predictable before the mainshockof the earthquake, because the preparation zone is proportional to the rupture area. One the other hand, the latter indicates that themagnitude is unpredictable, because the rupture consisting of sequence of tiny earthquakes is unknown to terminate. Since thisissue is still controversial, we would like to verify the two models using the methodology proposed by Lippiello et al. (Scientificreports, 2012). In the analysis, spatial occurrence rates of the foreshock and the aftershock are statistically compared. The results
show that both the rates are similar and the distribution of the rates versus the epicentral distance depends on the magnitude of
the mainshock. From the interpretation of these results, the nucleation model seems reliable.
キーワード: 地震, 前震, 余震
Keywords: Earthquake, Foreshock, Mainshock

1/1

SSS27-02

会場:103

時間:5 月 24 日 10:15-10:30

稀な現象の確率予測にも使える適切な採点法
Proper scoring systems available for probability forecasts targeting rare phenomena
林 豊 1∗HAYASHI, Yutaka1∗
1
1

気象研究所
Meteorological Research Institute

新たな手法を導入して予報を確実に改善できるためには、予報精度を適正に評価する採点法が存在し、新手法を導入した場合に出すであろう予報の採点結果（スコア）が従来の方法による予報のスコアを上回ることが必要条件である。
「適正な」採点ルールは、予報者が自分の判断通りに予報を出すことを奨励する必要がある（Murphy and Epstein, 1967）。降
水確率予報などの評価に用いられているブライアスコア（Brier, 1950）はこの数学的な条件を満たすが、カルバック・ラ
イブラー情報量（Kullback and Leibler, 1951）などはこの条件を満たさない。
「不適正な」採点ルールの下では、予報を細工した方が高いスコアを得られることがあるので、スコアを高める努力は、一般にはよりよい予報手法を目指す正直な行為とはいえない。このため、特に地震予知の分野で広く適用されている情報量等を尺度とした「不適正な」予報精度の測定や比較の意義の解釈には、慎重になるべきである。
ところで、ブライアスコアは予報値と実況値の平均自乗誤差の期待値で、現象の有と無に等しい重みを置いている。一方、稀な現象が対象の予報利用者にとっては、予報値が高い確率の場合や現象が実際に発生した場合の予報が重要な情報である。以下に、地震予測のように稀な現象を対象とした確率予報を適切に評価できる採点法を導出する。
予想値（予報者の真の信念に基づく確率）を p、予報値を f、実況値 i（現象が発生すれば 1、発生しなければ 0）、基準予報値（例えば統計値に基づく単純な予報値）を c、スコアを Si (f,c) と定義する。適正な採点ルールでは、期待値
Es (x,c)=pS1 (x,c)+(1-p)S0 (x,c) (式 1) について、Es (p,c) >Es (x,c) for all x ≠ p (式 2) が必要である。また、稀な現象をより重視する場合でも、基準値から考えられる予報の難易度について公正な条件 S0 (c,f) ≡ S1 (1-c,1-f) (式 3) と、難易度が高い予報ほど高いスコアが得られる条件∂ S1 /∂ f|c=const >=0, ∂ S1 /∂ c|c=const <=0 (式 4) の成立が必要である。
式 1,3 を満たす解は式 (5) を満たす。ここで A は任意の関数である。
S1 (f,c)=-(i-f)B’(f,c)+(i-x)B’(x,c)-B(f,c)+B(x,c), A=d2 B/df2 (式 5)
適当な境界値と A の関数形で式 5 を解き、Si が不等式 2 と 4 を満たせば、適正な採点ルールを満たし、かつ、都合のよい性質を持つ評価法を得られる。
例えば、基準値予報のスコア Si (c,c) ≡ 0、完全な予報のスコア Si (i,c) ≡ 1 の境界条件で解くと、A=-2 で Si (f)=1-(i-f)2(式 6) を得る。ここで、1-Si (f) はブライアスコアに一致する。
もう一例として、基準となる予報値でのスコア Si (c,c) ≡ 0、完全な予報の期待スコア Es,p ≡ 1 を境界条件として A が
f の 0 次式の解を求めると、A=-2 / c(1-c) から Si (f,c)={(i-c)2 -(i-f)2 } / c(1-c) (式 7) を得られる。これも適正な採点ルールを満たす。c=1/2 の時、(1-Si (f))/4 はブライアスコアであるから、式 7 は基準値予報の選択の任意性に関してブライアス
コアを拡張した評価式と解釈してよい。
以上のように、基準予報値を参照した採点式にすることで、地震予測のように稀な現象を対象とした確率予報の適正な採点法の一般解と、特に、稀な現象により重みを置いて適正に評価できるいわゆる拡張ブライアスコアと呼ぶべき解を導いた。
大会では、式の導出の詳細と、これら採点法の具体的な適用に向けての課題も議論したい。なお、津波警報のように現象の有無を対象とする二値予報の一種を適切に評価する採点法は、効用理論から導出し議論済みである（林,2014, JpGU）。
キーワード: 基準値予報, 地震予測, 拡張ブライアスコア, 情報利得, 確率予報, 適正な採点ルール
Keywords: benchmark forecast, earthquake prediction, extended Brier’s score, information gain, probability forecst, properscoring rule

1/1

SSS27-03

会場:103

時間:5 月 24 日 10:30-10:45

大地震前の地殻活動に伴う異常現象の物理的解釈と検知
Physical interpretation and detection of anomalies associated with crustal processes leading to large earthquakes
末広 潔 1∗ ; サックス セルウィン 2 ; 高波 鐵夫 2 ; ライデレク ポール 2 ; スミス デボラ 2 ; 波多野 恭弘 3
SUYEHIRO, Kiyoshi1∗ ; SACKS, Selwyn2 ; TAKANAMI, Tetsuo2 ; RYDELEK, Paul2 ; SMITH, Deborah2 ;HATANO, Takahiro3
1

海洋研究開発機構, 2 米国ワシントン市カーネギー研究所, 3 東京大学地震研究所
Japan Agency for Marine-Earth Science and Technology, 2 Carnegie Institution of Washington, 3 Earthquake Research Institute,
University of Tokyo1

自然地震の経験則、観測を説明するセルラーオートマタモデル（サックスとライデレク、１９９５）の考えに基づき、地震予知（短期予知）への観測的アプローチを提案する。このアプローチはプレサイスミックスリップの検出とは異なってダイラタンシー現象の捕捉を考える。これまでに報告されているマグニチュードに依存する地震活動静穏化（たとえば１９８２年浦河沖地震、１９９４年ノースリッジ地震、２００８年四川地震など）は、ダイラタンシー硬化過程の進行を、このモデルに組み込むと再現できる。しかし、この情報は直前の変化ではなく年オーダーの変化であるので、注意の喚起までである。この変化は、地震域を囲む大きな領域に現れる傾向にあることが特徴である。ダイラタンシーの考えを進めると、地震直前の観測異常はダイラタンシーによる微小な割れ目群に流体が移動する過程に加速的な作用が加わって現れると解釈できるが、物理モデルに照らして量的に検討したうえで、観測の方法として垂直歪成分観測を新たに注目すべき観測パラメタとなる可能性を示す。
キーワード: 地震活動静穏化, ダイラタンシー, セルラーオートマタ地震モデル
Keywords: seismic quiescence, dilatancy, cellular automata earthquake model

1/1

SSS27-04

会場:103

時間:5 月 24 日 11:00-11:15

体積歪みによる群発地震の予測
Predicting changing rates of swarm activity by volumetric strain
熊澤 貴雄 1∗ ; 尾形 良彦 1 ; 木村 一洋 2 ; 前田 憲二 2 ; 小林 昭夫 2
KUMAZAWA, Takao1∗ ; OGATA, Yosihiko1 ; KIMURA, Kazuhiro2 ; MAEDA, Kenji2 ; KOBAYASHI, Akio2
1
1

統計数理研究所, 2 気象研究所
The Institute of Statistical Mathematicsa, 2 Meteorological Research Institute

Near the eastern coast of Izu peninsula is an active submarine volcanic region in Japan, where magma intrusions have been
observed many times. The forecast of earthquake swarm activities and eruptions are serious concern particularly in nearby hot
spring resort areas. It is well known that temporal durations of the swarm activities have been correlated with early volumetricstrain changes at a certain observation station of about 20 km distance apart. Therefore the Earthquake Research Committee(2010) investigated some empirical statistical relations to predict sizes of the swarm activity. Here we looked at the backgroundseismicity rate changes during these swarm periods using the non-stationary ETAS model (Kumazawa and Ogata, 2013, 2014),
and have found the followings. The modified volumetric strain data, by removing the effect of earth tides, precipitation andcoseismic jumps, have significantly higher cross-correlations to the estimated background rates of the ETAS model than to theswarm rate-changes. Specifically, the background seismicity rate synchronizes clearer to the strain change by the lags arounda half day. These relations suggest an enhanced prediction of earthquakes in this region using volumetric strain measurements.
Hence we propose an extended ETAS model where the background rate is modulated by the volumetric strain data. Here wehave also found that the response function to the strain data can be exponential functions with the same decay rate, but that theirintersects are inversely proportional to distances between the volumetric strain-meter and the onset location of the swarm. Our
numerical results by the same proposed model show consistent outcomes for the various major swarms in this region.
キーワード: 非定常 ETAS モデル, 常時地震活動, 群発地震, 体積歪み, 予測
Keywords: nonstationary ETAS model, background seismicity, swarm, volumetric strain, prediction

1/1

SSS27-05

会場:103

時間:5 月 24 日 11:15-11:30

Seismo-ionospheric precursor monitoring system based on near-real-time spaceborne andground GPS observation
Seismo-ionospheric precursor monitoring system based on near-real-time spaceborne andground GPS observation
TSAI, Ho-fang1∗ ; LIU, Jann-yenq2 ; CHEN, Chia-hung1 ; SUN, Yang-yi2 ; HO, Tiffany3
TSAI, Ho-fang1∗ ; LIU, Jann-yenq2 ; CHEN, Chia-hung1 ; SUN, Yang-yi2 ; HO, Tiffany3
1

Department of Earth Sciences, National Cheng Kung University, Tainan, Taiwan, 2 Institute of Space Science, National Central
University, Jhongli City, Taiwan, 3 Taiwan Analysis Center for COSMIC, Central Weather Bureau, Taipei, Taiwan
1
Department of Earth Sciences, National Cheng Kung University, Tainan, Taiwan, 2 Institute of Space Science, National Central
University, Jhongli City, Taiwan, 3 Taiwan Analysis Center for COSMIC, Central Weather Bureau, Taipei, Taiwan
Recently, the global ionosphere map (GIM) has been used to study the seismo-ionospheric precursors (SIPs) intensively. In
order to shorten the data latency, an SIP monitoring system is built up based on the near-real-time GIM processing. The GIM
data is derived from the combination of the ground-based and spaceborne total electron content (TEC) observation by meansof the spherical harmonic function, where the data is retrieved from a global GPS observational network and the FORMOSAT3/COSMIC radio-occultation (RO) experiments. The temporal statistical analysis is developed to detect the SIPs at several
important metropolitans such as Tokyo. Some new finding and results are to be further discussed. The spatial analysis will be
introduced to finding the repeat, duration and distribution of worldwide SIPs to estimate the possibility of forthcoming large
earthquakes in the future.
キーワード: global ionosphere map, seismo-ionospheric precursor, total electron content
Keywords: global ionosphere map, seismo-ionospheric precursor, total electron content

1/1

SSS27-06

会場:103

時間:5 月 24 日 11:30-11:45

2014 年長野県北部地震に先行した長期的温泉水温変化
Long-term groundwater temperature change at a hot spring preceding the 2014 Naganoken Hokubu earthquake of M6.7
佃 為成 1∗TSUKUDA, Tameshige1∗
1
1

なし
None

2014 年 11 月 22 日の長野県北部地震 (長野県神城断層地震) は，1986-90 ごろからその地震像が描き出されていた（地震予知総合研究振興会，1990，地震テクトニクスに関する総合研究 第３章 北部フォッサマグナ周辺に発生可能な地
震像，83-142.）．その主な根拠は，
1) 歴史上の地震（1714 年 M6.3）の存在．
2) 活動度 A の活断層（神城断層）の存在．
3) 1714 年の地震以後，白馬付近では目立った地震活動はない．長期の空白域である．
4) この糸魚川・静岡構造線に沿った地帯は，連鎖的に地震が発生する（1714, 1858, 1890, 1918 年の各地震）．
5) 1986 年の地震（M5.9）が前駆的な地震活動である可能性がある．
6) 90 年間の三角測量（三辺測量）によると地殻歪が集中的に蓄積している．最大せん断歪は 4〜6x10−５ (5〜7x10−７ /yr)．
先の地震（1714 年）から蓄積したとすると，緩和された分を考慮しても，地震がもうすでに発生してもよいという確率は高い．東西方向の縮みが確認されている．歪蓄積の範囲から規模を推定すると M6〜7．
白馬地域を主なターゲットとした 1995-2008 年の東京大学地震研究所の共同研究プロジェクト「内陸直下地震の予知 - 地震研究所特定共同研究 A」では種々の観測がおこなわれた．その中で，白馬倉下の湯温泉の源泉（白馬鉱山 HR-1
号泉）に於いては 1998 年 10 月から現在まで水温の連続観測がおこなわれている．泉源からガス圧で自噴した温泉水は、
まず 4t タンクに貯蔵されてガス抜きされる。このタンクから 90t タンクへ向けてポンプアップが行なわれる．4t タンクからのオーバーフローの水を流す排水路のマンホール (湯溜) に白金抵抗体水温計を設置している．
水温グラフのだいたいの傾向に注目すると，地震の 5 年ほど前から，それまではほぼ一定 (若干下降 -0.17 ℃/year)
だった水温に明瞭な下降傾向が見える (-1.5 ℃/year)．これは，震源域の神城断層付近では地震の前からねじれ歪が生じ，断層の西側，下盤側に位置する倉下の湯直下では引っ張りが働き，岩盤膨張，圧力低下，上昇流体の量減少，そして水温低下となったと考えると理解できる．温泉は地震後自噴が停止した (12 月 18 日には回復)．これも地震による，より大きな岩盤膨張のためと理解できる．
キーワード: 2014 年長野県北部地震, 地震予知, 前兆, 水温, 岩盤膨張, 岩盤収縮
Keywords: the 2014 Nagano-ken Hokubu earthquake of M6.7, earthquake prediction, precursor, water temperature, dilatation,
contraction

1/1

SSS27-07

会場:103

時間:5 月 24 日 11:45-12:00

前震活動に基づく地震発生の経験的予測−長野県北中部地域への適用−
Empirical forecast of mainshocks based on foreshock activities - Application to the northcentral Nagano prefecture 前田 憲二 1∗ ; 弘瀬 冬樹 1
MAEDA, Kenji1∗ ; HIROSE, Fuyuki1
1
1

気象研究所
Meteorological Research Institute【はじめに】 2014 年 11 月 22 日 22 時 08 分に長野県北部の神城断層付近で M6.7 の地震が発生し（最大震度 6 弱）、
負傷者や住家全半壊などの被害を生じた。気象庁の一元化カタログによると、この地震の約 4 日前の 11 月 18 日から 19日にかけて、本震のごく近傍で M3 クラス以下の地震が 40 個程度発生しており、小規模ながら活発な前震活動があった
ことを示している。しかし、こういった前震活動と本震との因果関係は物理的には解明されておらず、本震発生前に前震を識別し、本震発生を高い精度で予測することは一般に困難である。一方で、続発的に地震が発生しやすい特定の地域
では、本震前に活発な地震活動がみられる場合があり、その統計的性質から本震発生を経験的に比較的効率よく予測できる場合がある。そのような地域として、Maeda（1996, BSSA）、前田・弘瀬（2012 年秋、地震学会）や Maeda & Hirose（2014 年春、合同大会）は日本海溝沿いの３領域や伊豆地域を指摘し、これらの領域における前震を基にした予測効率を明らかにしてきた。本発表では、これまで提案してきた手法を上記の長野県北部の地震を含む長野県北中部の領域に適用し、この地域ではどのようなまとまった地震活動を前震活動とみなせば、過去の事例から経験的に本震（ここでは M5以上とした）の予測に効率的であるかについて調査したので報告する。
【予測手法】 前震識別の手順は、１）震源カタログから本震とのマグニチュードの差が 1 以上の余震活動を除去し、
２）特定の大きさ（緯度 D°×経度 D°）のセグメントの中で、特定の規模（Mf）以上の地震が特定の期間（Tf 日）の間に特定の数（Nf）だけ発生した時、前震（群）の候補とみなし、３）その後特定の期間（Ta 日）内に本震が発生した場合に真の前震（群）であったと判定する。４）この前震（群）の定義において、D, Mf, Tf, Nf, Ta をパラメータとして、本震を予測するために効率のよい前震（群）のパラメータをグリッドサーチにより求める。予測効率の指標としては、本震は時空間的に一定の発生率を持ってランダムに発生すると仮定したモデルを基準とし、前震候補による予測時空間における本震発生率が他の時空間より高いとする予測モデルとの AIC の差（dAIC）を主に用い、その他、予知率（AR：予
測された本震の割合）や適中率（TR：真の前震の割合）、確率利得（PG：背景確率に対する予測時空間の発生確率の比）も補助的に用いた。
【データおよび解析結果】 1998 年から 2014 年までの気象庁の震源カタログから M1 以上、深さ 30km 以浅の地震に対し、D（0.1, 0.2, 0.3）, Mf(1, 1.5, 2, 2.5, 3), Tf(1, 2, 3, 5, 10), Nf(1, 2,...,20), Ta(1, 2,...,30) の 45000 ケースについてグリッドサーチを行った結果、長野県北中部地域（35.6° N − 37.1° N, 137.2° N − 139.0° E）の M5 以上の本震を予測するには、一辺 0.1°の大きさの矩形領域内に 1 日間に M2 以上の地震が 5 個発生した場合を前震候補とし、その後 5 日間にM5 以上の地震（余震は除く）が発生するという予測が予測効率として最も良いことが分かった。このパラメータ値を選択した場合、予知率は約 45 ％ (=5/11)、適中率は約 12% (=8/69)、PG=333、dAIC=66 であり、冒頭の長野県北部の地震も予測された地震に含まれる。言い換えると、長野県北中部においては、2014 年 11 月 22 日の本震の約 4 日前からみられたような（上記のパラメータで定義される）地震活動が M5.0 以上の本震に至る割合は約 12%であり、M5.0 以上の地
震の約 45%はそのような（上記のパラメータで定義される）前震活動を伴う、といえる。これまでの調査で予測効率が比較的高いことが分かっている日本海溝沿いの 3 領域〔予知率=約 38 ％ (=11/29)、 適中率=約 30 ％ (=13/44)〕や伊豆地域〔予知率=約 68 ％ (=44/65) 、適中率=約 23% (=46/196)〕に比べると、予測効率は、特に適中率の点であまり高くない
ことがわかる。このように、前震活動に基づく本震発生の予測効率は地域によって差がみられる。

キーワード: 地震予測, 予測効率, 前震, 長野県, 地震統計, 経験則
Keywords: earthquake prediction, performance, foreshocks, Nagano prefecture, statistics, empirical relation

1/1

SSS27-08

会場:103

時間:5 月 24 日 12:00-12:15

地震発生の１サイクルの間に生じる断層の非弾性変位
Estimation of inelastic displacement of a fault zone during an earthquake cycle
山本 清彦 1∗YAMAMOTO, Kiyohiko1∗
1
1

なし
none

１．はじめに：破砕帯／アスペリティーモデルは震源パラメターを断層帯の物理的性質として解釈することを目的に提案された (Yamamoto & Yabe, 2003, 2006 地震学会)．このモデルでは，アスペリティーを完全弾性体と見なして変形と破壊を扱うことによって，見かけの破壊エネルギーを見積もっている．さらに，破砕帯の回転によって生じる断層面に垂直な変位を一次近似で計算しているが，この近似は破壊エネルギーを過大評価することがわかった．ここでは，この垂直変位をアスペリティーの非弾性変形を考慮してより正確に再計算した．この結果に基づき，プレート境界の地震について，地震に伴う変位および再来時間と地震の規模との関係を検討した．
２．破砕帯／アスペリティーモデル：このモデルでは，アスペリティーを含む破砕帯を断層帯，断層とそれを囲む母岩との境界面を断層面とよぶ．アスペリティーは断層面間の変位が u∗c , u∗c =uc +uf c =t(ec +ef c ) になった時に破壊する．こ
こで uc と uf c は u∗c の弾性および非弾性成分であり，t は断層帯の厚さである．また，これらの変位はそれぞれ歪み ecと ef c に対応している．uc は臨界変位であり，アスペリティーの変位の弾性成分である．アスペリティーが破壊すると断層面には新たに ub が生じ地震を起こす．また，この時，断層帯の中にすべり面ができる．スリップ面での相対変位は
uc +ub になり，地震に伴う変位よりも大きい．
すべり面は断層帯中に伝播して拡大してゆくが，これに伴って，破砕帯中に回転が生じ断層面に垂直な変位 v を生じる．この変位は断層面に働く法線応力 sn に対して仕事をする．断層面の単位面積当たりの仕事を仕事密度 w とよぶと，仕事密度 w は， w=(sn v)}≈sn (1+2a)ec uc /2 と書ける．ここで a ＝ ef c /ec である．見かけの破壊エネルギーはこの仕事にほぼ等しい．震源の深さを約 10 から 20km とし，sn をその深さの静岩圧，岩石の剛性率を 30GPa とすると，破壊時の弾性歪みは約 10−2 である．これを使って得られる臨界変位と見掛けの破壊エネルギーの関係を大中と松浦（2002）によって纏められた実験結果と比較すると (1+2a)/2≈1 が得られる．これから，ecf ≈ec /2，すなわち，アスペリティーが破壊するときには弾性歪みの 1/2 の大きさの非弾性歪があることを示している．
３．プレート境界に発生する大きな地震の再来時間：1923 年関東地震と 1703 年元禄地震は相模トラフ沿いに起こった大地震として知られている．この間は 220 年である．前者の規模は M8.1 ± 0.2 とされるが，ここでは，M8.0 を採用
する．1703 年元禄地震の規模は 7.9 から 8.5 で，研究者によって異なる．M8.0 の地震を伴う平均断層変位は約 5.1m，臨界変位は約 2.56m である．元禄地震を基準にすると，断層帯の変位増加率として 11.6mm/year が求まる．非弾性成分が，常時，弾性成分の 1.5 倍あるとすると，全体の変位増加率は 17.4mm/year で Loveless, J.P. and Meade, B.J. (2010) のものにほぼ一致する．
プレート境界での弾性的な変位増加率を 11.6mm/year とすると，再来時間 400 年の地震の臨界変位は約 4.6m，断層変位の平均は約 9.2m，滑り面上の平均変位は約 16.1m になる．また，断層長は約 288km になり，相模トラフの全長にほぼ一致する．この地震の規模は約 M8.5 になる．

キーワード: 破砕帯／アスペリティーモデル, 非弾性変位, 破壊エネルギー, 臨界変位, 1923 年関東地震, プレート境界
Keywords: damage zone/asperity model, inelastic displacement, fracture energy, critical displacement, 1923 kanto earthquake,plate boundary

1/1

SSS27-09

会場:103

時間:5 月 24 日 12:15-12:30

発震機構オイラー回転による東日本大震災の本震判定
Euler Rotation of Focal Mechanism to determine the main shock of the East Japan Super
Earthquake, March 2011
新妻 信明 1∗NIITSUMA, Nobuaki1∗
1
1

静岡大学地球科学教室、仙台
Institute of Geosciences, Shizuoka University, Sendai

震源域の応力場変動を解明するために，観測された地震の発震機構の変化を定量的に検討する方法として，発震機構
回転角の算出法を考案し（http://www.niitsuma-geolab.net/ 特報 4），2011 年 3 月東日本大震災の本震域で起こった発震機
構変動から，応力場の変遷を解析した．
解析に使用したのは，気象庁がホームページで公開している，東日本大震災本震域で 2011 年に起こった地震の CMT
発震機構解である．
発震機構解として主応力軸（圧縮主応力 P 軸・引張主応力 T 軸・中間主応力 N 軸）の下半球方位が公開されている．
これらの主応力軸は互いに直交し，主応力軸間の角距離が 90°と一定しているので，主力軸方位の手系が一致していれば，二つの地震について 1 組の主応力軸方位を用いて，主応力軸方位の変化を球面上の 1 点であるオイラー極の回りの回転として表現可能である．二つの地震について主応力軸 PT・TN・NP の 3 組を用いて 3 組のオイラー極と回転角が算出できる．公開されている発震機構の主応力軸方位は度単位で公開されているので，算出されたオイラー極と主応力軸組について算出される回転角は必ずしも一致しない．そこで，主応力軸組の回転角差が最小となるオイラー極と回転角を採用した．最小回転角差は 0.1°程度である．
算出されたオイラー回転極が，海溝軸方位に直交する線よりも南側に位置すれば，発震機構座標頂が島弧側へ回転し，北側に位置すれば海溝側へ回転したことになる．島弧側への回転を「正」，海溝側への回転を「負」と，符号を付けて発震機構回転角に回転極性も表現できるように定義する．
東日本大震災の本震の発震機構を基準に，本震域で 2011 年に起こった地震の発震機構回転角を算出した．本震が基準なので本震の発震機構回転角は 0 となるが，2 月 16 日から開始して本震直前まで続いた 18 個の前震は全て正回転角であるのに対し，本震後の 35 個の余震は負回転角であることが判明した．発震機構回転角符号が本震を境に前震と余震で急変することは，それまで蓄積されてきた歪が本震によって解放されたことに対応する．日本海溝に沿う太平洋スラブ上面では，スラブ上面に沿う西傾斜の剪断応力とスラブ上面に直交する法線応力の合力が圧縮主応力 P 軸方位となる
ので海溝側に傾斜する．前震の発震機構が島弧側への正回転であることは，前震の圧縮主応力軸圧縮軸傾斜が本震の傾斜（36°）よりも小さいことを意味しており，最も破壊し難く法線応力が最大の破壊が本震であり，破壊し易く法線応力の小さい破壊が前震であることを示している．余震の発震機構が海溝側への負回転であるのは，本震によって剪断応
力は解放されたが，静岩圧はそのまま保持され応力場の主体となるため，圧縮主応力軸が垂直方向に向いて島弧方向に回転して正断層型になるためである．このように東日本大震災の発震機構回転方向の変化は，破壊力学的に説明できる．
3 月 9 日の最大前震 M7.3 の本震誤認が防災対策に深刻な影響を与えたが，発震機構回転角の符号に基づけば，前震と判定でき，本震が近付いていると予報可能であった．最大前震が前震であることは，翌 3 月 10 日の M5.2-6.8 地震 7 個の発震機構回転角絶対値がそれまでの前震と同様に 25°以内に収まっていることからも前震と判定できるが，符号付き
発震機構回転角を用いれば明確に判別できる．ちなみに，最大前震の発震機構回転角は+2.0°と本震の発震機構と殆ど同じであり，本震後の余震の発震機構回転は-25°よりもずっと大きな負回転である．
キーワード: オイラー極, 回転角, 発震機構, 手系, 本震判定, 東日本大震災
Keywords: Euler pole, rotation angle, focal mechanism, hand system, determination of main shock, East Japan Super Earthquake

1/1

SSS27-P01

会場:コンベンションホール

時間:5 月 24 日 18:15-19:30

スラブ全体の応力分布が 2004 年紀伊半島南東沖の地震後に一変した〜折り返し点通過、残り 50 年か〜
The stress distribution in the whole area of the slab changed completely after the 2004 off
Kii Peninsula earthquake
間瀬 博文 1∗
MASE, Hirofumi1∗
1
1

所属なし
none

(図面を参照のこと スラブ名・海底地形名等は便宜上命名した)(1) では 2004 年紀伊半島南東沖の地震と 1944 年東南海地震の関係を研究し以下の結果を得た。
東進する中国地方のマントルが南海スラブの縁を押す。スラブは右回転力とトラフ南側からの反作用の間に挟まれて圧縮されている。しかしトラフ南側は一様でなく杭のように抵抗する海山が二つあり反作用応力が集中している。2004
年の地震で「杭」の北側が破壊した後、「杭」へ力が伝わるルートが「北から」が「北西から」に転換した。
(2) では南海スラブの表面や内部で発生した地震の発震機構を利用して右回転力の証拠を探った。その結果紀伊水道南を回転中心とする (円弧)ar1〜ar5 の接線と圧縮軸が調和的な地震が広く分布し右回転を仕向ける力が実在することが判明した。
本稿は (1)(2) を統合することで得られる成果について報告したい。
図 2 のブロック矢印は (2) で得られた南海スラブの各領域で発生した地震の平均的な圧縮軸の方向である。縁を押す
圧縮力はスラブ内を伝達し応力集中する海山やその他に達する。その代表的伝達ルートをブロック矢印に続けて (破線矢印)A2〜D2 で表現したが A2 以外は海山に直結する「北から」圧縮ルートである (図 2)。これらは 2004 年に劇的に方向転換した。A2,B2 がそれぞれ A1,B1 に変化し「北西から」圧縮ルートを形成した (図 1)。A1 の根拠については 1944 年の震源と 3 個の余震が列を成して海山 WM と結ばれている (図 1) ので 1944 年以前の圧縮ルートと判断できたことによる。
浜名湖より南で ar4 と ar5 の間の範囲では 1944 年以降 2004 年まで中小地震の発生も少ない (3)。これは 1944 年の破壊により応力伝達ルートとして失格になった証と考えたい。1944 年での右回転の外周は ar5 より内側の亀裂 (b) をなぞった (4)。1944 年までの伝達ルートは C1,E1 のようであったろうから横ずれと潰れが発生した (図 1)。ただし浜名湖の真下は両スラブの分離・密着の始点終点で履歴により形状複雑なところへまた南海スラブが這い上がる。新たな破断・褶曲が避けられない。1944 年の地震前後で観測された上下変動 (天竜川東岸の隆起と名古屋周辺の沈降) の説明に浜名湖周辺域での前駆的すべりを仮定する (5) よりは本震破壊過程 (根幹はスラブの回転・這い上がり) そのものであろうと主張しておきたい。この上下変動とスラブの運動の効果・影響は調和的である。それでは 1944 年以降の伝達ルートはどうなったか。C1 は西へ反れて「北から」圧縮ルート C2.D2 に変化し、E1 は E2 の如く東方へ向きを変えた圧縮力に変化したはずである (図 2)。
(6) は非定常地殻変動 (東海スロースリップ) の水平変位を以下如く明らかにした。観測開始から 2000 年半ばまでは定常と見なせるがそれ以降 2004 年 8 月まで南東方向への変位が連続した。地震後は 2005 年を経て 2006 年までに南南東方向への変位に変わった。この事実は、2004 年の地震の前の E2 で表す圧縮力の存在と矛盾せず、さらにその地震後 E2 はE1 で表す圧縮力に変化したとすることに調和的である。
今後の見通しは「北西から」圧縮がしばらく続いた後「北から」圧縮に戻す地震が発生する。これが次回の東南海地震である。2004 年の地震は折り返し点で文字通りの中間点なら今回のサイクルは 120 年、残された時間は 50 年である。
(1) 間瀬博文/杭の如き応力集中を起こす南海トラフ沿い２つの海山/JpGU2015S-SS30 投稿中(2) 間瀬博文/発震機構が示す紀伊半島とその下のスラブの右回転傾向の証拠/JpGU2015S-SS31 投稿中(3)AIST/地下構造可視化システム/気象庁震源カタログ (4) 間瀬博文 (2014)/Jpgu2014/SSS29-P10(5) 木股・鷺谷 (2005)/水準測量データの再検討による 1994 年東南海地震プレスリップ/名大(6) 水藤・小沢 (2006)/12-1GPS 連続観測から見た東海スロースリップ/GSI

1/2

SSS27-P01

会場:コンベンションホール

2/2

時間:5 月 24 日 18:15-19:30

SSS27-P02

会場:コンベンションホール

時間:5 月 24 日 18:15-19:30

電界観測による地震予知
Predicting earthquakes just before by observing electric fields
高橋 耕三 1∗
TAKAHASHI, Kozo1∗
1
1

無し
None

1. 地震前兆電界
地震前兆電界が下記のメカニズムで発生する可能性が大きい。
(1) 震源域の地殻の結晶境界面に地震前に微小亀裂が発生し、間隙水が浸透する。
(2) この間隙水に結晶境界面のウラン化合物・ラジウム化合物・ラドンなどの放射性物質が溶け出す。
(3) 亀裂により間隙水が地表と繋がり、放射性物質を含む水が震源域の地表に噴出する。
(4) 地表に出た放射性物質により、下層大気が電離し、電気伝導度が局所的に増加する。
(5) 地表と電離層間の数百万 Volt の電位差と電気伝導度の局所的増加により、宇宙線シャワーの軌跡に沿った電流が、震源域で局所的に一時的に増大する。
(6) この電流は間歇的脈流であり、広帯域の地震前兆電波が輻射される。
上記のメカニズム成立の基本となる (1) の微小亀裂の発生は観測されていないが、(3) の放射性物質及び（５）の電流による地震雲及び (6) の広帯域電波は観測されており、地震前兆の微小亀裂発生の可能性及び上記のメカニズムの妥当性は大きい。
2. 上記（５）の地表・電離層間電流発生のメカニズム
積乱雲の雲頂は約 100 MV に達し、電離層にたいし高電圧となり、電離層から電子・負イオンが雲頂に流れ込み、電離層と地表の電位差は数百万 Volt となる。雲頂は下記のメカニズムで高電圧になると推定される（付図参照）。
(I) 中緯度では、積乱雲の -10 ℃以下の領域で水滴は氷晶となっている。固体の融点は表面の方が内部よりも低い。氷晶の場合 -10 ℃付近では氷晶は水膜で覆われている。氷晶内には正孔 (ph) と自由電子 (fe) が同数あり、ph は結晶から抜けだせないが fe は水膜に移動できるため、負に帯電した水膜で氷晶は覆われる。積乱雲の中では氷晶は激しく衝突し合っており、-10 ℃以下では、衝突の多くは弾性的である。衝突の際の速度変化は小さい氷晶の方が大きい氷晶のそれよりも大きい。そのため、氷晶間の衝突により、小さい氷晶の水膜は大きい氷晶に移り、小さい氷晶は、小さくなるとともに正に帯電し、上昇気流により雲頂に運ばれ、雲頂は高電圧となる。
一方、大きい氷晶は、大きくなるとともに負に帯電し、地表に落下する。
(II) 低緯度では水滴は凍らないため、メカニズムは異なる。積乱雲内部の電界は、上向きで約 1000 ｖ/m になっている。
このため、水滴の上部は負に、下部は正に偏極する。積乱雲内では上昇気流のため、小さい水滴は大きい水滴よりも高速になっており、小さい水滴は、大きい水滴の下部に衝突し、小さい水滴の上部の負の電荷が大きい水滴の下部の正の電荷と中和し、小さい水滴は正に帯電し、中緯度の場合同様に、雲頂を高電圧にする。
(III) 火山の噴煙でも雷の発生が観測されている。噴煙の場合は、火山噴出物（火山灰・礫・岩塊）の摩擦電気による帯電であり、帯電した噴出物は上記の（II) と同じ理由により、噴煙の上部を正の高電圧にする。この高電圧は観測されており、上記のメカニズム（Ｉ）・（II）が共に妥当なことを示している。
３. 前兆電界による地震予知（Refer to (1),(2)）
震源域は微小亀裂の発生領域と一致すると考えられるため、前兆電界の発生領域が震源域となる。規模は震源域の広さにほぼ比例するため、前兆電界の発生領域の広さから推定できる。発生日時は、巨大地震の場合、前兆電界が観測され始めてから約一週間後に発生している。
References(1) ”Detectable magnitude of earthquakes and epicentral distance with precursor at 1.5kHz, Kozo Takahashi, et al. Japan
Geoscience Union Meeting 2007, S148-017(2) ”Short term prediction of shallow great earthquakes by observing electric fields” Kozo Takahashi, et al. Japan Geoscience
Union Meeting 2009, S151-005

キーワード: 地震予知, 地震前兆電界, 雷発生のメカニズム, 中緯度の雷, 低緯度の雷, 噴煙中の雷
Keywords: earthquake prediction, precursory seismic slectric fields, mechanism of generating thunder, thunder in middle-

1/2

SSS27-P02

会場:コンベンションホール

latitude, thunder in low-latitude, thunder in smoke of volcano

2/2

時間:5 月 24 日 18:15-19:30

SSS27-P03

会場:コンベンションホール

時間:5 月 24 日 18:15-19:30

繰り返し回数が非常に少ない地震系列に対する長期発生予測の成績
Long term predictability for the earthquake recurring a few times
田中 昌之 1∗ ; 岡田 正実 1 ; 内田 直希 2
TANAKA, Masayuki1∗ ; OKADA, Masami1 ; UCHIDA, Naoki2
1
1

気象研究所, 2 東北大学大学院 理学研究科
Meteorological Research Institute, 2 Graduate School of Science Tohoku University

１．はじめに相似地震（小繰り返し地震）は、プレート境界の小さな固着域の破壊で発生すると考えられている（Igarashi et al., 2003;
Uchida et al., 2003）。他の地震との識別が容易で、比較的短期間に繰り返すことから、予測実験にも使用されてきた（Okada
et al., 2012）。一方、地震調査委員会が長期評価の対象としている大規模地震や、気象研究所（2014）などが調査した中規模地震では、発生間隔が長く、観測回数が非常に少ないケースが多い。繰り返し回数が非常に少ない地震系列に対する長期的発生確率予測の信頼度を見るために、相似地震を用いて、Bayes 統計などで予測した成績がデータ数によってどのように変化するかを調査・検討する。
２．データと予測方法
相似地震のカタログは、東北大学で作成したもので、筆者達が 2008 年に事前予測の実験に使用したものである。この2008 年予測実験 は、127 系列を対象とし、各系列は 1993 年〜予測時点（2008.1.1）に 5 個以上の相似地震を含み、平均マグニチュードは 2.75 以上であった。今回は、この実験結果に加え、予測時点（2008.1.1）に近い地震を 2 個、3 個、4 個、または 5 個を抽出したものを用いて、2008 年の発生確率を計算し、データ数による成績変化を調べる。
予測方法は、ベイズ統計対数正規分布モデル（LN-Bayes）、小標本論対数正規分布モデル（LN-SST）及び指数分布モデル（Exp-Pin）である（Okada et al., 2012）。指数分布モデルは、データから得られた平均発生間隔をそのまま発生間隔
分布のパラメータとして使用する Plug-in 方式である。
３．予測成績
予測成績は、2008 年に観測された相似地震と比較して求める。平均対数尤度（Mean log-likelihood）および Brier スコ
アは表のとおりである。平均対数尤度は大きい方が、Brier スコアは小さい方が優れた予測である。表中の all は 2008
年予測実験の結果であり、他は予測（2008.1.1）直前のイベント（2〜5 個）のみから発生確率を計算した場合である。
LN-SST では、発生間隔のデータが 1 個（イベント 2 個）のみの場合に発生確率の計算ができない。
いずれの予測方法でも 、予測に使用する地震が少なくなるにつれて、両方のスコア（成績）が悪くなる。しかし、成績悪化はさほど急速ではない。予測方法を成績のよいものから並べると、LN-SST、LN-Bayes、Exp-Pin である。一方、2006 年から 2010 年に実施した 4 実験全体（528 予測）の成績では、LN-SST と LN-Bayes は同程度であり、優劣はほとんどないと考えられる。Exp-Pin は、成績が最も劣るのが通例である。
なお、データ期間、予測期間、及び平均マグニチュードなどの条件を変更した場合の成績変化についても調査を進めており、連合大会当日にはそれらの結果も紹介する予定である。

Table 1. Number of events and the forecast scores for 2008.
—————————————————————————————————————
Mean log-likelihood
Brier score————— ——————————————— ———————————————
events
LN-Bayes LN-SST Exp-Pin
LN-Bayes LN-SST Exp-Pin————— ——————————————— ———————————————
all-0.531-0.497-0.638
0.178
0.167
0.223
5
-0.546-0.511-0.644
0.182
0.174
0.226
4
-0.566-0.521-0.658
0.190
0.176
0.233
3
-0.611-0.556-0.675
0.204
0.184
0.240
2
-0.671−-0.809
0.228
−
0.254
—————————————————————————————————————
・平均対数尤度と Brier スコア
予測発生確率を P とし、該当するイベントが発生したときは Ev=1、発生しなかったときは Ev=0 とする。平均対数尤
度は Ev∗log(P)+(1-Ev)∗log(1-P) の平均 で、Brier スコアは (P-Ev)ˆ2 の平均 である。 log は自然対数。
キーワード: 繰り返し地震, 発生予測, ベイズ定理, 小標本論, 平均対数尤度, Brier スコア
Keywords: recurrent earthquake, forecast, Bayesian approach, Small sample theory, Mean log-likelihood, Brier score

1/1

SSS27-P04

会場:コンベンションホール

時間:5 月 24 日 18:15-19:30

地震活動静穏化・活発化域と大地震の震源域
Seismic quiescence and activation that precede large earthquakes and their source regions
吉川 澄夫 1∗ ; 林元 直樹 2 ; 明田川 保 3
YOSHIKAWA, Sumio1∗ ; HAYASHIMOTO, Naoki2 ; AKETAGAWA, Tamotsu3
1
1

気象研, 2 気象研, 3 沖縄気象台
MRI, 2 MRI, 3 OMO

筆者らは近年開発した地震活動静穏化・活発化解析手法（eMAP）を用い大地震に先行する静穏化現象の検出と定量化を試みている．これまでに静穏化域の大きさと先行時間の双方に地震規模に関するスケーリング則が認められたことか
ら中期的な地震発生予測の可能性を示した（例えば吉川・他，2014）．今後実用的見地から精度向上が必要となることから，今回は地震活動の静穏化現象の発生場所・先行時間と本震の震源域との関係を詳細に評価していくこととする．
調査には従来と同様気象庁震源カタログ（1983 年?2011 年）からマグニチュード７クラス（M6.7 以上）以上で震度 5弱以上を観測した地震（深さ 200km 以浅）など全 26 個を対象に選んだ．これまでの調査の結果では 15 個の地震について本震発生前の地震活動の静穏化が認められており，相対震源距離（破壊開始点から静穏化域中心までの距離）と静穏化域の大きさ（長径）との間に比例関係があることが明らかになっている． ここではこれら 15 個の静穏化検出事例における調査結果を示す．この結果によると静穏化域を活発化域が一部を取り囲む，もしくは接するドーナツパターン(Mogi,1969) が 12 例ほど見受けられた．これは活発化域が出現することにより静穏化域を特定できる場合が多いことを反映するものである．いずれの場合も事前の予想に反して静穏化域がそのまま将来の震源域になるということは殆どな
く，その大きさや場所を決定づけるものではなかった．しかしドーナツパターンを形成する場合，活発化域に震源域が重複する場合が少なくとも 9 例見られた．これはドーナツパターンが将来の震源域を予測するのに貢献する可能性を示す．
地震空白域は次の２通りに分類される（Mogi,1979）．すなわち第１種地震空白域は一連の大規模地震帯の中で未破壊域もしくはアスペリティとして残された地域であるのに対し，第２種地震空白域は地震活動の一時的静穏化の出現した地域である．静穏化域が将来の震源域と一致しないということは，静穏化現象が未破壊域やアスペリティのような高強度域で起こるのではなく非地震性滑りなどが生じやすい低強度域で起きることを示す．一方で同時期に出現するドーナツパターンにおいて静穏化域周囲で発生する活発化域は周辺のアスペリティにおける地震活動を反映するものと考えられる．

キーワード: 地震活動, 静穏化現象, 地震予測
Keywords: seismicity, quiescence, earthquake forecast

1/1

SSS27-P05

会場:コンベンションホール

時間:5 月 24 日 18:15-19:30

1946 年南海地震前の紀伊半島沿岸部における目撃証言
On the witness testimonies before the 1946 Nankai earthquake on the Pacific coast of Kii
peninsula, Japan
梅田 康弘 1∗ ; 板場 智史 1
UMEDA, Yasuhiro1∗ ; ITABA, Satoshi1
11

産総研
GSJ, AIST

１．調査の目的
1946 年南海地震（以下では本震と呼ぶ）の前に起こったことを可能な限り詳しく調べている．これまでに四国沿岸部においては，本震直前までの地殻の上下変動，地下水（主に井戸水）や海水位の変化に関する目撃証言を収集してきた．今回は紀伊半島における本震前の異常現象に関する目撃証言の調査結果を報告する．また四国における調査結果とも比較議論する．
２．調査の方法
文献調査が主である．紀伊半島の全体的な報告としては，本震直後に調査が行われた以下の 3 点による．
A．中央気象台（1947）「 昭和 21 年 12 月 21 日南海道大地震調査概報」．
B．水路局（1948）水路要報 201 号「地変・被害編」
C．同胞援護会和歌山県支部（1948）「昭和紀伊洪浪の記」
その他は各地における体験談などの資料による．これらを収集し，まとめた資料としては「D．和歌山県災害史，和歌山県（1963）」など，他にもいくつかある．証言は，地盤や海面に関するものの他，発光現象，生物と漁獲，天候など多岐にわたるが，本稿では地殻の上下変動に関連すると思われる現象として，井戸水の変化，地盤・海面変動などを取り上げる．
３．井戸水の減少
本震前に井戸水が減少したのは 7 カ所で 9 件の証言がある．東海岸の尾鷲から串本，紀伊由良を経て淡路由良まで分布している．水位低下の証言数は四国の 18 カ所に比べると多くは無いが，地震前に水位が低下した井戸があったことは確かなようである．逆に，四国ではほとんどなかった井戸水の濁りは紀伊半島では 7 カ所ある．本震前に水位が上昇したという証言は，四国の場合と同様に無い．
４．地盤・海面変動など
鳥羽市から那智勝浦にかけて地盤沈下が報告されているが，本震前の現象か，2 年前の 1944 年東南海地震あるいはその余効変動によるものかはわからない．串本では「地盤が次第に沈下し東海岸にあっては砂浜が殆ど無くなり，満潮のときは波が県道（現在の国道 42 号）まで達していた」との証言もある．
海面変動に関するは 3 か所あるが，四国での 38 カ所に比べるとわずかである．目撃されるような大きな海水位の変化は無かったらしい．海鳴りや振動（地震か）も 2 か所で 3 件の証言があるが，東南海地震の余震とも考えられる．
５．異常なしについて
海南市役所と塩屋村役場（現在の御坊市塩屋）では「異常らしいものを認めなかった」と証言されている．特に塩屋村役場では「・・・強震の起こる前には地響き海鳴りなぞの徴候があると昔から言ひ伝へられてゐるがか様な前兆は全然感じられず，・・・」とあり，地震前の異常に関する言い伝えがあることは興味深い．
６．議論と結果
紀伊半島全体としては四国に比べて本震前の異常現象に関する証言数は少ない．高知県の須崎湾や宇佐港で目撃された 2m〜3m もの海水位低下は紀伊半島ではない．直前の異常現象はあったとしても僅かだったか，気付かれない程度だったと思われる．
地盤地下に関しては，水路要報 201 号「地変・被害編」によれば，紀伊半島の東海岸では 1944 年東南海地震時に0.6m 程度の沈降があったと報告されている．小林・ほか（2002）は那智勝浦町浦神の潮位記録を解析し，本震の 2，3 日
前から 10cm 程度の沈降があった可能性を指摘している．4 節に記した串本の証言「地盤が次第に沈下し，・・・」は，東南海地震時あるいはそれ以降の余効変動なのか，それとは別に本震直前に起きた現象なのかはわからないが，小林・ほかの指摘を考えると，直前の地盤沈下も含まれている可能性がある．
四国で井戸水が減少した領域は，本震前の隆起域に対応していたが，串本での証言は，当地が本震前の沈下域に相当することから，四国の場合とは逆である．しかし，すぐ近くの潮岬では伝聞として土地の隆起も証言されており，さらに精査が必要である．

キーワード: 1946 年南海地震, 目撃証言, 井戸水, 海水位変化, 地盤沈下

1/2

SSS27-P05

会場:コンベンションホール

時間:5 月 24 日 18:15-19:30

Keywords: 1946 Nankai earthquake, witness testimony, well water, sea level change, land subsidence

2/2

SSS27-P06

会場:コンベンションホール

時間:5 月 24 日 18:15-19:30

トンガ・ケルマディック海溝で発生するプレート境界型地震と潮汐との関係
Relation between tidal triggering effect and interplate seismicity along the Tonga-Kermadec
trench
弘瀬 冬樹 1∗ ; 上垣内 修 2 ; 前田 憲二 1
HIROSE, Fuyuki1∗ ; KAMIGAICHI, Osamu2 ; MAEDA, Kenji11
1

気象研究所, 2 気象庁
Meteorological Research Institute, 2 Japan Meteorological Agency

トンガ・ケルマディック海溝は世界でも有数の地震多発帯である．プレート収束速度は南から北にかけて増加し，それに比例して地震発生率は高くなっている［Ide, 2013, Nature Geo.］．この地域では，Mw7.5 以上の大規模な地震活動（プレート境界型以外の地震も含む）は 1980 年前後と 2010 年前後に比較的集中している．前者の期間については，Tanaka, etal., [2002, GRL] が 1977 年〜2000 年の GCMT 解データを用いて，1982 年 12 月トンガ地震（Mw7.5）前に p 値（Schuster,
1897, PRSL，地震活動と潮汐との相関を表す指標で 0-1 の間の値を取る．通例として 0.05 以下であれば有意に相関が高いと判断される）が低下し，本震後に増加していることを指摘した．後者の期間については，2009 年 3 月に同規模（Mw7.6）のプレート境界型地震が 1982 年の本震付近で発生している．
本研究では，先行研究［Tanaka, et al., 2002, GRL］で用いられたデータの期間を延長し，2009 年本震前後の p 値の時間変化が 1982 年と同様の傾向を示すか調査した．用いたデータは GCMT 解のプレート境界型地震（すべり角 60-120°，
深さ 70 km 以浅，走向 150-230°，1977 年〜2013 年）である．地殻の理論潮汐応答は，固体地球潮汐と海洋潮汐荷重効
果の和で表現され，前者は earthtide mod［小沢, 1974, 測地学会誌; 中井, 1979, 緯度観測所彙報 上垣内, 2015, 私信］を，
後者は Gotic2 [Matsumoto, et al., 2001], をベースに改造したプログラム［Kamigaichi, 1998, PMG; 上垣内, 2015, 私信］をそれぞれ用いて，各イベントの位置における歪テンソル 6 成分を算出した．このとき，グリーン関数の計算に用いる地球モデルは PREM とした［上垣内, 2015, 私信］．歪テンソル 6 成分から断層面上のΔ CFF（摩擦係数は 0.4 とした）を算
出し，イベント発生時刻の位相を決定した．イベント 50 個ずつを単位とし，1 イベントずつずらして p 値の時間変化を算出した．結果は以下の通りである．
A. 1982 年 12 月本震前に徐々に低下し本震後に増加
B. 2009 年 3 月本震前に徐々に低下し（ただし，最小で 0.1），本震後に増加
C. p 値が 0.05 以下となった時期は全部で 5 回（1982 年 12 月，1988 年 1 月，1993 年 6 月，1998 年 4 月，2000 年 8 月）
あるが，1982 年 12 月を以外の 4 回については，対応する大きめの地震（Mw7.0 以上）はない．
p 値の時間変化については，2004 年スマトラ沖地震 Mw9.0（とその最大余震 Mw8.6）や 2011 年東北地方太平洋沖地
震 Mw9.0 の前に低下し，その後増加するという傾向がみられることから［Tanaka, et al., 2010, 2012, GRL］
，地震予測に有効なツールとなると期待されている．しかしながら，p 値を用いた予測を行う際は上記 C のように空振りとなることも考慮して慎重に行う必要がある．
キーワード: 固体地球潮汐, 海洋潮汐荷重効果, Δ CFF, p 値, トンガ・ケルマディック海溝
Keywords: Earth tide, Ocean tidal loading effect, delta CFF, p-value, Tonga-Kermadec trench

1/1