談 話 室分場拾遺Ⅸ 柳瀬氏から託されたデータから－地震と漁獲－Ⅰ 寺田寅彦
柳瀬良介氏は平成 13 年度、伊豆分場長で定年となり、静岡県職員を退職しました 。その際、柳瀬氏は「伊東分場に勤務していた時に気になって収集した」と言って、私にあるデータを託されました。
伊豆半島では昭和 49 年の伊豆半島沖地震以降、主に伊東周辺で群発地震が発生し、平成元年の手石海丘噴火につながり、その後も群発地震が発生していました。データ
は定置網の日々の漁獲量変動に群発地震の影響が表れていないか、あるいは群発地震の前兆現象が日々の漁獲量変動に表れていないかとの観点から収集されたものでした。図 1 は平成 5 年 1 月に発生した群発地震前後の川奈定置網の魚種毎の日々の漁獲量変動です。破線で地震の日別発生回数を、実線でそれぞれの魚種の日々漁獲量を表しています。サンマでは地震回数と同期しているように入網しています。それ
図 1 柳瀬良介氏の資料に対し、イボダイやキントキ
ダイ、ヤリイカ、ムツ、タチウオといった魚種では群発地震の前にまとまって獲れ1

1

柳瀬良介氏は昭和 53 年度に宮城県水産試験場から静岡県水産試験場に移られ、主に伊豆分場、栽培漁業センターで栽培漁業研究を担当されていました。平成 13 年度をもって県を定年退職
し、神奈川県大磯の自宅で過ごしていましたが、残念なことに病に侵され平成 19 年にご逝去されました。

－17－

ており、群発地震の前兆であるかのようです。しかし、大半の魚種では関係はないようです。柳瀬氏も地震と定置網の漁獲量に明確な関係を示すことができないために、発表する機会がなかったのでは思われます。
伊豆半島では昭和 49 年以降だけでなく、戦前にも群発地震が発生しており、昔から特定の魚種の漁獲量が地震と関係しているのではないかと興味を持たれていました。特に当場が収集している定置網の漁獲量データは多種多様な魚種の漁獲量が日々連続して得られるので、地震の前兆や影響を見るにはうってつけのデータです。今号から数回にわたり、伊豆沿岸を舞台に行われた「地震と漁獲」研究を振り返ってみたいと思います。

－18－

相模湾

駿河湾

研究史Ⅰ 地震と漁獲 寺田寅彦 1-3)
最初にこの問題に科学的なメスを入れたのは、科学者・随筆家として有名な寺田寅彦氏でした。
寺田は東京大学教授に兼務とい内浦湾
う形で数々の研究所に所属していました。
「地震と漁獲」研究 1-3)
は、地震研究所所員でありながら、水産講習所の嘱託であった寺田の面目躍如の研究でした。
図 2 伊豆半島北部と北伊豆地震震源
昭和2年に寺田の東京大学の教え子である木村喜之助氏と宇田道隆氏が水産講習所に奉職、二人は競うように伊豆半島をはさむ駿河湾と相模湾で研究を開始し、成果を発表しました。木村は沼津内浦湾を、宇田は相模湾をそれぞれ研究フィールドとし、定置網の日別漁獲データを収集していました。日本水産学会は昭和 7 年に設立され、同年発行された学会誌第1 巻第 1 号のトップを飾ったのは木村の論文 4)で、内浦湾の定置網に入網したクロマグロ、キハダマグロの漁獲資料から得た体重組成の推移で成長を考察したものでした。木村は水産物理談話会報 15 号（1931 未見）に大正 13～昭和 4 年の内浦湾
淡島定置網の様々な魚種の日別漁獲量の図を載せたところ、地震と漁獲の関係について興味を持っており、信頼できる漁獲データを探していた寺田から昭和 5 年のデータを要望されました。
それは昭和 5 年が伊豆半島で地震の多い年だったからです 5)。2 月 10 日から 4 月
10 日に伊東沖で群発地震が発生し、一時鎮静したものの 5 月 8 日から再活性化し、8 月に収束、その後、11 月 7 日から網代の西 10km の地点で群発地震が発生し、11図 3 昭和 5 年における伊豆周辺の群発地震回数と重寺漁場のマアジ漁獲量Proceedings of the Imperial Academy, 8, Terada(1932)からの転載
月 26 日の早朝に本震である北伊豆地震が発生しており、寺田はこの日別地震データとの比較に耐えうる漁獲データを探していました。
寺田の要望に応えた木村が収集した漁獲データのうち、重寺漁場のマアジが地震の発生回数と良い対応がありました（図 31)）
。ピークで見ると 3 期ある群発地震群とマアジの漁獲は同調しており、地震と漁獲は関係あることを想像させます。寺田
はこの結果に気を良くし、大正 13～昭和 4 年のマアジ漁獲データと地震との関係を検討しました。
寺田が採用した科学的メスは次のような手順でした 1-3,6,7)。
① 漁獲量 n を、y=nlog(3/2)/log2 で変換する。この変換によって漁獲量が倍増した時に変換漁獲量ｙは 1.5 倍になる。
② 変換漁獲量と地震回数を旬毎に集計する。1 月は定置網の操業が行われないので、1 年間で 33 旬（N）のデータが得られる。
③ ある年の変換漁獲量（f）と地震回数（s）の極大値（Ｍ）と極小値（ｍ）の回数を数える。Ｍｆ、Ｍｓ、ｍｆ、ｍｓの 4 種類のデータを得る。
④ Ｍｆ、Ｍｓ、ｍｆ、ｍｓから次の確率を計算する。
漁獲の極大と地震の極大が一致する確率 Ｐ１=Ｍｆ・Ｍｓ/Ｎ2
漁獲の極小と地震の極小が一致する確率 Ｐ２=ｍｆ・ｍｓ/Ｎ2
漁獲の極大と地震の極小が一致する確率 Ｐ３=Ｍｆ・ｍｓ/Ｎ2－19－漁獲の極小と地震の極大が一致する確率 Ｐ４=ｍｆ・Ｍｓ/Ｎ2⑤ ある年の変換漁獲量と地震回数から、直接、漁獲の極大と地震の極大が一致する回数を数え、Ｎで割り、実際の確率Ｑ１を求める。漁獲の極小と地震の極小が一致する実際の確率Ｑ２、漁獲の極大と地震の極小が一致する実際の確率Ｑ３、漁獲の極小と地震の極大が一致する実際の確率Ｑ４も同様に求める。
⑥ ＰｉとＱｉを比較する、また、（Ｐ１＋Ｐ２）/（Ｐ３＋Ｐ４）、（Ｑ１＋Ｑ２）/（Ｑ３＋Ｑ４）を算出する。変換漁獲量と地震回数に正の相関があるなら、Ｑ１＞Ｐ１、Ｑ２＞Ｐ２、Ｑ３＜Ｐ３、Ｑ４＜Ｐ４、（Ｑ１＋Ｑ２）/（Ｑ３＋Ｑ４）＞（Ｐ１＋Ｐ２）/（Ｐ３＋Ｐ４）となる。
一読しただけでは
表 1 寺田が採用した方法の解説
理解しにくいものの、
地震
表 1 のように分割表
計
非極大・非で考えると理解しや極大回数 極小回数
極小回数
すくなります。実際
極大回数 Q１⇔P１ Q3⇔P３
Mf/Nの確率Ｑｉと周辺度
漁
極小回数
Q4⇔P ４ Q2⇔P２
ｍf/N
数（Ｍｆ、Ｍｓ、ｍ
獲 非極大・非ｆ、ｍｓ）から期待極小回数される確率Ｐｉとで
計
Ms/N
ｍs/N
1
どちらが大きいか比べているわけです。
この解析の結果、マ表 2 昭和 5 年の独立性の検定結果
アジの漁獲と地震は地震関係があるらしいと1930年
計
極大回数 極小回数結論づけています。
漁 極大回数213
しかし、寺田が採獲 極小回数
21517
用した方法は確率の計
41620
大小を比べているだ
フィッシャーの正確確率検定
P＝0.0877
けでデータ数の情報を使っていません。
従って、統計的な検定ではないことは明らかです。統計的検定を行うとすれば、データは表 1、2 のような分割表にまとめられるので、独立性の検定を行うことが考えられます。独立性の検定 8)とは二つの属性間に関係があるかどうかを検定するもので、この場合は地震と漁獲との間に関係があるかどうかを検定します。実測度数と周辺度数からの期待値との食い違いを利用するχ2 検定、あるいは実測度数が生起する確率を計算するフィッシャーの正確確率検定（期待値が小さい場合）のどち－20－

らかを使用します。
寺田の報告 2)から大正 13～昭和 5 年のデータを復元して、独立性の検定を行ってみました。復元したデータに矛盾が生じることが多く、寺田がどのようにして前掲
手順③と⑤を進めたか不明な点が多々あります。データを当てはめた 2×2 分割表にフィッシャーの正確確率検定 9)を行ったところ、帰無仮説（地震と漁獲は独立＝関係はない）はどの年も有意水準 5％で棄却されませんでした。残念ながら統計的検定ではマアジの漁獲と地震は関係がないとなります。一例として昭和 5 年の検定結果を表 2 に示しましたが、有意水準 10％であれば、帰無仮説は棄却されます。
昭和 5 年のデータは図 3 に示したように、視覚的にはよく対応していましたが、地震と漁獲変動の極大値と極小値の一致に注目した解析では統計的には有意になりませんでした。マアジも含めて定置網に入網する魚種は、地震でなくても移動回遊の結果、ある程度決まった時期に漁獲されます。通常の漁期を超えた地震の影響の抽出方法の検討、並びに採用する統計的検定の検討が必要と思われます。
また、寺田 2)は地震が漁獲に影響を与える可能性を 3 点指摘しました。第一は、地震の衝撃またはそれらに関連する刺激が、魚によって直接感じられ、何らかの方法で、漁場に近づくように誘導する可能性です。第二は、地震の衝撃は、魚の餌となるプランクトンの生息層に何らかの影響を与える可能性です 。第三は、地震は沿岸の海水の化学的性質に影響を及ぼし、地下水の擾乱を引き起こし、プランクトンや魚に間接的に影響を与える可能性です。
文献1) Terada T. (1932) Earthquakes and Fisheries, Proceedings of Imperial
Academy、8、83-86.
2) Terada T. (1932) On some Probable Influence of Earthquakes upon Fisheries、
Bulletin Earthquake Research Institute、10、393-401.
3) Terada T. (1933) Earthquakes and Fisheries Ⅱ、 Bulletin Earthquake
Research Institute、11、713-716.
4) 木村喜之助（1932）重寺漁場の漁獲に表はれしクロシビ、キハダの成長度、日本水産学会誌、1（1）、1-4.
5) https://ja.wikipedia.org/wiki/北伊豆地震
6) Kimura K. (1936) A New Coefficient of Association between Two Attributes、
日本水産学会誌、5（3）、146-148．
7) 宇田道隆（1960）海洋漁場学、恒星社厚生閣、東京、347pp.
8) http://aoki2.si.gunma-u.ac.jp/LaTeX/ExactTest.pdf.
9) http://aoki2.si.gunma-u.ac.jp/exact/exact.html.
（長谷川雅俊）－21－