伊万里湾における数値シミュレーションモデルを活用した
2017 年カレニア赤潮発生期の窒素・リンの負荷量の推定
Estimation of Loads of Nitrogen and Phosphorus in Imari Bay during the 2017 Karenia mikimotoi red tide Utilizing
Numerical Simulation Model
三浦 浩*・伊藤 靖**・内田 智***・岩永俊介****Hiroshi MIURA, Yasushi ITO, Satoshi UCHIDA, and Shunsuke IWANAGA*(一財)漁港漁場漁村総合研究所

第 2 調査研究部

次長**(一財)漁港漁場漁村総合研究所

第 2 調査研究部部長

***
****

長崎県水産部参事監

長崎県水産部漁港漁場課

課長補佐

In July, 2017, a lot of cultured fish were killed by a large-scale red tide caused by Karenia mikimotoi inImari Bay which extends over Nagasaki Prefecture and Saga Prefecture. The damage cost Nagasaki
Prefecture over 600 million yen. In consequence of this, the local fishery operators request administrative
bodies to take drastic measures such as bottom sediment improvement for mitigation of damages from
red tide, in addition to conventional measures including monitoring of red tide and taking preventive
measures at the time of outbreaks thereof. To obtain basic data that are helpful for consideration of theabove matter, this study shows organized previous knowledge and the results of field survey. Then, weconstructed a numerical simulation model of the inside of Imari Bay to estimate the loads of nutrient salts.
Key Words : Imari Bay, red tide, Karenia mikimotoi, numerical simulation, nutrient salt,loads

溶出を抑えることで赤潮の発生及び拡大の抑制を期待するものである．一般的に，赤潮の発生は栄養塩，日射，海水温，流れ等の要因が複雑に絡み合うことから，底質改善
長崎県と佐賀県にまたがる伊万里湾では，1978 年から による抜本的な赤潮の予防は難しい課題であると考えら2017 年の 40 年間に 125 件(3 件/年)の赤潮が発生し，漁業 れている．
そこで，本研究では，上記の検討に資する基礎資料を得
被害は計 10 件発生しており，全て有害プランクトンを優占種とする有害赤潮によるものである．内訳はカレニア赤 るため，既往知見の整理と現地調査を行い，伊万里湾内の潮が 8 件，コクロディニウム赤潮が 2 件となっている 1)． 数値シミュレーションモデルを構築し，2017 年のカレニ
また，2017 年 7 月に発生した大規模なカレニア赤潮によ ア赤潮発生期における栄養塩の負荷量を算定した．
り，大量の養殖魚が斃死し，長崎県では約 6 億円を超える甚大な被害が生じた．
2．材料と方法
これまで伊万里湾では，大きな漁業被害に遭遇する度に，赤潮被害軽減に地域が一体となり，取組をおこなっている．
2.1 調査場所
長崎県では，発生初期から，赤潮モニタリング調査を継続すると共に餌止め等，被害防止に向けた指導を行うととも
調査場所は伊万里湾(東西約 20km×南北約 30km の範囲)に，赤潮発生時の被害軽減策として，監視体制の強化を図 である(図-1)．伊万里湾は，佐賀県と長崎県の県境に存在った．さらに，カレニア等有害赤潮が引き起こす漁業被害 する湾で，断層によって地層が落ち込み，海水が入り込んの防除を目的に，2017 年 11 月に「伊万里湾赤潮対策ガイ でできたものである．湾口幅は 4.2km，面積は 120km2，最ドライン」を策定するための検討会を設置し，養殖業者， 大水深 56ｍ，平均水深 23ｍである．大小多数の島々が点
漁業協同組合，市役所，長崎県水産部，国立研究開発法人 在し，外海からの風や波をさえぎるため，穏やかな内海と水産研究教育機構で協議を重ね，2018 年 4 月に「伊万里湾 なっている．湾内水は，東から日比水道，青島水道，津島
赤潮対策ガイドライン」を策定した 1)．
水道により外海に通じている．伊万里湾への流入河川としこうした取組の一方で，抜本的な対策として底質改善の ては，2 級河川の伊万里川，伊万里川，志佐川をはじめ，取組が地元漁業者等から要望されている．これは底質改善 大小約 20 の河川が流入する．湾奥部には干潟が形成され，により赤潮の発生要因の一つである栄養塩の底質からの カブトガニをはじめ，多種多様な生物が生息している．こ1．はじめに― 45 ―

うした海域特性を利用して，湾口から湾央部では，ぶり類，2.4 底質調査
まだい，ふぐ類，くろまぐろ等魚類養殖が，湾央から湾奥部ではまがき，あこやがい等の貝類養殖が盛んに行われて2018 年 8 月 6 日に伊万里湾内 5 地点(赤潮発生域 3 地
いる．本研究では，伊万里湾の海域特性をふまえ，湾内を 点，養殖場周辺 2 地点)において，底質調査を実施した(図湾口東，湾口西，湾央東，湾央西，湾奥の 5 区分した．
-2)．採泥にはスミス・マッキンタイヤ型採泥器(1/20 ㎡)を用いた．同一地点で 3 回採泥し，表層を採取した後，バットに移し入れ，小石・貝殻・動植物等の異物を取り除き，十分に均一化したものを分析試料とした．採取した試料は，粒度組成，化学的酸素要求量，全窒素(T-N)，全りん(T-P),
全硫化物(T-S)について分析を行った．

2.5 溶出速度試験
底泥からの栄養塩の溶出速度を算定するために，室内試験を実施した．試料の採取は，底質調査を実施した地点のうち，①③⑤の 3 地点において，20cmφの柱状採泥をダイバーにより行い，実験に用いる試料とした．

図-1 対象海域

2.2 資料収集整理
当該水域における調査・研究成果等を収集・整理するとともに，赤潮対策として不足する情報をとりまとめた．収集・整理の対象は，地形，水深，気象・海象，水質(流入も含む)，底質，赤潮発生状況，養殖業の実態並びにカレニ
ア ミキモトイに関する赤潮対策とその効果等とした．
図-2 底質調査点

2.3 水質調査2019 年5 月～7 月に伊万里湾の沖合域の1 点(N33°27′30″，E129°38′45″)において，採水器によるサンプリングと水温・塩分観測を実施した(図-1)．なお，サンプリ
ング及び水質観測は長崎県総合水産試験場の御協力により実施した．上記の採水試料を入手し，有機体窒素(TOC)，全窒素(T-N)，全リン(T-P)，アンモニア態窒素(NH4-N)，亜硝酸態窒素(NO2-N)，硝酸態窒素(NO3-N)，リン酸態リン(PO4-P)，クロロフィル a の各項目の分析を行った．分析検体数は，1 地点×3 層(海面下 1m，1／2 水深，海底上 1m)×3 回(5 月，6 月，7 月)の全 9 検体とした．分析結果は，水質モデルの外海境界データとして利用し，水質の再現性向上を図った．

図-3 溶出試験装置― 46 ―

採取した試料は，速やかに試験室に搬入し，試験を行った．試験装置は，図-3 に示したものに準拠し，直上水の水温，DO は底質調査時に観測した値を参考に設定した．試験は温度を一定に保つことが可能な恒温室内で行い，開始時及び一定時間後(0,1,2,4,8 時間後の合計 5 回)に分析試料を採水し，アンモニア態窒素(NH4-N)，亜硝酸態窒素(NO2N)，硝酸態窒素(NO3-N)，リン酸態リン(PO4-P)について分析を行った．また，20 検体のうち 5 検体は底泥を伴わない直上水のみ(コントロール)の分析検体とした．各時間の分
析結果より濃度変化の傾きを求め，溶出速度を算定した．

図-4 流動モデルが解析する主な現象と出力項目

2.6 養殖負荷量の設定
伊万里湾内で養殖されているふぐ類，まぐろ類，まだい，ぶりの主要 4 魚種について，養殖施設からの負荷を算定するための基本情報として，養殖漁業関係者への聞き取りを行った．聞き取り内容は，(1)養殖場の位置，(2)養殖場の魚種・養殖量，(3)幼魚導入時期・水揚げ時期，(4)餌の種類・給餌量・給餌頻度とした．
負荷量の算定にあたっては，1 尾，1 年あたりの養殖場における炭素・窒素・りんの負荷量を算定し，総養殖尾数(推定)を乗じることで，伊万里湾全体における養殖魚種別の負荷量とした．また，養殖場からの負荷量の設定位置は，新松浦市漁協への聞き取り結果に基づいて設定した．
[給餌量*]×(残餌率+糞尿率**)×(1−餌料の水分含有率***) ×(餌量の𝐶∙𝑁∙𝑃含有率***)[総負荷量]=[1 尾あたりの負荷量]×[総養殖尾数]*魚種毎の 1 年間の成長量，増肉係数，養殖期間を設定**魚種の餌料種類に応じた残餌率および糞尿率を設定***投餌する餌料種類に応じたの水分と C,N,P 比を設定

2.7 数値シミュレーション
伊万里湾内の負荷量を算定するため，2.2 資料収集整理～2.6 養殖負荷量の設定の情報 2)-16)を用い，数値シミュレーションモデルを構築して以下の検討を行った．検討期間
は大規模赤潮が発生した期間を含む 2017 年 5 月 1 日～
2017 年 10 月 1 日とした．適用する数値シミュレーション
モデルは，流動モデルと水質モデルから構成される．この2 つのモデルは，それぞれに独立したモデル(非結合モデル)であり，流動モデルで計算される移流・拡散量，水温・塩分を入力条件として水質計算を行うものである．

図-5 水平格子分割及び鉛直層分割(1) 流動計算
流動モデルが解析する主な現象と出力項目を図-4 に示す．用いる流動モデルは，一般鉛直座標系を用いたσ座標モデルとした．一般鉛直座標モデルでは図-5 に示すように，鉛直分割率を変化させ，水平方向に層数を変化させる― 47 ―

図-6 水質モデルの概要

表-1 「流動モデル」の入力条件に係る収集データ

表-2 「水質モデル」の入力条件に係る収集データなど，自由度に富んだ鉛直格子の設定が可能である．水平格子間隔は 100m～500m の可変格子，鉛直格子は全 21 層の多重σ座標系とした．層厚は，伊万里湾の密度鉛直分布を表現できるよう 15m 以浅の上層を細かく(0.5～1m)設定した．モデルの基礎方程式は，運動方程式，水温・塩分の保存式，乱流統計量の保存式である．流動計算で使用した各種パラメータを表-1 に示す．

― 48 ―

(2) 水質計算
プランクトンまでを考慮した低次生態系モデルを用いて水質計算を行った．水質モデルは炭素(Ｃ)，窒素連続の式，静水圧近似およびブシネスク近似が施され(Ｎ)，リン(Ｐ)および酸素(Ｏ)を指標元素とした物質循環モデルであり，低次生態系における物質循環を考慮した水質変化を予測するものである(図-6)．ここで流動モデルによって計算された物質の交換量および水温・塩分を境界条件として水質計算を行った．
植物プランクトンについては，カレニアミキモトイを想定した渦鞭毛藻類とした．底泥による酸素消費，栄養塩の溶出については，室内実験結果，既往の知見を参考に設定した．観測・分析結果を用いて，水質計算の再現性を検証した．
図-8 底質分析結果(3) 発生源別の窒素・りん負荷量の試算
カレニア赤潮の初期発生が確認された水深5ｍ以深の水深層に対し，①底泥からの溶出による負荷量，②給餌養殖による負荷，③他のエリアからの水平方向からの負荷，④上層からの負荷量について，エリア毎に集計を行った(図-7)．
集計期間は，カレニア赤潮の初期発生期にあたる2017年 6 月～7 月とした．

3.3 溶出速度試験
本試験結果を図-9 に，溶出量の他海域との比較を図-10に示す．無機態窒素は 8～38 mg/m2/day，無機態リンは-1.8～3.0 mg/m2/day の範囲であり，他海域と比較して中程度の値であった．試験結果の特徴は以下に示すとおりである．
・NH4-N 以外の項目は，時間経過による有意な濃度上昇はみられない・リンは好気的環境では溶出は少なく，貧酸素状態での溶出量が増加するものと考えられた

図-7 負荷量のエリア別整理

3．結果および考察
3.1 水質調査結果
水質調査結果で得られた TOC，全窒素，全りんの実測値に基づいて，外海境界の水質濃度の設定を行った．
図-9 溶出速度試験結果

3.2 底質調査結果
既往の調査結果と同程度の底質濃度であった．①と②を比較すると，養殖施設南西側の②の方が硫化物・有機物が多く，採取場所でばらつきがあることが示唆された．佐賀県側の④は湾奥の⑤と同様に硫化物・有機物の濃度が高い傾向を示した(図-8)．

― 49 ―

表-3 伊万里湾全体での養殖に伴う負荷量

ｸﾛﾏｸﾞﾛ
ﾄﾗﾌｸﾞ
ﾏﾀﾞｲ
ﾌﾞﾘ
計

炭素負荷量
(kgＣ/day)
737.6
136.0
126.9
534.4
1534.9

窒素負荷量(kgＮ/day)
81.4
22.7
21.1
89.1
214.3

リン負荷量
(kgＰ/day)
18.2
5.1
4.7
19.9
47.9

3.5 数値シミュレーション(1) 流動計算とその再現性検証
流動モデルを構築し，それによる計算を実施し，潮位，調和定数，潮流楕円，流速の平面分布，表層水温の時系列変化，水温・塩分の鉛直プロファイルを指標として計算値と観測値を比較した．これらの結果，計算値は観測値を概ね再現していることが示された．

(2) 水質計算と再現性検証
図-10 伊万里湾の溶出量と他の湾の比較

3.4 養殖場からの負荷量算定結果
伊万里湾内の養殖に伴う負荷量の算定結果，1 尾あたりの平均負荷量は，クロマグロが最も高く，次いで，ブリ，マダイ，トラフグの順であった(表-1)．また，総養殖尾数を推定したところ，トラフグが最も多く，次いで，ブリ，マダイ，クロマグロの順であった(表-2)．これらから総負荷量を推定した結果，海域全体では，炭素負荷量が1534.9kgC/day，窒素負荷量が 214.3kgN/day，りん負荷量が 47.9 kgP/day と推定された(表-3)．内訳をみると，ク
ロマグロとブリが相対的に高く，トラフグやマダイの約 4倍の窒素，リン負荷量であることが明らかになった．
表-1 給餌に伴う 1 尾あたりの平均負荷量

ｸﾛﾏｸﾞﾛ
ﾄﾗﾌｸﾞ
ﾏﾀﾞｲ
ﾌﾞﾘ

炭素負荷量
(mgC/尾 day)37479.5
155.2
269.6776.7

窒素負荷量(mgN/尾 day)
4135.225.9
44.9129.5

リン負荷量(mgP/尾 day)
924.55.810.028.9

表-2 伊万里湾全体での総養殖尾数の推定生産量(t)
出荷時重量(kg)
生産尾数
養殖期間(年)
歩留まり
総養殖尾数

ｸﾛﾏｸﾞﾛ
320.5
60.05,342
3.5
0.95
19,680ﾄﾗﾌｸﾞ
571.0
1.2496,522
1.5
0.85876,215

ﾏﾀﾞｲ
423.5
1.5
282,333
1.7
1.00
470,556ﾌﾞﾘ
1668.5
5.0
333,700
2.0
0.97
688,041

水質モデルにより，窒素・りん濃度及びクロロフィル aを指標として計算値と観測値を比較した．これらの結果，計算値は観測値を概ね再現した．一方で，クロロフィル a
は赤潮発生初期の再現は可能であったが，ブルーム期の再現ができない時期・地点も一部でみられた．

(3) 発生源別の窒素・りん負荷量の試算
計算結果から得られた窒素の負荷量として，発生源別の窒素負荷量を図-11 に示す．
伊万里湾全域・全層では，「河川から」の負荷が最も高く 43％を占め，「底泥からの溶出」とあわせて約 70％を占めた．一方で，降雨や養殖による負荷量はいずれも 3％以下と少なかった．カレニア赤潮の増殖層と考えられている5ｍ以深では，底泥の割合が 63％と高かった．
カレニア赤潮の初期発生水域と考えられる湾奥や湾央東海域の 5m 以深では「底泥からの溶出」がそれぞれ 26％と 34％であった．
計算結果から得られたリンの負荷量として，発生源別のリン負荷量を図-12 に示す．伊万里湾全域・全層では，「河川から」の負荷が最も高く59％を占め，「底泥からの溶出」とあわせて，70％以上を占めた．一方で，降雨や養殖による負荷量はいずれも 2％以下と少なかった．カレニア赤潮の増殖層と考えられている 5ｍ以深では，底泥の割合が57％と高かった．カレニア赤潮の初期発生水域と考えられる湾奥や湾央東海域の 5m 以深では「底泥からの溶出」がそれぞれ 39％と 49％であった．

― 50 ―

5m 以深（ｶﾚﾆｱ増殖層）

全層

全層

エリア

DIN総量：574トン

底泥からの溶出

DIN総量：187トン

養殖漁場から

全域

他の海域から

河川
42.8%

底泥
25.1%

上層から

上層
19.6%

養殖
0.8%

5m 以深

他海域
15.7%

養殖漁場から

底泥
62.9%

他海域
28.0%

湾口西

湾口東

DIN総量：287トン
降雨
2.0%

養殖
1.1%

河川 底泥
12.3% 13.4%

降雨
6.6%
他海域
17.7%

DIN総量：124トン

養殖
2.3%

上層
2.4%

底泥
28.8%

養殖
2.0%

他海域
35.5%

底泥
55.5%
降雨8.5%

DIN総量：65トン

降雨
1.7%

底泥
15.6%

河川82.7%
養殖
0.5%

DIN総量：13トン

DIN総量：70トン

底泥34.2%

底泥
25.7%

上層
45.1%
養殖
0.2%

上層
65.6%

底泥
74.2%

底泥
97.7%

DIN総量：211トン

底泥
49.3%

他海域
14.7%

他海域
22.7%

養殖
0.7%

増殖層

河川
26.9%

養殖
1.0%

養殖
1.4%

DIN総量：11トン

湾奥

DIN総量：19トン

河川
19.3%

底泥
45.6%

降雨
5.4%

上層
33.7%

降雨から

湾央東
DIN総量：91トン

他海域
71.2%

ｶﾚﾆｱ

他の海域から

増殖層

河川から

湾央西

DIN総量：29トン

河川
47.5%

5m 以深

ｶﾚﾆｱ

養殖
1.7%

降雨
3.2%

全層

底泥からの溶出他海域
29.2%

図-11 DIN を例にしたエリア別・要因別の濃度増分(2017 年 6～7 月平均，単位はトン/月)
5m 以深（ｶﾚﾆｱ増殖層）全層
エリア

DIP総量：61トン

全層

底泥からの溶出養殖漁場から

DIP総量：22トン

他の海域から

全域

上層から

底泥
23.8%

養殖
1.6%

5m 以深

上層
34.8%

底泥からの溶出養殖漁場から

河川
58.8%

底泥
56.5%

他海域
15.6%

湾口西

湾口東

河川
27.1%

底泥
17.0%

底泥
42.2%

河川
53.8%

河川
24.0%
降雨0.3%

他海域
52.7%

降雨0.4%

5m 以深
ｶﾚﾆｱ

DIP総量：9トン

上層
41.2%

湾央東

DIP総量：11トン

降雨0.4%

養殖
1.8%

DIP総量：8トン

DIP総量：1トン
養殖
6.0%

上層
7.8%

底泥
38.8%

他海域
39.0%
養殖
5.9%

DIP総量：29トン

底泥
94.0%

養殖
1.3%

底泥
52.0%

降雨0.1%河川
33.4%

底泥
39.4%

他海域
34.4%

養殖
3.7%

増殖層
他海域
14.1%

湾奥
DIP総量：2トン

底泥
14.1%

養殖
3.1%

降雨
0.2%

降雨から
河川から

湾央西
DIP総量：3トン

DIP総量：22トン

他の海域から増殖層

他海域
5.5% 養殖
3.2%

降雨
0.2%

全層

ｶﾚﾆｱ

他海域
1.2%

底泥
63.9%

河川
85.8%

養殖
1.0%

DIP総量：2トン

底泥49.0%

上層
50.7%
養殖
0.3%

DIP総量：8トン

底泥
39.0%

上層41.6%他海域
19.3%

図-12 DIP を例にしたエリア別・要因別の濃度増分(2017 年 6～7 月平均，単位はトン/月)

4．おわりに
本研究では，伊万里湾における大規模赤潮が発生した2017年 6～7 月の流動特性および水質動態をふまえ，カレニア
赤潮の初期発生域における増殖を抑制することを検討するために，湾内の栄養塩濃度分布の推定を行った．その結果，海域別，水深帯別の負荷量を算定するとともに，カレニア赤潮の初期発生場所と考えられている湾奥や湾央東の海域における底泥からの溶出量について数量化を試みた．今回得られた知見をふまえ，今後も引き続き，基礎データを収集し，計算データに基づいて，被害軽減のための対策について検討を行うことが望まれる．

― 51 ―

5．謝辞

10) 鈴木誠二・多田彰秀・西田渉・坂口正人：汚濁流入負荷を考慮した伊万里湾における流動特性及び水質動態に関する研究,

本研究は，当研究所が平成 29 年～令和元年に長崎県よ土木学会論文集 B2(海岸工学)，vol.67,No.2, pp.I_831I_835,2011.
り受託した「伊万里湾漁場環境改善対策検討調査業務委託」
等の成果の一部をとりまとめたものであり，現地調査の実 11) 高木雅子・村瀬茂世・寺田和光・松原英隆: 博多湾底質からの窒素,りんの溶出に関する基礎研究(第 5 報), 福岡市保環施に際しては，新松浦漁業協同組合，長崎県総合水産試験
研報,23,1998.
場の協力の下で実施した．また，調査・検討結果等について専門分野の知見を有する学識経験者を委員とする検討 12) 西田修三・中谷祐介・島田恭佑・入江政安: 降水の水質特性と大阪湾への栄養塩負荷の影響解析, 海岸工学論文集,第 55
委員会を 3 ヵ年で計 7 回開催し，検討委員会での議論を経巻,pp.1061-1065，2008.て，成果に反映した．ここに記して深謝申しあげます．
13) 船越茂雄: 伊勢・三河湾における漁業生産による窒素,リンの回収, 愛知県水産試験場報告第 14 号, pp.1-6,2008.

参考文献

14) 山口峰生: Gymnodinium nagasakiense の赤潮発生機構と発生

1) 伊万里湾赤潮対策検討会議：伊万里湾赤潮対策ガイドライン．

予知に関する生理生態学的研究, 南西水研研

pp.1-17, 2018.

報,No.27,pp.251-394,1994.

2) 飯塚昭二: 大村湾における Gymnodinium type-’65 年型種赤15) 山田達夫・橋本俊也・末永慶寛・一見和彦・多田邦尚: 魚類

潮の発生機構, 日本プランクトン学会報第 19 巻第 1 号,1972.

(ハマチ)養殖場における残餌の沈降と拡散, 日本水産学

3) 今村正裕・松梨史郎：湾奥部における水-底泥間の窒素・リン会,75(4),pp.666-673,2009.のフラックス，
海岸工学論文集,第 44 巻，,pp.1081-1085,1997.

16) Yoshitsugu Koizumi,Takuji Uchida and Tsuneo Honjo:

4) 神山孝史・辻野睦・玉井恭一: 夏季成層期の播磨灘海底におDiurnal vertical migration of Gymnodium mikimotoi

ける栄養塩類溶出量, 南西水研研報, No.31,pp.33-43,1998.

during a red tide in Hoketsu Bay,Japan, Journal ofPlankton Research Vol.18 no.2,pp289-294,1996.

5) 北辻さほ・紫加田知幸・鬼塚剛・青木一弘: 八代海における赤潮発生メカニズムの解析, 九州海域での有害赤潮・貧酸素
水塊発生機構解明と予察･被害防止等技術開発報告書,2017.
6) 小泉喜嗣・高島景・神薗真人・江藤卓也: 西部瀬戸内海における Gymnodinium mikimotoi の増殖域の環境特性と分布拡大

関連情報1) 29 漁港環第 4 号 伊万里湾漁場環境改善調査業務委託報告書

機構, 海の研究 Vol.3,No.2,1994.

平成 30 年 3 月 一般財団法人 漁港漁場漁村総合研究所

7) 郡山益美・荒木啓輔・伊藤祐二・片野俊也・石谷哲寛：有明2) 30 漁港環第 1 号 伊万里湾漁場環境改善対策検討調査業務委

海奥部底泥における栄養塩フラックスの季節変化，土木学会

託報告書 平成31 年 3 月 一般財団法人 漁港漁場漁村総合

論文集 B2(海岸工学)Vol.68, No.2, pp.Ⅰ_1056-Ⅰ_1060，

研究所

2012.

3) 01 漁港環第 1 号 伊万里湾漁場環境改善対策検討調査業務委

8) 鈴木誠二・西田渉・平智樹・野口正人：伊万里湾の流動構造に及ぼす風の影響に関する研究．海岸工学論文集，第 55 巻，pp.396-400，2008.
9) 鈴木誠二・多田彰秀・平智樹・坂口誠人・松本健太郎・西田
渉: 現地観測及び流動水質モデルを用いた伊万里湾における流動・水質特性の把握,土木学会論文集 B2(海岸工学)，vol.66,No.1, 1021- 1025,2010.

― 52 ―

託報告書 令和 2 年 3 月 一般財団法人 漁港漁場漁村総合研究所