福 岡 水 海 技セ 研 報

第 31号

2021年3月

Bull.Fukuoka.Fisheries.Mar.Technol.Res.Cent.No31.March 2021

たる流し漁業におけるケンサキイカの漁獲特性
長本篤（水産海洋技術センター）

たる流し漁業では，ケンサキイカの漁獲に漁具の揺れや水温が影響していることが経験的に知られている。このような漁獲特性の把握は，効率的な操業方法の検討や CPUE 標準化による資源変動の推定を行う上で重要である。本報告では同漁業の漁獲特性を把握するため，本県筑前海区に来遊する春季及び夏季成熟群を対象に，漁獲尾数に影響を与える海洋環境要因について検討した。海洋環境要因を説明変数とした一般化線形モデルによる解析の結果，漁獲尾数に大きな影響を与える最良モデルの変数は，底層水温と有義波高であった。また，漁獲尾数は有義波高約1.6m までは比例して増加した。このことから，有義波高データを用いることにより，効率的な操業方法の検討や CPUE 標準化による資源変動推定の可能性が示唆された。
キーワード：ケンサキイカ，たる流し，有義波高，CPUE 標準化作業を1日に6～8回繰り返す。漁具は回次あたり約1時間

ケンサキイカ日本海・東シナ海系群の資源水準は2000
1）

年以降低位で推移し ，福岡県筑前海区におけるケンサ波に揺られるが，その揺れがケンサキイカの漁獲に影響

2）キイカの漁獲量は1992年以降減少している 。

していることが経験的に知られている。

水産資源の評価では，資源量指標値として単位努力あ

本報告では，たる流し漁業におけるケンサキイカの漁たり漁獲量（以下，CPUE）が用いられる。一般に CPUE獲特性把握のため，漁業者が記帳した操業日誌及び水温，は資源密度に比例していると考えられており，CPUE の波高などの海況データを用い，ケンサキイカの CPUE

年傾向をみるだけで相対資源量の増減傾向を把握するこに影響を及ぼす海洋環境要因について検討した。また，

3）とが可能である 。しかし，CPUE は季節や海区，漁具

波による漁具の揺れ及び漁獲への影響についても検討しなど資源密度以外の様々な要因も含んでおり，資源変動た。

を正確に把握するためには CPUE 標準化によりこれらの影響を取り除く必要がある3）。

方法

他海域では CPUE 標準化によりケンサキイカ釣り漁業の CPUE に影響を与える要因について明らかされて本県筑前海区に来遊するケンサキイカは，春季成熟群，4）いる 。本県においても，ケンサキイカの漁獲特性を把夏季成熟群，秋季未熟群の3群の存在が知られている6）。

握し，CPUE 標準化により資源密度以外の要因を取り除本研究では春季及び夏季成熟群を対象に，漁獲尾数及びき，より精度の高い資源変動傾向を把握する必要がある。

漁獲量と海洋環境要因との関係，波高とたる流し漁具の本県筑前海区のケンサキイカは主にたる流し漁業を含挙動との関係を検討した。各々の使用データ及び解析方

むいかつり漁業で漁獲されており，1982年から夜いかつ法は次のとおりとした。

り漁業と昼いかつり漁業に加え，たる流し漁業による漁獲が始まり 5），特に主要漁協ではその漁獲割合が高くな１．漁獲尾数と海洋環境要因との関係っている。たる流し漁業では，図1のように直径30cm 程たる流し漁業における操業回次ごとのケンサキイカ漁

度の発泡スチロール製の尺玉（タル）に6本程度のいか獲尾数と海洋環境要因との関係を，操業データ及び海況針を付けたテグスを取り付け，海底付近に垂下する漁具データを用いて検討した。

を用いる。1隻あたり20個程度の漁具を約100m の間隔（１）解析に用いたデータで順次投入し，約1時間潮に流して操業する。全ての漁操業データは，協力漁業者１名が2019年4月13日～8月

具の投入が終わると最初に投入した漁具から引き上げ，9日に記帳した操業1回次及び最終回次の緯度経度，操業掛かったイカを外して再び漁具を投入する。この一連の回次ごとの開始，終了時刻と漁獲尾数を用いた。各- 29 -

長本

MuMIn パッケージの Dredge 関数を用いて総当たり法により赤池情報量規準（AIC）が最小となる最良モデルを選択した。
（式1）漁獲尾数=月+表層水温+底層水温+表層塩分+底層塩分+表層流速+底層流速+表層流向区分+底層流向区分+有義波高+offset（log（操業回次））
２．たる流し漁具の挙動と波高との関係
操業中のたる流し漁具の挙動と波高との関係を，深度データ及び波高データを用いて検討した。
深度は，協力漁業者１名が2020年4月9日～6月18日に
図1のとおり漁具の最下部から約10m の位置（以下，下部）に水圧計（JFE アドバンテック社製 DEFI2-D20）を取り付け1分間隔で観測した。取得した海中の深度データから漁具投入及び回収中など深度が急激に変化するデ
図1ータを除き，深度がほぼ一定になる時刻を操業中と判断

たる流し漁業の漁具及び観測機器設置位置した。
操業回次の操業位置は，1回次及び最終回次の緯度経度

波高データは，国土交通省港湾局によって観測され，を各操業回次の開始時刻と終了時刻の中間時刻で按分し港湾空港技術研究所で処理された全国港湾海洋波浪情報て求めた。

網のうち，玄界灘（N33°56′02″，E130°28′05″）

海況データとして，九州大学応用力学研究所が開発，の2020年4～6月の有義波高データ（速報値）から操業時運用している海況予測モデル DREAMS_Squid（ 水平300m間帯の値を抽出し操業回次ごとに平均したものを用いメッシュ，深度2～4m 間隔，1時間間隔，以下，DR_S）た。に観測データを同化した再解析値から，操業回次の位置と時刻に最も近い水温，塩分及び潮流の東西，南北成分３．漁獲量と波高との関係を抽出した。さらに，流向，流速を東西，南北成分からたる流し漁業におけるケンサキイカの漁獲量と波高と算出した。このうち解析には，タルが潮流等の影響を受の関係を，操業日誌データ及び波高データを用いて検討ける層として1m 深（以下，表層），たる流し漁業がケした。

ンサキイカを狙う底層として DR_S の最深層値から概ね操業日誌データは，漁業者46名が2007～2019年の4～84m 上層（以下，底層）のデータを使用した。

月に記帳した1日あたりの漁獲量及び緯度経度2分メッシ
ュの漁区に区分された操業場所のデータ（計4,782件）

波高データは，気象庁が公開している沿岸代表点の波を用いた。

浪推定値のうち，玄界灘（N34°15′，E130°00′）の午前９時の有義波高データを用いた。

波高データには，気象庁が沿岸波浪数値予報モデルから計算した有義波高 GPV（0.05°メッシュ）を用いた。

（２）解析方法
統計解析パッケージ R（Ver.3.5.2）の一般化線形モ

解析では，操業した漁区の中心の緯度経度に最も近い午デル（GLM 関数）を用いて解析した。今回の操業デー前9時のデータを使用した。

タには，ケンサキイカが獲れないゼロキャッチデータは存在しなかった。

結果

一般化線形モデルとして，式1のとおり漁獲尾数を応答変数に，月，水温（表層，底層），塩分（表層，底層），１．漁獲尾数と海洋環境要因との関係

流速（表層，底層），流向区分（表層，底層），有義波高

今回の解析に用いた操業データの操業場所を図2，月を説明変数としたフルモデルを作成した。流向区分は真別漁場別の操業回数を図3に示す。漁場は，4～8月の沖方位0°から90°ごとに4分割した流向区分 A ～ D のカノ島以南の漁場（以下，漁場 A）と4，5月の沖ノ島以北テゴリカル変数とした。漁獲尾数は正の整数であるため，の漁場（以下，漁場 B）に大別された。全データ計255

誤差構造は負の二項分布を仮定した7）。フルモデルから- 30 -たる流し漁業におけるケンサキイカの漁獲特性

表1

操業回次別操業開始時刻及び件数

操業回次時刻1

24:00 ～5:00 32

1

沖ノ島

5:00 ～6:00

漁場A

6:00 ～

7:00

30

57:00 ～8:00

229

58:00 ～9:00

323

7719漁場B

53操業場所

5691516411:00 ～ 12:00

56212:00 ～ 13:00

2213:00 ～ 14:00

12136

表2373737353226

PIntercept

-0.66

0.45

0.14

底層水温

0.16

0.03

<0.01

表層流速-0.76

0.42

<0.1

底層流向区分B

0.16

0.29

0.60

底層流向区分C

-0.40

0.28

0.15

底層流向区分D

-0.38

0.11

<0.01

有義波高

0.81

0.16

<0.01

件のうち，漁場 A のデータが224件と約9割を占めたこ

表3とから，今回の解析では漁場 A のデータを用いた。
操業回次別の開始時刻及び件数を表1に示す。1日当たりの操業回数は2～9回で，午前4～5時に操業を開始し，遅くても13～14時に終了していた。
操業及び海況データを用いて一般化線形モデルによる1

最良モデル説明変数の偏回帰係数及び標準誤差

偏回帰係数 標準誤差

月別漁場別操業回数

810:00 ～ 11:00

計

図3749:00 ～ 10:00

図24

説明変数及び ΔAIC

説明変数

ΔAIC

底 層 水 温38

表 層 流 速1

底層流向区分

6

有 義 波 高21

解析を行った結果，フルモデルの説明変数のうち，月，表層水温，底層水温，表層塩分，底層塩分で多重共線性

底層水温，有義波高，底層流向区分，表層流速の順で大が認められたことから，底層水温，表層流速，底層流速，きく，特に底層水温，有義波高が大きかった。なお，最表層流向区分，底層流向区分，有義波高を変数とするモ良モデルでは多重共線性，線形性に問題はなかった。

デルを作成した。総当たり法により選択された最良モデ

0.1℃ごとに区分した底層水温とその区分内の漁獲尾ルの変数は，式2のとおり底層水温，表層流速，底層流数との関係を図4に示す。底層水温は15.2～21.4℃，漁向区分，有義波高であった。最良モデルの変数は，表2

獲尾数は2.5～42尾/回次であった。漁獲尾数と底層水温のとおり Intercept，底層流向区分 B，底層流向区分 C をの相関係数は0.607となり，正の相関がみられた。（無相除き有意であった。

関検定

（式2）

p<0.01，n=34）0.1m ごとに区分した有義波高とその区分内の漁獲尾漁獲尾数=底層水温+表層流速+底層流向区分数との関係を図5に示す。有義波高は0～1.6m，漁獲尾+有義波高+offset（log（操業回次））

数は5～34尾/回次であった。漁獲尾数と有義波高の相関係数は0.673となり，正の相関がみられた。（無相関検定
p<0.05，n=13）

最良モデルから各説明変数を除いたモデルと最良モデ
ルの AIC の差（以下，ΔAIC）を表3に示す。ΔAIC は，- 31 -

長本

図4

底層水温と漁獲尾数の関係

図6

図5

下部の深度の推移

有義波高と漁獲尾数の関係

２．たる流し漁具の挙動と波高との関係
各操業回次のうち，有義波高が最大であった2020年5月21日の1回次（有義波高2.2m）と最小であった2020年6月4日の3回次（同0.2m）における下部の深度の推移を
図6に示す。5月21日の下部の深度は50.6～53.1m，6月4日の下部の深度は79.4～80.0m の範囲で推移した。
各操業回次における下部深度の最大値と最小値の差（以下，変動幅）と有義波高との関係を図7に示す。変

図7

下部の深度の変動幅と有義波高との関係

動差は0.1～2.7m，有義波高は0.2～2.2m であった。変
動幅と有義波高の相関係数は0.736となり，正の相関がみられた（無相関検定

p<0.01，n=87）。

３．漁獲量と波高との関係
0.1m ごとに区分した有義波高とその区分内の漁獲量との関係を図8に示す。漁獲量は7～31kg/日，有義波高は0～3.1m であった。漁獲量は，有義波高0m で最も少なく7kg/日，有義波高2m で最も多く31kg/日であった。
有義波高1.6m までは7～29kg/日となり，漁獲量と有義波高との相関係数は0.794となり正の相関がみられた（無相関検定

p<0.01，n=17）。有義波高が1.7m 以上では

図812～31kg/日となり，ばらつきが大きくなった（無相関検定有意差なし，n=14）。

- 32 -

有義波高と漁獲量の関係

たる流し漁業におけるケンサキイカの漁獲特性考

カイカが疑似餌を捕捉するときの巻き上げ速度が速いと察

きは触腕で，遅いときは触腕以外で疑似餌に掛かる傾向
一般化線形モデルによる解析の結果，漁獲尾数に大きがあることが報告されている12）。このことから，有義波

な影響を与える変数として底層水温と有義波高が選択さ

高1.7m 以上の条件では，いか針の動きが大きくなるこれた。とや速くなることにより，ケンサキイカの漁獲に悪影響を与えている可能性が考えられた。

底層水温と漁獲尾数の関係をみると，4～8月の漁場 A
では底層水温に比例して漁獲尾数が増加した。筑前海区

今回の解析により，たる流し漁業では有義波高が漁獲におけるケンサキイカは，6～8月にかけて沖合から沿岸尾数に大きく影響することが示唆されたことから，漁業域まで広く分布するが，次第に沿岸域への来遊が少なく種類によっては，波浪による漁具の挙動や船体の揺れがなり9月以降はみられなくなる 。また，ケンサキイカがCPUE 標準化の際に密度効果以外の重要な要因になると漁獲される水温範囲は18～24℃であり，20℃付近が盛漁考えられる。今後は，複数の漁業者の操業データを利用期の水温である 。漁場 A は，水深40～80m の沿岸域でしてデータ数を増やし，解析の精度を向上させる必要が4～8月にかけて底層水温が上昇する海域であり，底層水ある。

8）9）温の上昇に伴いケンサキイカが漁場に来遊したと考えらたる流し漁業では対象となるケンサキイカ以外にブリれる。本研究では，漁場 A で水温が上昇する期間（4～類，ハタ類などの魚類も混獲される。魚類はいか針に掛8月）のみ解析を行ったが，CPUE 標準化により年変動かったケンサキイカを捕食することもあるが，直接いかをみる場合は，底層水温の長期的な変動を考慮する必要針に掛かることもあり，魚種により掛かりやすいいか針がある。の動きがあると考えられる。いか針の動きは波高以外に有義波高と漁獲尾数の関係をみると，有義波高が高く内部波など潮流の影響も受けると考えられることから，なると漁獲尾数も多くなった。たる流し漁業では，漁具よりケンサキイカの選択性を高めるためには，潮流等のを約1時間潮に流して操業するため，漁具は操業中に波影響も考慮した針の動きと漁獲尾数との関係の把握が必浪の影響を受けている。また，下部の深度の変動幅と有要である。この関係解明により，波高が低いときでもい義波高の関係をみると，有義波高が高くなると下部の深か針が動く効率的な漁具の開発も期待される。

度の変動幅も大きくなった。テグスといか針を結ぶ糸の謝長さは約50cm であることから，漁具の動きが有義波高

辞に比例して大きくなることにより，いか針が大きく動くと考えられる。イカ釣り機に最も重要なことは釣糸を引

本論文の執筆にあたり操業データを提供していただい

く速度に変化をつけて，釣り針を餌の動きのようにみせた宗像漁業協同組合の桑村勝士代表理事組合長及び長年るシャクリ機能である10）。このことから，有義波高に比にわたり操業日誌の記帳に協力して頂いた漁業者の皆様

例していか針の動きが大きくなり，シャクリのような機に感謝の意を表する。また，本報告の解析にあたり，国能が働くことにより漁獲尾数が増加したものと考えられ

立研究開発法人水産研究・教育機構水産資源研究所大下た。

誠二副部長には有益な御助言を頂き感謝の意を表する。

図8で示した有義波高と漁獲量との関係をみると，有義波高と漁獲尾数との関係（図5）と同様に有義波高が~文

献1.6m の間は有義波高が高くなると漁獲量も多くなった。
一方，有義波高が1.7m 以上になると漁獲量のばらつき

１）依田真里，高橋素光．平成29（2017）年度ケンが大きくなり，漁獲量は有義波高0~1.6m のときより少サキイカ日本海・東シナ海系群の資源評価．平成29なくなることがあった。操業日誌データは，複数の漁業年度我が国周辺水域の漁業資源評価（魚種別系群別者が記帳したものであるため，漁船の規模により波高が資源評価・TAC 種以外）第3分冊，水産庁増殖推

高くなると操業を中止して帰港する漁業者が含まれ，１

進部，東京． 2018；2110-2125．

日あたりの漁獲量が少なくなることが考えられる。アカ２）ケンサキイカ共同研究報告書．ケンサキイカ研究会，

イカ釣りにおける釣り落としと船体動揺との関係では，

2019;55-76．

凪の時には船体の動揺が相対的に大きい位置の脱落率が３）庄野宏．CPUE 標準化に用いられる統計学的アプ相対的に低くなるが，時化の時には船体の動揺が大きい

ローチに関する総説．水産海洋研究 2004；68（2）船首尾で脱落率が極端に高まる場合がある11）。また，ア：106-120．

- 33 -

長本；83-91．

４）金元保之，柴田泰宙．島根県沿岸におけるケンサキ
イカの CPUE に影響を与える要因の検討と将来予

９）古田－漁場と分布－．昭和49年度福岡県福岡水産試験場

測．日本水産学会誌 2020；86（5）：371-385．
５）秋元聡．筑前海におけるケンサキイカ資源と利用

研究業務報告 1976；50-56．

実態．福岡県水産海洋技術センター研究報告 1992

10）浜出雄三．イカを騙す秘訣－自動イカ釣り機の開発．「イカの春秋」（奥谷喬司編著）成山堂書店，東；18：7-11．
６）山田久典．筑前海におけるケンサキイカについて

英明，小川嘉彦，森脇晋平，岡島義和．日本

京．1995；207-215．

海西部沿岸域におけるケンサキイカ・ブドウイカの

11）山下秀幸，黒坂浩平，越智洋介，小河道夫．

生物学的特性．日本海西部沿岸に生息する“シロイ

アカイカ釣りにおける釣り落としと船体動揺との関カ“（ケンサキイカ・ブドウイカ）に関する共同研

係．日本水産学会誌 2008；74（4）：697-699．

究報告書 1986；（1）：29-50．
７）下野嘉子．特集

統計解析（再？）入門

R を用12）黒坂浩平，山下秀幸，越智洋介，小河赤松友成，稲田博史，渡部俊広．自動イカ釣道夫，いた一般化線形モデル（回帰係数編）：カウントデ機の釣具ライン巻き上げ速度がアカイカの擬餌針捕ータを例に．雑草研究 2010；55（4）：287-294．

捉行動に及ぼす影響．日本水産学会誌 2009；75（1）８）古田久典．筑前海におけるケンサキイカについて．

昭和48年度福岡県福岡水産試験場研究業務報告 1975

- 34 -：83-85．