水産海洋研究79（3） 141–148，2015

宗谷岬沖潮流カレンダーによる樽流し漁業活動の変化

Bull. Jpn. Soc. Fish. Oceanogr.

宗谷岬沖潮流カレンダーによるミズダコ樽流し漁業活動の変化
佐野稔 1 †，坂東忠男 2，本前伸一 3，江淵直人 4

Change in drift-ﬁshing activities due to the introduction of a tidal current
calendar in a coastal area oﬀ Cape Soya in Hokkaido, Japan
Minoru SANO1 † , Tadao BANDO2, Shinichi MOTOMAE3 and Naoto EBUCHI3宗谷岬沖潮流カレンダーは 2009 年から北海道宗谷地区で配布され，ミズダコ樽流し漁業で活用されている．本研究では，
本カレンダーで用いている「樽流しに適した表層流速度」を検証し，同時に本カレンダーの導入による漁業活動の変化を明らかにすることを目的とした．宗谷地区の短波海洋レーダの観測値と樽流し漁業の 1 日あたり 1 隻あたり漁獲量を用いた回帰樹解析の結果，宗谷岬沖の海域において 0.1–1.0 kt の表層流の発生時間が長くなれば操業時間も長くなり，漁獲量も増加した．2007 年の表層流速は，2003–2006 年のレーダのデータから推定した潮汐流と平均流をもとに推定し，この予測値と 2007 年の実測値の間の相関係数は 0.79 を示した．樽流し出漁日数が年間 100 日以上と 50 日以上の漁業者の出漁日数は，カレンダー導入後減少し，樽流しに適さない日に出漁しなくなったものと推定した．

A tidal current calendar for the coastal area off Cape Soya in Hokkaido, Japan, has been used by ﬁshermen of the driftﬁshery for giant Paciﬁc octopus, Enteroctopus doﬂeini, since 2009. In the present study, the accuracy of suitable surface current velocities described in the calendar was veriﬁed and changes in drift-ﬁshing activities resulting from theintroduction of the calendar were clariﬁed. A regression tree analysis of surface velocity measurements recorded byhigh-frequency ocean radar stations in the Soya region and the catch（kg・boat­1・day­1）showed that the catch increased when the surface velocity was in the range of 0.1 to 1.0 kt. Surface velocities were simulated based on tidal
currents and the mean velocity of the Soya Warm Current. A correlation coefﬁcient of 0.79 was obtained between theobserved velocity and the simulated velocity. After the introduction of the calendar, the proportion of actual ﬁshingdays to possible ﬁshing days（Fishing Effort Index）for ﬁshermen who ﬁshed ≥50 days per year decreased. Most ﬁshermen did not go ﬁshing on days considered unsuitable for drift ﬁshing after the introduction of the calendar.
Key words: drift ﬁshery, Enteroctopus doﬂeini, HF radar, Soya/La Pérouse Strait, surface current, tidal current calendar

はじめに
宗谷海峡は宗谷岬とサハリン最南端のクリリオン岬との間にあり，対馬暖流から派生した宗谷暖流が日本海から北海2014 年 10 月 10 日受付，2015 年 7 月 13 日受理
1
北海道立総合研究機構稚内水産試験場
Wakkanai Fisheries Research Institute, Hokkaido Research
Organization, Wakkanai, Hokkaido 097–0001, Japan2
宗谷漁業協同組合
Soya Fisheries Cooperative Association, Wakkanai, Hokkaido 098–
6755, Japan3
根室地区水産技術普及指導所標津支所
Nemuro Fisheries Extension Ofﬁce, Shibetsu, Hokkaido 086–1364,Japan
4
北海道大学低温科学研究所
Institute of Low Temperature Science, Hokkaido University, Sapporo,
Hokkaido 060–0819, Japan†sano-minoru@hro.or.jp

道沿岸のオホーツク海へと流れている．この海域はミズダ
コ Enteroctopus doﬂeini の好 漁 場であり，年間 1,600–4,200 t漁獲される．この海域においてミズダコは主に樽流し（北海道水産業改良普及職員協議会，2006）で漁獲され，その年間漁獲量は宗谷漁業協同組合のミズダコ漁獲量の約 8 割を占める．樽流しでは，浮き（樽）に糸と擬餌針のついた
1 組の漁具を 1 隻あたり15–20 個用いる．これら漁具を潮の流れにまかせて浮標させ，擬餌針が海底を引きずられるよ
うにすることで，ミズダコが擬餌針に飛びつき漁獲される．
そのため，潮の流れが漁獲の成否を決める重要な要素であり，流れが速すぎると擬餌針が海底から離れてしまい，遅すぎると擬餌針がほとんど動かず，ミズダコが漁獲されなくなる．佐野ほか（2012）は，樽流し漁船の航跡の GPS 情報（緯度，経度，時刻）から算出した樽流し漁具の漂流速度とミズダコ漁獲量の関係を解析し，宗谷海峡での最適な— 141 —

佐野稔，坂東忠男，本前伸一，江淵直人

Figure 1. Example of the tidal current calendar for the coastal area off Cape Soya.（a）; calendar month showing thedata for June 2014.（b）; enlarged view of a day on the calendar. Velocity category shows the potential suitability
for drift ﬁshing for giant Paciﬁc octopus on a given day.

漁具漂流速度は 0.1–1.0 kt であり，漁獲量を増やすためには
樽流し時間を増やすことが必要であることを明らかにした．
一方で，宗谷海峡では北海道大学低温科学研究所の短波
海洋レーダによる宗谷暖流の観測が 2003 年より行われている（Ebuchi et al., 2006）．稚内水産試験場，北海道大学低温科学研究所，および稚内地区水産技術普及指導所は共同で，この観測データを用いてミズダコの漁場と重なる北緯45 度 35 分から 45 度 40 分，東経 141 度 50 分から 142 度 10 分の範囲の海域について 1 時間ごとの表層流の発生予測を行い，予測結果を記した宗谷岬沖潮流カレンダー（Fig. 1,
http://www.fishexp.hro.or.jp/cont/wakkanai/section/zoushoku/

f1hig400000002l5.html, 2014 年 8 月 19 日）を 2009 年 か ら 宗
谷漁業協同組合へ提供している．このカレンダーでは，樽流し漁具の漂流速度の解析結果（佐野ほか，2012）をもとに，樽流し漁業において遅い表層流速度を 0–0.2 kt，最適な速度を 0.2–1.0 kt，やや速い速度を 1.0–1.3 kt，速い速度
を 1.3 kt 以上とした．ここで，最適な表層流速度を最適な漁具漂流速度の 0.1–1.0 kt ではなく，0.2–1.0 kt としたのは樽流し漁具の漂流速度が遅くなる時間帯を強調するためである．このカレンダーは宗谷漁業協同組合のカレンダーとして全組合員に配布されている．樽流し漁業者から潮流予測が当たるという好評を得ており，宗谷海峡の樽流し漁業— 142 —

宗谷岬沖潮流カレンダーによる樽流し漁業活動の変化に貢献している（北海道新聞，2011 年 6 月 1 日）．
この宗谷岬沖潮流カレンダーの表層流速度は短波レーダの流速データを用いているが，最適としている表層流の速度は，漁船の GPS データから算出した樽流し漁具の最適な漂流速度（佐野ほか，2012）を根拠としているため，短波レーダによる表層流速度の実用性とその予測値の科学的根拠を明確にする必要がある．そして，カレンダーの提供から5 年が経過しており，宗谷地区の新しい暦として定着している（北海道新聞，2011 年 6 月3日）にもかかわらず，カレンダーの導入による樽流し漁業活動の変化は定量的には解析されていない．これらの検証を行うことは，今後他の海域でも表層流予測情報の提供による漁業活動の効率化を進めていくうえで参考になると思われる．そこで，本研究では最適な表層流速度を検証し，表層流速度予測値の精度を明示するとともに，宗谷岬沖潮流カレンダーの提供による樽流し
漁業活動の変化を定量的に明らかにすることを目的とした．
材料と方法
海洋レーダ
北海道大学低温科学研究所は，宗谷海峡周辺域および雄武­紋別沖海域沿岸に 5 局の海洋レーダ局を設置し，2003 年
8 月より当海域における海洋表層流速の毎時観測を継続して行っている（Ebuchi et al., 2006）．海洋レーダは，陸上に設置したアンテナから短波帯の電波を海面に向けて発射
し，海面からの後方散乱される信号を受信するものである．受信電波のドップラーシフトを通して，海面直下の流速の視線方向成分が観測できる（例えば，土木学会海岸工学委員会，2001）．同一の海面を複数の方向から観測した
視線方向速度をベクトル合成することにより，流速・流向を 求 め る こ と が 可 能 で あ る．宗 谷 海 峡 に お い て は，ノシャップ岬，宗谷大岬，猿払の 3 局を配置することにより，宗谷海峡からオホーツク側の海域の表層流を，空間分解
能 3 km で毎時観測している（Fig. 2, http://wwwoc.lowtem.
hokudai.ac.jp/hf-radar/, 2014 年 9 月 3 日）
．Ebuchi et al.（2006）
は，宗谷海峡域の海洋レーダで観測した表層流速が，漂流ブイ，船舶搭載超音波流速計データと残差の標準偏差20 cm・s­1（0.39 kt）程度で一致することを示している．
樽流しに最適な表層流速度の検証
宗谷海峡において，樽流し漁業の漁具の漂流速度が 0.1–
1.0 kt の 時 に ミ ズ ダ コ の 漁 獲 量 が 多 く な る（佐 野 ほ か，
2012）．そこで，0.1–1.0 kt の表層流の 1 日あたりの発生時間が長ければ宗谷漁業協同組合の樽流し漁業の 1 日あたり
1 隻あたりの漁獲量（kg・boat­1・day­1）も増えるという仮説を立てて，宗谷海峡の 0.1–1.0 kt の表層流発生時間数とこの漁獲量との関係を検討した．宗谷岬沖潮流カレン
ダーを 2009 年から漁業者へ提供しているため，カレンダーが提供される前の 2007 年のデータを用いて解析した．解析に用いた表層流速の観測値は，樽流し漁場と重なる北緯

Figure 2. Main ﬁshing grounds of giant Paciﬁc octopus withHigh Frequency（HF）radar station coverage（NS; Noshappu radar station, SY; Soya radar station, SR; Sarufutsu radar
station）and the area of simulated tidal current in Soya/LaPérouse Strait. Numerals on bathymetric contours indicate
depth in meter.

45 度 35 分から 45 度 40 分，東経 141 度 50 分から 142 度 10 分の範囲（Fig. 2）の平均値を用いた．樽流し漁業の出漁は日の出から日没までなので，4–16 時の間に 0.1–1.0 kt の表層流速が観測された合計時間数を日ごとに算出した．樽流
し漁業の日別ミズダコ漁獲量，日別出漁隻数の集計には，
宗谷漁業協同組合の水揚げ伝票を用いた．出漁隻数は，その日に出荷記録のあった漁業者数とした．1 隻あたり 1 日あたり漁獲量は，宗谷漁業協同組合全体の日別漁獲量の合計をその日の出漁隻数で除した値とした．
日ごとの樽流し漁業の操業状況と表層流速との関係を解析する際に，漁業協同組合内部の取り決めによる出漁禁止
日や，天候により操業できない時化日を除く必要がある．
宗谷漁業協同組合の樽流し漁業は，毎年 4 月1日から 12 月24日まで行われ，その期間中に宗谷漁業協同組合による出漁禁止日が設定されている．そのため，同組合の資料である「休漁日および沖止票」により組合が設定した出漁禁止日を把握して，解析から除いた．時化による休漁日については明確に定義された資料がないため，推定する必要があ— 143 —

佐野稔，坂東忠男，本前伸一，江淵直人

る．時化の指標には波高が適切であると思われるが，宗谷
海峡の日毎の波高を知ることは容易ではない．そこで，気象庁で観測した 2007 年 4–12 月の宗谷岬における風速の日平
均 （http://www.data.jma.go.jp/obd/stats/etrn/index.php, 2014 年
8 月20日）を用いて時化を判断する指標とした．時化による休漁日と判断する風速範囲を明らかにするために，1日あたりの出漁隻数を被説明変数，風速の日平均を説明変数にした回帰樹解析を行い出漁隻数が最も少なくなる時の風速を明らかにし，この風速が認められた日を推定時化日とした．
0.1–1.0 kt の表層流発生時間数が多くなれば 1 日あたり 1隻あたり漁獲量も高くなるという仮説を検証するために，
2007 年 4 月 1 日から 12 月 24 日で組合が設定した出漁禁止
日，風速に基づく推定時化日を除いた日について，表層流発生時間数を説明変数，1 日あたり 1 隻あたり漁獲量を被説明変数にした回帰樹解析を行った．一連の統計解析には，R 2.10.0（Ihaka and Gentleman, 1996）を用いた．
表層流速度の推算
過去の海洋レーダで観測されたデータをもとに，将来の表層流の予測を試みた．ここでは 2003 年 8 月から 2006 年 12月までのデータを用いて，以下の比較的簡便な手法を用いて 2007 年の表層流速を推算し，その推算値を 2007 年の実際の海洋レーダデータと比較した．宗谷海峡の対象海域における海洋表層の流れを
潮汐流

+

平均流

+

吹送流

の 3 つの成分からなるものと仮定した．
潮汐流は，2003–2006 年に海洋レーダで観測された毎時の表層流（東西成分，南北成分）データに調和解析（柳，1993）を行い，得られた主要 8 分潮（K1, O1, P1, Q1, M2,
S2, N2, K2）の調和定数を使って計算した．この海域では，
日周潮流（主に K1, O1）が卓越しており，潮流成分は最大で 60 cm・s­1（1.2 kt）程度以上の大きさを持つことが知られている（Odamaki, 1994; Ebuchi et al., 2006）．
宗谷暖流を含む平均流には，季節変動のみを考慮した気候値を作成して使用した．2003–2006 年の海洋レーダの毎時の表層流（東西成分，南北成分）データに 25 時間の移動平均をかけて潮汐流成分を取り除いた後，日平均データを作成した．それを各年の同じ日付で平均した後，15 日の移動平均をかけて平滑化して気候値とした．対象海域での宗谷暖流は，9 月から 10 月に最も強くなり，その最大流速は 60 cm・s–1（1.2 kt）に達する．また，12 月から 1 月には，10 cm・s­1（0.2 kt）以下の最小値をとる．宗谷暖流の経年変化は，今回の推算では考慮していない．しかし，Ebuchi et al.（2009）によれば，2003 年以降の海洋レーダの観測データを見る限り，対象海域周辺での経年変動は小さ
く，ほぼ同じ季節変動が繰り返されているとしている．
吹送流は，海上風の数 % の大きさを持つで，秋季から冬季の強風時には，潮汐流や平均流と同程度の大きさを持つことが予想される．しかしながら，海上風速・風向を 1
年以上の期間にわたって予測することは現実的には不可能なので，今回の潮流推算では吹送流成分は無視することとした．ただし，上記の平均流の気候値には，平均的な風系による吹送流成分の気候値が含まれている．
宗谷岬沖潮流カレンダーの提供前と後の樽流し漁業の操業状況の比較
樽流し漁業の操業状況の指標として，漁船ごとの年間出漁日数と 1 日あたり 1 隻あたり漁獲量の年間平均値を用い，
宗谷岬沖潮流カレンダーの提供前と後で比較を行った．
2006–2013 年の宗谷漁業協同組合の水揚げ伝票から，漁船ごとの出漁日，日別漁獲量を集計した．組合が設定した出漁禁止日を解析から除くために，宗谷漁業協同組合の資料である 2006–2013 年の「休漁日および沖止票」から組合が設定した出漁禁止日を把握した．さらに，時化による休漁日を除くために，気象庁で観測した 2006 年 4 月から 2013年 12 月までの宗谷岬における風速の日平均（http://www.
data.jma.go.jp/obd/stats/etrn/index.php, 2014 年 8 月 20 日）を用いて，「樽流しに最適な潮流速度の検証」で推定した時化日と思われる風速にもとづいて時化日を推定した．これ
らを用いて，年間出漁可能日数に対する実際の出漁日数の割合（Fishing days ratio: FD%）を求め，カレンダーの導入によって漁業者が自主的に出漁を見合わせた日数の割合の増加を推定した．式はFD % =

D× 100
d − (M + S )である．ここで D は各漁船の年間出漁日数，d は各年の 4
月 1 日から 12 月 24 日の期間で 267 日間（定数），M は出漁禁止（組合が設定した休漁日）の日数，S は推定時化日数であり，2006–2013 年の各年について漁船ごとに算出した．
宗谷漁業協同組合の樽流し漁船は 200 隻近くあり，漁業者により樽流し漁業への依存度は変わり，各漁船の年間出漁日数は異なる．加えて，この期間に廃業，休業や新たに着業した漁業者も含まれる．そこで，2006–2013 年まで毎年出漁した漁業者を抽出し，専業と思われる漁業者として年間最低出漁日数が 100 日以上の漁業者をグループ A と定義した．次にこれ以外の漁業者について，比較的樽流し漁業への依存度が高い兼業漁業者として年間最低出漁日数が50 日以上の漁業者をグループ B に，残りの漁業者すなわち樽流し漁業への依存度が低い兼業漁業者としてグループ Cと 定 義 し た．各 漁 業 者 に つ い て，カ レ ン ダ ー 導 入 前 の2006–2008 年の FD% の年平均値と，導入後の 2009–2013 年の年平均値の差をもとめ，グループごとにカレンダー導入
前と後で漁業者対応のある 1 標本 t 検定を行った．さらに，
1 日あたり 1 隻あたり漁獲量についても，同様にグループ— 144 —

宗谷岬沖潮流カレンダーによる樽流し漁業活動の変化

Figure 3. （a）Daily number of drift ﬁshing boats afﬁliated with
the Soya Fishery Cooperative Society ﬁshing for giant Pacific octopus and（b）mean wind velocity（m・s­1）at Cape
Soya from April to December 2007. Black bars: numbers ofdrift-ﬁshery boats; Gray bars: days on which a moratorium
on drift ﬁshing was declared by the Soya Fishery Cooperative Society.

Figure 5. （a）Daily catch（kg・boat­1・day­1）for the driftﬁshing boats afﬁliated with the Soya Fishery Cooperative
Society catching giant Paciﬁc octopus, and（b）length oftime（h）that the surface current was in the range of 0.1 to
1.0 kt in the study area in Soya/La Pérouse Strait from April
to December 2007. Black bars: daily catches for the driftﬁshing boats.

Figure 4. Regression tree for daily number of drift-ﬁshingboats afﬁliated with the Soya Fishery Cooperative Society
catching giant Paciﬁc octopus. The explanatory variable ismean wind velocity（m・s­1）at Cape Soya; n indicates thenumber of day in the nodes. Numbers above n values indicate number of drift-ﬁshing boats in the nodes.

Figure 6. Regression tree for daily catch（kg・boat­1・day­1）in the giant Paciﬁc octopus ﬁshery. The explanatory variable
is the length of time（h）that the surface current was in therange of 0.1 to 1.0 kt in the study area of Soya/La Pérouse
Strait from April to December 2007. n indicates the numberof day in the nodes. Numbers above n values show the mean
daily catch for the drift-ﬁshing boats in the nodes.

ごとに 2006–2008 年の年平均値と 2009–2013 年の年平均値の間で差を求め，漁業者対応のある 1 標本 t 検定を行った．

数が少なくなり，平均風速 9.05 m・s­1 以上の日で出漁隻数が平均 2.8 隻となった．そのため，平均風速 9.05 m・s­1以上の日を推定時化日とした．
1 日あたり 1 隻あたり漁獲量が比較的高くなる期間と，宗谷岬沖に 0.1–1.0 kt の表層流が発生する時間数が多い期間が一致する傾向を示した（Fig. 5）．この漁獲量は日による変動はあるものの，4 月から徐々に増加して 9 月ごろ最も高くなり，10 月以降に減少した．9–11 月では 0.1–1.0 kt の表層流発生時間数が 5 時間を超える日で漁獲量が多くなった．
回帰樹分析の結果（Fig. 6）, 0.1–1.0 kt の表層流発生時間数
が 5.5 時間以上であれば，漁獲量は 121.5 kg・boat­1・day­1となり，この発生時間数が短くなれば漁獲量も減少し，0.5

結 果
樽流しに最適な表層流速度の検証
2007 年 4–12 月の出漁禁止日を除いた日について，出漁隻数が 0 の日 の多くは，日平均風速が約 10 m・s­1 以上であった（Fig. 3）．出漁隻数は 4 月では多い日で約 40 隻であり，5 月以降では約 70–80 隻となった．出漁がほとんど認められない日は 9–12 月に頻繁に認められ，そのときの風
速は約 10 m・s­1 以上であった．回帰樹分析の結果（Fig. 4）
では，平均風速が大きくなるにつれて樽流し漁船の出漁隻— 145 —

佐野稔，坂東忠男，本前伸一，江淵直人

Figure 7. Comparison between observed and simulated velocities in the area of the current in Soya/La Pérouse Straitfrom 2003 to 2006. Observed velocities were recorded by the HF radar;（a）observed velocities from 1 May to 1
October in 2007,（b）simulated velocities from 1 May to 1 October in 2007.

時間未満では 27.1 kg・boat­1・day­1 と低かった．
表層流速度の推算
2007 年の実際の観測データ（Fig. 7a）と，2003 年 8 月から2006 年 12 月までの海洋レーダ観測データを使って推算した表層流速（Fig. 7b）を比較した．観測データは，Fig. 2に示した対象海域内での平均値を用いた．観測データは，
時間スケールの短い風の変動による吹送流成分を含むた
め，変動が大きいが，全体的な変動のパターンは推算値とおおむね一致していた．特に，強い流れが現れる時間帯が日付とともに変わっていく様子や，大潮・小潮の 14 日周期で流速が変動している様子などがよく再現された．
2007 年の毎時の東西流速の観測値と推定値の関係を Fig.
8 に示した．吹送流によると思われる観測値の変動が大き
いため，点のばらつきは比較的大きいが，全体的な傾向はよく一致していた．残差の平均値は 8 cm・s­1（0.16 kt），
標準偏差は 28 cm・s­1（0.54 kt），相関係数は 0.79 であった．
宗谷岬沖潮流カレンダーの提供前と後の樽流し漁業の操業状況の比較
樽流しによりミズダコを水揚げした漁船数は，年変動があるもの 163–196 隻の範囲であった（Table 1）．操業許可日
数は，宗谷岬沖潮流カレンダーを提供する前（2006–2008年）の 245–246 日間に比べて，提供後（2009–2013 年）には出漁禁止日（組合が設定した休漁日）が増加したことにより 237–244 日間に短くなった．推定された時化の日数は年によるばらつきがあり，46–66 日の範囲を示し，最大 20日間の差があった．結果として，操業許可日数から時化日
数を除いた出漁可能日数は，カレンダー提供前の 2006–

Figure 8. Relationships between observed and simulated velocities in the area of the current in the Soya/La Pérouse Strait
in 2007. Observed velocities were recorded by the HF radar.

2008 年には 180–200 日間に対し，提供後の 2009–2013 年には 176–185 日間に短くなった．
樽流し漁業への依存度が高いグループ A の FD% の平均
値は，カレンダー導入前で 68.6% であったが，カレンダー
導入後には 64.6% まで有意に減少していた（対応のある 1
標本 t 検定：p=0.04; Fig. 9）．抽出された樽流し漁業者の半数を占め，樽流し漁業への依存度が中程度のグループ B に— 146 —

宗谷岬沖潮流カレンダーによる樽流し漁業活動の変化

Table 1. Number of ﬁshing boats, total catch, and possible ﬁshing days from 2006 to 2013.

Year

No. ofﬁshingboats

Totalcatch（t）

Permitted
ﬁshingdays1）

Stormydays onﬁshingdays2）

Possibleﬁshingdays3）

2006
2007
2008
2009
2010
2011
2012
2013

194
163
178
194
188
196
190
185

1,413
1,407
1,533
1,309
1,329
1,470
1,359
1,425

246
246
245
244
243
238
237
237

46
66
49
64
58
62
61
53

200
180
196
180
185
176
176
184

Figure 10. Difference in catch （kg・boat­1・day­1） before（2006–2008）and after（2009–2013）introduction of the tidal calendar to the drift ﬁshery in the coastal waters off Cape
Soya. Group A（n=6）: ﬁshermen with ≥100 d of actual ﬁshing per year from 2006 to 2013. Group B（n=40）: ﬁshermenwith ≥50 d of actual ﬁshing per year, exclusive of Group A.
Group C（n=32）: ﬁshermen with <50 d of actual ﬁshing peryear from 2006 to 2013. p: the probability of paired-sample
t-test. Error bars show standard deviations.

1）Days exclusive of the moratoria on ﬁshing.
Wind velocity ≥ 9.05 m・s­1 on permitted ﬁshing days.
3）
Wind velocity<9.05 m・s­1 on permitted ﬁshing days.
2）

有 意 差 は な か っ た が（導 入 前：103.4 kg・boat­1・day­1，
導入後：93.7 kg・boat­1・day­1; p=0.06），導入後は導入前の 91% を示した．

Figure 9. Difference in mean ﬁshing days in percentage before（2006–2008）and after（2009–2013）introduction of the tidal calendar to the drift ﬁshery in the coastal waters off Cape
Soya. Fishing days: Proportion of actual ﬁshing days to possible ﬁshing days（%）. Group A（n=6）: ﬁshermen with≥100 d of actual ﬁshing per year from 2006 to 2013. GroupB（n=40）: ﬁshermen with ≥50 d of actual ﬁshing per year,
exclusive of Group A. Group C（n=32）: ﬁshermen with<50 d of actual ﬁshing per year from 2006 to 2013. p: theprobability of paired-sample t-test. Error bars show standard
deviations.

おいても，FD% の平均値はカレンダー導入前の 56.9% から，導入後の 48.7% まで有意に減少した（p<0.001）．グルー
プ C では，カレンダー導入前後で FD% の平均値に有意差は認められなかった（p=0.89）．
1 日あたり 1 隻あたり漁獲量の平均値は，グループ A, B
において宗谷岬沖潮流カレンダー導入後に有意に減少した（それぞれ p=0.04, <0.001; Fig. 10）．平均値はカレンダー導入前にグループ A で 180.1 kg・boat­1・day­1，グループ B
で 158.1 kg・boat­1・day­1 であったが，導入後にはそれぞれ 166.6 kg・boat­1・day­1, 131.1 kg・boat­1・day­1 に 低 下した．グループ C では，カレンダー導入の前後で平均値に考 察
本論文の目的である，比較的簡便な方法によって，ミズダ
コ樽流し漁に最適流速が発生する時間帯を，将来の 1 年以上の期間にわたって予測する観点でいえば，本研究で提案した表層流速推定モデルと実測値との差は，実用的には問題のない精度と考えられた．
宗谷海峡において樽流し漁船の GPS データと漁獲量の関係から算出したミズダコ樽流し漁具の最適な漂流速度
は，0.1–1.0 kt と推定されている（佐野ほか，2012）．本研
究では，短波レーダによる 0.1–1.0 kt の表層流速度帯の発生時間数と樽流しの 1 日あたり 1 隻あたり漁獲量の関係を検証した結果，この発生時間数が長くなると漁獲量が多くなることを示した（Fig. 6）．したがって，短波レーダによるデータにおいても，0.1–1.0 kt の表層流速は，樽流し漁業に適した流速と判断される．宗谷海峡内の表層流速は潮汐の影響により日周期的に変化しており，その範囲は0–2.0 kt で あ る（Ebuchi et al., 2006）. 0.1–1.0 kt の 表 層 流 発
生時間も周期的に増減しており（Fig. 5），宗谷岬沖潮流カレンダーを導入する前には，漁業者は 0.1–1.0 kt の表層流が発生する時間帯を潮汐表と経験などから予測して出漁していたと思われる．
0.1–1.0 kt の表層流が 10 時間以上となる日は，4–5 月で多かったが（Fig. 5），1 日あたり 1 隻あたり漁獲量の最大値
は 4 月でもっとも小さく，その後徐々に大きくなり 9 月に最 大 と な っ た（400 kg・boat­1・day­1 以 上）．そ の た め，
この最大値の季節変化は 0.1–1.0 kt の表層流発生時間だけ— 147 —

佐野稔，坂東忠男，本前伸一，江淵直人

では説明できない．宗谷海峡において樽流し漁場は 4–6 月には宗谷岬周辺の浅所に形成され，7–10 月にかけて深所へ移り，11–12 月には再び宗谷岬周辺の浅所へと変化する（佐野ほか，2012）．一方で，樽流し漁業の漁獲物の体重組
成は 5–6 月に体重 4–5 kg にモードがあり体重の大きい方への裾広がりの分布型であったものが，7–8 月には体重の大きい個体が漁獲されなくなり，9 月以降再び体重の大きい
方 へ 徐 々 に 裾 広 が り と な る 分 布 型 に な る（佐 野 ほ か，
2011）
．つまり，ミズダコの分布密度や，樽流し漁具で漁獲される大きさも季節的に変化することから，1 日あたり
1 隻あたり漁獲量の月別最大値が季節変化するのは，ミズ
ダコ資源の時空間的な動態やミズダコの漁獲サイズの変化も反映しているものと思われる．
本カレンダーの導入により，大半の樽流し漁業者は，0.1–1.0 kt の表層流がほとんど発生しない日の出漁を止めていたと思われる．樽流し漁業への依存度が比較的高い漁業 者 の グ ル ー プ A, B の FD% の 平 均 値 は そ れ ぞ れ 64.6,
48.7% で，導入前に比べ­4.0 および­8.2 ポイントとなった（Fig. 9）. 2009–2013 年の出漁可能日数の平均値は 180.2 日であったことから，カレンダーが未導入だった 2006–2008 年に比べてグループ A では 7.2 日，グループ B では 14.8 日短
くなったこととなる．一方，年間最低出漁日数が 50 日未満のグループ C では，カレンダー導入後も FD% の減少は認められなかった（Fig. 9）．このグループ C は，他の漁業を優先したため，カレンダーの導入による出漁日数の減少の効果が顕著ではなかったものと思われる．
カレンダーの導入により，いずれの漁業者グループでも樽流し漁業の 1 日あたり 1 隻あたり漁獲量は上昇しなかった（Fig. 10）．グループ A, B は，最適な表層流発生時間数が少ない日，すなわち漁獲が期待できない日を避けて出漁していたので，1 日あたり 1 隻あたり漁獲量が少ない日が減り，結果的に 1 日あたり 1 隻あたり漁獲量は上昇することが期待された．しかし，実際には低下していたことから，2009 年以降ミズダコ資源の水準そのものが低下するなどの別な要因に起因したのかもしれず，本カレンダーの導入に伴って，無駄な出漁の削減には役立つものの，現段階では 1 日あたり 1 隻あたり漁獲量の増加への効果は明らかではない．
本研究により宗谷岬沖潮流カレンダーは，宗谷海峡において樽流しに適さない日を示して，樽流し漁業者の無駄な出漁を避けるようにうながす効果があると考えられた．さらに稚内水産試験場で個別に配布した本カレンダーは曳き網漁業や遊漁などに活用されていることが，カレンダーを配布する際の聞き取り調査により確認されているため，このような表層流の予測情報は，他の沿岸漁業の燃油の節約
や計画的な操業にも貢献できるようである．さらに，このような表層流予測情報は，海上で作業を行う定置網の漁業者からも情報提供の要望が寄せられており，網揚げ作業の可否を判断するうえでも貢献しそうである．そのため，他の海域でも当海域のように潮汐情報よりも進んだ表層流予測情報を提供することで，沿岸漁業の効率化に貢献できると思われる．ただし，本計算手法では，潮汐と海流の影響が強いことが必要であり，これらの影響が弱く，風の影響が強い海域では予測精度が低下するだろう．そのため，漁業現場へ表層流予測情報を提供する際には，始めに継続的な流速の観測を行って海域の特性を明らかにして，表層流の予測精度が沿岸漁業の実用に耐えうるものかどうかを検討する必要があるだろう．
謝 辞
本研究を行うにあたり，宗谷漁業協同組合の皆様には多大なご協力を頂いた．心よりここに厚くお礼申し上げる．本
研究は，農林水産委託研究事業「新たな農林水産施策を推進する実用技術開発事業」の「宗谷海峡の空間情報統合によるミズダコ資源管理システムの開発（平成 18–21 年度）」を活用した成果である．ここに記して謝意を表す．
引用文献

土木学会海岸工学委員会（2001）「陸上設置型レーダによる沿岸海洋観測」．土木学会，東京，212 pp.
Ebuchi, N., Y. Fukamachi, K. Ohshima, K. Shirasawa, M. Ishikawa, T.
Takatsuka, T. Daibo and M. Wakatsuchi（2006）Observation of theSoya Warm Current using HF ocean radar. J. Oceanogr., 62, 47–61.
Ebuchi, N., Y. Fukamachi, K. I. Ohshima and M. Wakatsuchi（2009）Subinertial and seasonal and variations in the Soya Warm Current revealed by HF radars, coastal tide gauges, and bottom-mounted
ADCP. J. Oceanogr., 65, 31–43.
北海道水産業改良普及職員協議会（2006）北海道の漁業図鑑 写
真でみる沿岸漁業最前線．北海道水産業改良普及員職員協議会，札幌．CD-R.
Ihaka, R. and R. Gentleman（1996）R: a language for data analysis andgraphics. J. Comp. Graph. Stat., 5, 299–314.
Odamaki, M.（1994）Tides and tidal current along the Okhotsk coast of
Hokkaido. J. Oceanogr., 50, 265–279.
佐野 稔・坂東忠男・三原行雄（2011）宗谷海峡におけるミズダ
コの成熟状態の季節変化．日本水産学会誌，77, 616–624.
佐野 稔・坂東忠男・江淵直人・高柳志朗（2012）宗谷海峡のミ
ズダコ樽流し漁業における漁具の漂流速度と漁獲量の関係．
水産海洋研究，76, 123–130.
柳 哲雄（1993）「海洋観測データの処理法」．恒星社厚生閣，東京，
113 pp.

— 148 —