平成24年度

北海道沿岸における漁場環境の変動機構解明の試み―鵡川出水前後の水質変化に関する観測・解析手法―(独)土木研究所 寒地土木研究所 水産土木チーム○須藤 賢哉大橋 正臣
山本 潤

河川流入がある沿岸漁場は漁業資源が豊富であるが、台風等による河川出水があった場合には水産生物の生息環境や餌料環境に大きな影響を及ぼすことがある。出水前後の状況を正確に把握するため、適切な観測・解析手法を用いた沿岸漁場環境の変動機構の解明が求められる。
本稿は鵡川海域において台風1112号による河川出水を対象に現地観測と数値解析を組み合わせ、漁場環境の変動機構の解明を試みた。特に生態系モデルによる現況再現と予測計算を行い、出水による影響について、基礎生産量を指標として考察を行った。
キーワード：沿岸漁場環境、河川出水、基礎生産量

1. はじめにに、数値解析を導入し、現地観測結果を補う必要がある。前報で須藤ら 3) は閉鎖性海域を対象に貧酸素化や富栄養化の問題に対する適用事例を示しているが、本報では開放性海域における河川出水前後の水域環境変化についての検討手法を示すものである。対象としたフィールドは鵡川海域として、台風 1112 号による鵡川の河川出水について現地観測と数値解析により沿岸漁場環境の変動機構の解明を試みたものである。

沿岸漁業の生産量は 1985 年前後に 220 万トンを上回っていたが、漁場環境の悪化や資源の低迷等により減少傾向となり、現在は約 130 万トン(2010 年)となっている
1)。現在でも海面漁業生産量の約 3 割を占め、我が国の水産物供給に大きく貢献している。沿岸漁場は水産生物の生息場や育成場として重要であると共に、漁業・海運等の利用、港湾・漁港建設等の開発、陸域からの負荷の影響を受けやすい水域環境でもある。中でも河川が流入する水域は、陸域からの栄養塩供給(図-1)により漁業資源が豊富な沿岸漁場といわれている 2)。今後も漁業生産
2. 水域環境把握のための現地観測と数値解析の維持には、河口周辺の水域環境に配慮する必要がある。
水域環境を把握するためには、現地観測と数値解
この河口周辺の水域環境を理解するためには、土砂や析を組み合わせることが有効である。特に数値解
栄養塩に着目して、物理場・水質・生物等の基礎データ析を導入するメリットは、現地観測だけでは網羅でを取得し解析することが基本となる。しかし、現地観測きない範囲・時刻の水域環境の再現・試算ができることは費用や観測期間の制約から実施規模には限界がである。よって、現地観測で数値解析の初期条件や境界あるため、適切な現地観測計画を立案するととも条件の取得を目的とした観測計画の立案が必要である。
現地観測と数値解析の基本的項目を以下に示す。
都市からの雑排水
森林からの土砂(1)項目１：現地観測による水域環境の把握①連続観測・移動観測：流速・波高の定点観測および流速の移動観測により波浪・流況を把握する。
農業由来の肥料、糞尿等②多点観測：垂下式観測（CTD）により水温・塩分・濁度・クロロフィルａの鉛直分布を把握する。採水・水
土砂、栄養塩類の供給
質分析によりSS･栄養塩・クロロフィルａ等の水質を把握する。採水の深度は躍層位置に留意して設定する。
漁場環境の変動
河川出水（土砂、栄養塩類）の影響解明
沿岸漁場③その他：河川流量・SS等のL-Q式・気象データ（気温・風速・日射量等）を把握する。
図-1 河川が流入する開放性海域の水域環境変化
Kenya Sudo , Masami Ohashi , Jun Yamamoto3. 鵡川海域における河川出水の影響

苫小牧港東港区

北 海 道
鵡川橋
鵡川海域
鵡川鵡川漁港

Ｃ地点

沙流川
Ａ地点

22.1km12.6km(2) 項目２：数値解析による水域環境の再現・予測①海洋モデル4)：現地観測により得られた水温・塩分等に基づいて3次元的・時系列的な解析により物理特性を再現する。このモデルは、非圧縮粘性流体・ブシネ
スク近似・静水圧近似を基礎としたモデルであり、マルチレベルモデル・POM等が一般に使われている。
②低次生態系モデル5)：このモデルは、生態系を構成する生物(動植物プランクトン等)と非生物(栄養塩等)を構成要素とし、主に食物連鎖に基づいてこれらの相互作用を計算するものである。これにより栄養塩、クロロフィルａ等の3次元分布とその時間変化を求める。

数値解析範囲
Ｂ地点

鵡川海域は一級河川である鵡川・沙流川が流入する海域で、ウバガイ（ホッキガイ）やシシャモの漁獲量が北図-2 現地観測位置と数値解析範囲
海道で上位に位置し豊かな漁場が形成されている。一方、台風等による出水時には多量の浮遊土砂を含む濁水が周
表-1 現地観測内容
辺海域に流れ込み、台風0115号に伴う鵡川・沙流川の出調査項目
方法
観測内容
水では、当時行われていたホタテ養殖事業が被害を受け
6)波浪・波高計，連続観測(8/5～9/15)、それ以後は養殖事業が中止された。この原因はホタ流況観測
流速計ADCP
流速移動観測(7/29)
テガイが濁水を吸い込み擬糞として排出することで体力水質(11 測点)
CTD垂下
水温,塩分,濁度,chl-aを消耗し、死に至ったと考えられている7)。これはホタ採水,
SS,COD,NH4-N, NO2-N,NO3-N,各測点2～4層 PO4-P,SiO2-Si,T-Fe,chl-a
テガイと摂餌機構が同じ二枚貝（当海域の代表種である細胞数植物プランクトン 採水(水質と同)ウバガイ）においても同様の懸念がある。
透明度板
透明度
透明度
この他にも出水の影響は、濁水・淡水・栄養塩の移流・拡散によって基礎生産や餌料環境に影響を及ぼし、 ※水質･植物プランクトン･透明度は平水時 8/3 出水時 9/10 観測
メリットとデメリットの両面があると考えらる。当海域
観測点は場を代表する１地点とした。数値解析の境界条の既往研究としては、山下ら8)が濁水の挙動を検討した
件として調査船にADCPを装着し、図-2に示す数値解析の例はあるが、基礎生産やそれに基づく水産生物の餌料環地形図範囲（12.6×22.1km）を操行しながら流速を観測境に着目した研究は少ない。よって、出水による漁場環した。
境の変動機構を解明するためには、基礎生産量の変化について検討する必要がある。

5. 現地観測結果
4. 現地観測方法(1)河川出水に伴う濁水の挙動
図-3 には 2011 年 9 月 1 日～15 日の鵡川の河川流量
現地観測計画は以下の通りとした。観測時期は台風の（鵡川橋）とともに、沖合で観測された有義波高および発生時期である夏季から秋季とし、平水時（8月）、台流速の時系列変化を示す。降水は 9 月 2 日頃から始まり風通過後（出水時）を対象とした。実際には2011年9月
河川流量は 3 日頃から急増し、4 日に時間最大流量の台風12号に伴う鵡川の出水を捉えることが出来た。
642m3/s を記録した。台風の接近に伴い 3 日から有義波測点の配置については、既往知見として当海域の漂砂高・流速が増加し、8 日頃には比較的静穏となっている。
機構として沿岸漂砂が西向き（襟裳岬から苫小牧に向か
表層流速は東西成分（東を正）と南北成分（北を正）をう流れ）であり、なおかつ夏季出水時の衛星写真から濁示している。出水前後での流速は西と北の成分が卓越し水が西寄りに広がっていることを考慮し、鵡川河口よりている。
西側を観測箇所とした。また，ウバガイの生息水深と過9)
図-4 には 2011 年 7 月 29 日の流速移動観測による沿岸去の近隣海域における同様の調査事例 より水深15m以浅流速ベクトルを示す。流速ベクトル値は水深2～5m の平に測点を配置した（以上図-2）。観測項目（表-1）は波均値で示している。Ａ地点～Ｂ地点(岸沖方向)の流速は、浪・流況・水質（栄養塩など）・植物プランクトン・透海岸線に沿った北西方向のベクトルが卓越しており、約明度とした。波浪・流況観測では、水深30m付近までは等深線が汀線にほぼ平行で沿岸方向の変化は小さいため、 0.2m/s である。これに基づき数値解析では境界条件のKenya Sudo , Masami Ohashi , Jun Yamamoto

河川流量(鵡川橋)
流速(東西成分)表層
流速(南北成分)表層有義波高80060020-2001.5-400

沿岸流の向き：北西

沙流川

32.5

400200

鵡川

9/3

9/4

9/5

9/6

9/7

St.1-4 →●●St.2-3St.2-1 St.2-2
●

0.5

9/2

St.1-3

St.1-1 St.1-2●

1-600-800

-1000
-1200
(南、西) 9/1

3.5

有義波高(m)

河川流量(m3 /s）、流速(mm/s)

(北、東)
1200
1000

●

●

2.4km

09/8 9/9 9/10 9/11 9/12 9/13 9/14 9/15
日付

図-3 鵡川の河川流量と沖合の波高・流速N

2011年9月14日空撮
St.3-1

St.3-2 ・・・●●

写真-1 濁水の流出状況（2011年9月14日13:15撮影）

40
表層（水深1m）

浮遊物質量[SS](㎎/L）

3530

平水時
出水時2520151050

図-4 流速移動観測による表層付近の沿岸流速ベクトル

St.1-1 St.1-2 St.1-3 St.1-4 St.2-1 St.2-2 St.2-3 St.3-1 St.3-2 St.3-3 St.3-4
測点

図-5 SS の分布（平水時・出水時）

恒流成分は 0.2m/s とした。
透明度板により目視観測した透明度の３倍として示して
写真-1 には 9 月 14 日に撮影された鵡川・沙流川の出水時における濁水の流出状況を示す。濁水は沿岸流に影いる。また、植物プランクトンの増殖に必要とされる栄養塩の閾値は窒素で 0.014mg/L、リンで 0.003mg/L とし響されて北西(写真下方向)に偏向して広がっており、その範囲は数 km に及んでいる。河川流量のピークから 10（例えば高橋ら 10））、それ未満は栄養塩枯渇と考える。
平水時（図-6(1a)(1b)(1c)）では、水温・塩分の躍層日以上が経過しても海域に濁りが広がっており、多量の土砂が陸域から供給されていることが分かる。このようについては、St1-2では水深4m付近、St.2-3では水深5m付に当海域では長期間の濁水の流入があることが分かる。
近、St.3-4では水深7m付近にあり、躍層の上下で水塊構
図-5 には表層（水深 1m）における測点毎の浮遊物質造が異なっている。また、補償深度は海底面より下であ（SS）の濃度を示す｡ここで平水時は 8 月 3 日、出水時り光量の不足は無い。次に栄養塩のうちリンについては、は 9 月 10 日のデータを示している。全体的に出水時が上述した閾値を全測点・全層で上回っており枯渇は無い。
大きく、最大値は河口に近い St.1-1・St.1-2 の 37mg/L
一方、窒素についてはDINとして見た場合にSt.1-2ではで、河口から離れるほど SS 濃度が小さい。また、河口全層で閾値を上回っているが、St.2-3・St.3-4では躍層から離れるほど出水前後の差異が小さくなる傾向にある。 の上部水塊で枯渇し、これによりクロロフィルａが低い。
(2)河川出水が栄養塩と基礎生産に及ぼす影響
躍層の下部水塊では、DINの枯渇が無いことからクロロ
図-6 には、St.1-2・St.2-3・St.3-4 の 3 点における河フィルaが上部水塊よりも高い。よって、平水時では躍川出水前後の水温・塩分・濁度・クロロフィルａ（以上層を境に上部水塊で窒素枯渇となり基礎生産の制限要素
は CTD 計測による、図中の実線）および窒素・リン・ケとなっている。
イ素・鉄・クロロフィルａ・SS（以上は水質分析による、一方、出水時（図-6(2a)(2b)(2c)）は、河川出水の影図中のプロット）の鉛直分布を示す。ここで窒素につい響により表層の塩分の低下が見られる。また栄養塩の供ては、植物プランクトンが栄養塩として利用可能な硝酸給があり、窒素・リンともに全測点・全層で閾値を上回
態窒素・亜硝酸態窒素・アンモニア態窒素の合計であるり窒素では枯渇が解消されている。しかし、SSの供給が溶存態無機窒素を DIN として示している。リンについてあることから、平水時に比べて濁度が増加し、特に表層は溶存態無機リンとして PO4-P を示している。濁りの目で顕著である。これにより補償深度はSt.1-2で0.6m、安としては濁度と SS を示している。植物プランクトン
St.2-3で4.5m、St.3-4で6.0mとなり、底層では十分な光の増殖に必要な光量が確保できる水深である補償深度は、 量が確保されていない。このため、クロロフィルaは表Kenya Sudo , Masami Ohashi , Jun Yamamoto

0.00

0.2

0.30.40.00.50012133440.10.20.35

補償深度 =6.0m

670.40.5

平水時平均値, 490,419↓500,000

991010234St.1-1

補償深度 = 0.6m

50123450102030
4050
水温(℃),塩分(‰),濁度(FTU),SS(mg/L)(1a)平水時(St.1-2)(2a)出水時(St.1-2)

NO3-N,NH4-N,PO4-P,DIN(mg/L)
0.01 0.02 0.03 0.04 0.05

0.00

SiO4-Si,T-Fe(mg/L)
0.2
0.3
0.4

0.0

0.0
00.1

0.01

PO4-P(mg/L)
0.02 0.03 0.04

0.05

NO3-N,NH4-N,T-Fe,DIN(mg/L)
0.5

0.1

0.2

0.3

0.4

0.501223344

水深（ｍ）

156

補償深度 = 4.5m

56778899

補償深度 = 27.0m

10100123chl-a(μg/L)

0102030

454050

水温(℃),塩分(PSU),濁度(FTU),SS(mg/L)

012345
SiO4-Si(mg/L),chl-a(μg/L)

0100.00

0.2

0.3

4050

0.01

PO4-P(mg/L)
0.02 0.03 0.04

0.05

NO3-N,NH4-N,T-Fe,DIN(mg/L)

SiO4-Si,T-Fe(mg/L)
0.130(2b)出水時(St.2-3)

NO3-N,NH4-N,PO4-P,DIN(mg/L)
0.01 0.02 0.03 0.04 0.05

0.020

水温(℃),塩分(PSU),濁度(FTU),SS(mg/L)(1b)平水時(St.2-3)

0.00

0.40.00.501122334455

66水深（ｍ）

078

0.1

0.2

0.3

0.4

0.5

補償深度 = 6.0m

7899101011111212

補償深度=33.0m

St.1-2

St.1-3

St.1-4St.2-1

St.2-2測点

St.2-3St.3-1

St.3-2

St.3-3St.3-4図-7 測点毎の植物プランクトン細胞数の分布

010
2030
4050
水温(℃),塩分(‰),濁度(FTU),SS(mg/L)

0.00

出水時平均値,67,348 ↓6chl-a(μg/L)

水深（ｍ）

1,000,000

81

出水時

070

平水時

植物プランクトン細胞数は鉛直方向平均値

1,500,000

58

水深（ｍ）

2,000,000

0.05

2水深（ｍ）水深（ｍ）

0.1

PO4-P(mg/L)
0.02 0.03 0.04

NO3-N,NH4-N,DIN(mg/L)

SiO4-Si,T-Fe(mg/L)

0.00.01

植物プランクトン総数(個/L)NO3-N,NH4-N,PO4-P,DIN(mg/L)
0.01 0.02 0.03 0.04 0.05

0.00

層で大きく、底層では小さい。よって、出水時では栄養塩は充分足りているが、SS の影響により底層で光量不足となり基礎生産の制限要素となっている。
図-7 には測点毎の植物プランクトン細胞数（鉛直方向平均値）を示している。全体的には平水時が出水時より細胞数が多く、全測点の平均値で比較すると平水時（490 千個/L）が出水時（67 千個/L）の 7.3 倍である。
また、河口からの距離が大きいほど平水時と出水時のプランクトン細胞数の差の割合が小さい傾向がある。これは河口から離れるほど濁水の影響が小さくなるためと思われる。
(3)漁場環境の変動機構
以上の現地観測結果に基づき、漁場環境の変動機構の概要をまとめると図-8 のとおりとなる。
河川流量が少ない(＝栄養塩類供給が少ない)平水時には、植物プランクトンが栄養塩（窒素・リン等）を摂取して光合成することで増殖し、基礎生産量の指標であり光合成色素であるクロロフィルａが増える(図中①、以下同様)。栄養塩が豊富な下部水塊では基礎生産が維持されるが、上部水塊では植物プランクトンの消費により枯渇する(②)。植物プランクトンは、底生生物である二枚貝の主要な餌料として摂餌される(③)。
次に、河川流量が増える出水時には、栄養塩とともに多量の土砂が海域に流れ込む。粒径の小さい浮遊土砂は濁水となって浮遊物質濃度（SS）を高め(④)、太陽光を遮断して光環境が悪化する(⑤)。これが進むと植物プランクトンの光合成が阻害されて基礎生産量が低下して二枚貝の餌料環境が悪化する(⑥)。また、濁水は二枚貝の体力を消耗させて斃死に至る場合があるなど生息環境を悪化させる(⑦)。

131314140101023chl-a(μg/L)
2030454050

水温(℃),塩分(PSU),濁度(FTU),SS(mg/L)

001234SiO4-Si(mg/L),chl-a(μg/L)
102030405平水時
50

水温(℃),塩分(PSU),濁度(FTU),SS(mg/L)
栄養塩類：窒素･リン等(1c)平水時(St.3-4)
○
●
▲

NH4-N
NO3-N
chl-a（採水）
水温補償深度□◇×

出水時
太陽光(2c)出水時(St.3-4)
DIN
PO4-P
SS
塩分濁度△
◆

SｉO2-Si
T-Fe

魚類②摂取による栄養塩枯渇①植物プランクトン増殖→クロロフィルａ増加

④浮遊土砂
SS増大
光合成の阻害

摂餌
動物プランクトン増殖

光合成

排泄･沈降

Chl-a(Obs.)※chl-a（採水）は吸光光度法，chl-a(Obs.)は蛍光光度計による測定

二枚貝

③摂餌

⑤濁水による太陽光の遮断
××××××××××××

太陽光

枯死･沈降

河川
出水沈降
⑥植物プランクトンの枯死→基礎生産の低下→餌料環境の悪化

枯死･沈降←貝類が摂餌 ←沈降堆積したデトリタス

図-6 水温・塩分・栄養塩・クロロフィルａの鉛直分布

Kenya Sudo , Masami Ohashi , Jun Yamamoto

図-8 漁場環境の変動機構⑦二枚貝の斃死

6. 数値解析(1) 物理特性の再現11)
前述の現地観測結果に基づいて、数値解析により物理特性の再現を試みた。
鵡川海域の物理特性の再現には、マルチレベルモデル4)を用いた。計算は前述の図-2に示す12.6km×22.1kmの範囲で行い、水平方向を格子数221×126で100mメッシュ、鉛直方向を19層(可変)とした。計算期間は前述の図-3に示す2011年9月1日から9月15日までとし、台風12号による出水の影響を網羅した。開境界条件のうち、潮位は苫小牧港の推算潮位を用い、恒流成分は前述の図-4に示した流速移動観測結果に基づいて沿岸方向に0.2m/secの流れを与えた。陸側境界条件は河川流量の実測値を与えた。
また、SSはL-Q式により与えた。
図-9には、前述の現地観測結果(図-6)と同様に河口から西に沖合に向かってSt.1-2・St.2-3・St.3-4の3点における水温・塩分の計算値と観測値の比較結果を、出水時(9月10 日) について示す。水温については計算値と観測値はほぼ一致している。塩分については河口に近いSt.1-2では計算値と観測値の一致が見られる。しかし、
St.2-3とSt.3-4の表層付近の塩分低下を計算が再現しき
図-10 SSの再現結果画像れていないと考えられ、今後の課題である。解析範囲全(上図:2011年9月4日4:00･下図:9月14日12:00)
体の傾向としては、より河口に近い箇所について良好に再現していると考えられる。
(2) 基礎生産量の試算12)
図-10には、SSに基づいた濁水の挙動の数値解析によ
前述の現地観測結果（図-6）に基づき、表-2に示す条る再現結果画像を示す。計算期間のうち、河川出水のピ
件で低次生態系モデル5) を用いてクロロフィルａの現況ークである9月4日4:00および前述の写真-1の空撮を行っ再現および基礎生産量の試算を行った。測点はSt.2-3とた時刻に近い9月14日12:00の画像を示している。SSの最大値を70mg/Lとし、色彩が赤色に近いほどSSの値が高く、 し、平水時から出水時に変化する過程の9月1日から10日の9日間（図-3参照）を計算期間とした。ここで9月1日
緑色から青色に近づくほどSSの値が低くなることを示しの初期値として平水時(8月3日)に観測したクロロフィルている。9月4日では鵡川・沙流川の河口付近で高いSSのａを用い、栄養塩および消散係数(透明度に基づく)は出赤色を示し、沖合に向かって緑色から青色に変化してSS水時(9月10日)の観測値を用いた。
が低下しており、前述の現地観測結果（図-5）の傾向にはじめにクロロフィルａの再現計算結果を図-11に示合っている。また、沿岸流により濁水が北西方向に偏向す。躍層のある5～7m付近を境に、平水時に表層で窒素して拡散する様子が再現されている。9月14日では9月4制限の影響を受けていたクロロフィルaが、出水時に底日に比べて河口付近のSSの値は低いが、9月4日と同様の層で光制限を受けた状態へと変化する過程を計算値が良SSの分布傾向を示している。前述の写真-1と同様に濁水好に再現している。の拡散範囲が数kmの沖合まで広がっており、広範囲に渡次に、表-2に示す計算条件に基づいて基礎生産量を試り二枚貝の生息環境に影響を及ぼしていることが示唆さ算する。基礎生産量は、クロロフィルaを日当り(/day)
れる。
単位水塊当たり（/m3）の炭素量（mg-C）に変換して表示する。ここで変換式はC=50×chl-aを用いた。計算ケースは、(A1)現況再現、(A2)現況A1に対して濁水対策を行った場合、(B1)現況のまま濁水による光制限がさらに9日間継続した場合、(B2)出水時B1に対して濁水対策を
2011.9.10 06:0000:00
2011.9.10 06:00
00:00
2011.9.10 08:00
00:00
St.1-2
St.2-3St.3-4行った場合の計4ケースである。ただし、B1･B2のクロロ
フィルａの初期値は9月10日とする。図-12 (a)は計算ケ【凡例】：水温(観測)○ ： 水温(計算)：塩分(観測)○ ：塩分(計算)ースA1･A2の結果を示す。濁水の影響がある現況(A1)で
図-9 水温･塩分の現況再現結果(2011年9月10日)
は、水深6m以深で基礎生産が負になっており、前述のWater temperature(℃)

Water temperature(℃)
51015200255101520

Water temperature(℃)

25002224446810Depth(m)0Depth(m)Depth(m)00681012510152025

Salinity (psu)3035152025681214010101214514051015 20

25

Salinity (psu)

30 350510 15 20 25 30 35
Salinity (psu)

Kenya Sudo , Masami Ohashi , Jun Yamamoto

表-2 計算条件 (St.2-3)
表-2 計算条件 (St.2-3)
chl-aの
初期値

Case名

栄養塩

A1(現況再現) 平水時（8/3）
A2
B1
B2

出水時(9/10)の栄養塩
出水時（9/10）
を供給出水時（9/10）

平水時（8/3）

7. おわりに濁水対策 消散係数※1の有無
ks

計算時間※2無し

1.133

9日(9/1～9/10)

有り

0.189

9日(9/1～9/10)

無し

1.133

9日

有り

0.189

9日

※1：消散係数は透明度(Ds)を用いて算出した．ks=1.7/Ds※2：9日間の計算を行い初期値からの1日あたりの変化量を基礎生産量として表示

01chl-a(μg/L)
2340
1
2
水深（ｍ）

34
5678910

観測値(平水時8/3)
初期値(計算点；平水時8/3）
観測値（出水時9/10）
再現計算結果(A1；9日後)基礎生産量(mg-C/m3/day)-100102030012345678910
A1:現況再現
A2水深（ｍ）

水深（ｍ）

図-11 クロロフィルａの再現結果(St.2-3)
基礎生産量(mg-C/m3/day)-100102030012345678910B1B2

本報では、漁場環境の変動機構を現地観測で捉え、数値解析で再現・試算することで、陸域の影響を適切に評価することが可能となった。以下にまとめを示す。
①研究の目的と現地条件を考慮した適切な現地観測計画の立案により、河川出水時の沿岸漁場環境の変動機構を解明した。
②現地観測結果を適切に再現する数値解析により、濁水対策を行った場合の基礎生産量の試算を行った。この
結果、濁水対策が二枚貝の餌料環境を改善し、生産性を向上させる可能性があることが示された。
今後は、河川出水時における沿岸漁場環境の変動機構解明を継続するとともに、二枚貝等の個体群動態モデルを構築する予定である。将来的には本方法である現地観測や数値解析の手法を、問題を抱える他の沿岸漁場や港湾・漁港の現場へ適用し、解決に繋げたいと考える。
参考文献
１）水産庁(2012)：平成 23 年度水産白書
２）真岡東雄(1993)：鹿島灘砂浜域における漁業生産の特徴と問題点、水産海洋研究、第 57 巻 第 2 号、pp.32-38
３) 須藤賢哉・大橋正臣・山本潤(2012)：水産水域環境の諸問題に対応した水質変動予測手法について、北海道開発技術研究発表会論文、環 29４）平野敏行監修(1998）：沿岸の環境圏、フジ・テクノシス
テム、pp.674-697
５) 中田喜三郎 (1993)：生態系モデル―定式化と未知パラメ(a)初期値：平水時8/3(b)初期値：出水時9/10

図-12 基礎生産量の試算結果(St.2-3)ータの推定法―，Journal of Advanced Marine Technology
Conference、 Vol.8, pp.99-138
６）北海道新聞社（2001）:9 月 18 日付け北海道新聞朝刊苫小図-6(2b)で示したような底層での光制限の影響を再現している。濁水対策を行った場合(A2)は基礎生産量が改善し全層で正となり、最大で約20mg-C/m3/dayとなって光環境の改善を示している。図-12(b) には計算ケースB1・B2の結果を示す。現況の濁水の影響が継続した場合(B1)では水深2m以深で基礎生産が負であるが、濁水対策を行った場合(B2)は全層で正となり、基礎生産量が改善している。
このことから河川出水が長期間続いた場合、中・底層で光環境が悪化し、沿岸漁場の基礎生産量を低下させていると考えられる。濁水対策を行うことで基礎生産が回復し、二枚貝の餌料環境が改善して生産性が向上する可能性があることが示された。
以上より、基礎生産量を指標とした出水前後の現地観測とこれを適切に再現する数値解析モデルにより、現況の基礎生産量を把握するとともに、現地における濁水対策の有効性について定量的に試算・検討することが可能であると考える。

Kenya Sudo , Masami Ohashi , Jun Yamamoto

牧版、23 面７）瀬戸雅文･服部志穂･新居久也(2005):ホタテガイの濾水浄化能力に関する基礎的研究、海洋開発論文集、第 21 巻、pp.247-252
８）山下俊彦・宮下将典・山崎真一・渡邊康玄 (2000)：河川から供給される物質の河口沿岸域での挙動、海岸工学論
文集，第 47 巻、pp.1026-1030
９) 北海道開発局室蘭開発建設部苫小牧港湾建設事務所(1997):平成８年度苫小牧港海域調査検討業務報告書、118p.
10）高橋正征・古谷 研・石丸 隆(1996)：生物海洋学２「粒状物質の一次生成」、東海大学出版会、174p．
11) Masami Ohashi, Jun Yamamoto, Kenya Sudo(2012) ： A
Numerical Study on Influence of Primary Productivity byRiver Flood in Coastal Region, Pacific Congress onMarine Science and Technology (投稿中)12）大橋正臣･山本潤･須藤賢哉･水垣滋･門谷茂･田中仁(2012):
鵡川沿岸の漁場環境に及ぼす河川出水の影響､土木学会論文集 B2(海岸工学)､Vol.68､No.2、pp.1121-1125