平成23年度

水産水域環境の諸問題に対応した
水質変動予測手法について（独）土木研究所寒地土木研究所

水産土木チーム

○須藤大橋
山本賢哉
正臣
潤

水産利用している沿岸漁場や漁港水面で水域環境の改善を図る場合、対象水域の水質悪化メカニズムを把握し、物理特性・水質特性を定量的に評価する必要がある。この評価の手法は、詳細な現地観測を行った上で、数値シミュレーションを用いて時間的・空間的な水質変動予測を行うものである。これにより水質改善目標及び対策の設定と対策後の予測評価が可能となる。
本報告は、水産水域環境の悪化が懸念される海域（閉鎖性海域、河川流入負荷のある海域）における水質変動予測手法について、現地観測と数値シミュレーションの両面から検討した結果を述べるものである。
キーワード：水産水域環境、水質変動予測、密度流計算、生態系モデル

1． はじめにを進めてきた 3)4)5)。本報告は、閉鎖性海域、河川流入負荷のある海域など水産水域環境の悪化が懸念される水域における水質変動予測手法について、現地への適用事例とその結果を併せて総体的に説明するものである。

近年我が国では、磯焼け等による藻場喪失や内湾域の水質悪化、海水温の上昇など海域環境の変化に伴い、漁業生産量がピーク時の半分程度に減少しており、水産資源の回復が大きな課題となっている。そのような情勢の2．水産水域環境における問題点
中で、北海道は漁業生産量が国内の約１／４を占め 1)、漁業生産拠点として水産物の安定供給に大きく貢献して漁業生産拠点となる内湾や漁港などの閉鎖性水域ではいる。北海道において将来にわたり安定的な漁業生産を海水の交換・循環が少ないため、水質・底質の悪化が生維持するためには、沿岸漁場の豊かな水域環境に配慮しじやすい。閉鎖性の水産水域環境における問題点は、水つつも、内湾や漁港を活用した蓄養、中間育成、養殖な域内への汚濁原因物質の蓄積と、それに伴う溶存酸素ど水面の高度利用をより一層推進していくことが必要で（ＤＯ）の消費に着目することで説明できる。図-1 に2)ある 。これらの整備は、施設周辺の水域環境の維持・水質・底質悪化のしくみの概略を示す。閉鎖性水域に河川水や養殖等に由来する汚濁負荷があると（図中①②、
保全と一体で進める必要があり、そのためには適切な水以下同様）、海水中の窒素・リン等の栄養塩の濃度が高質特性の評価が不可欠である。
まるとともに底質中に有機物が蓄積する。これは富栄養一方で、漁業生産拠点となる内湾や漁港では、外海との海水交換が不足して水域内に汚濁原因物質が滞留しやすく、水質悪化に伴う水産有用種の生育不良や大量斃死等の生産性低下に繋がる問題が見られる。加えて、河川流入負荷の大きい水域では、都市由来の生活雑排水、農業由来の肥料・糞尿、河岸・森林由来の土砂などが上流から海域に流れ込み、河口周辺海域の水質が大きな影響を受ける。将来にわたって良好な水産水域環境を保全し、生産性の維持・向上を図っていくためには、各水域における流動機構・水質変動機構を解明した上で、将来の事業効果を事前に予測し、評価しておくことが必要となる。
このような水産水域環境の諸問題を解決するために、我々は閉鎖性海域や河川流入負荷のある海域を対象とし
図-1 水質・底質悪化のしくみの概略て、水域の水質特性の解明と水質変動予測について研究Kenya Sudo , Masami Ohashi , Jun Yamamoto

水域環境改善の対策を行う上では必要不可欠な手法と位化と呼ばれ、過密養殖あるいは上流域に都市や農業地帯置づけることができる。
がある河川の流入により生じやすい。例えば後述の風蓮
水質変動予測では、現地観測により個々の海域の水質湖の事例 3)では河川からの汚濁負荷によりＣＯＤ（化学6)特性を明らかにした上で、数値シミュレーションにより的酸素要求量）が環境基準を超えている 。富栄養化が物理特性、栄養塩やＤＯ等の水質を計算するという手順
進むと、植物プランクトンとそれを捕食する動物プランをとる。計算の再現性を高めるためには、海域特クトンが大量発生する（③）。それらの死骸が有機物と性の正確な把握が不可欠であるため、流況・水質して底質中に蓄積し（④）、それが好気性細菌により酸化分解される過程で溶存酸素が大量に消費される（⑤）。 等に関する詳細な現地調査を行う。
水質変動予測手法の手順は以下のとおりである。
溶存酸素は動物プランクトンや魚介類によっても消費さ(1) 手順１：現地調査による海域特性の把握れる。すると溶存酸素が極端に少なくなり、貧酸素水塊①連続観測：定点設置型の観測機器により流速・波高・が形成される（⑥）。この貧酸素化による魚介類の窒息
水温・塩分・ＤＯ・クロロフィル等を計測する。
が、内湾での養殖や漁港水面での蓄養における大量斃死
②多点観測：垂下式機器により水温・塩分・クロロフィの直接的な原因である。例えば後述の能取湖 5)では貧酸
ル等の鉛直分布を計測する。
素水塊によるホタテ稚貝の大量斃死が生じている。貧酸③試料採取：採水・採泥により栄養塩・ＣＯＤ等に関し素化が進んだ底質においては、蓄積された有機物が嫌気て水質・底質を分析する。
性細菌（硫酸還元菌）により分解される過程で硫化水素④気象データ収集：気温・風速・日射量等を把握する。
が発生する（⑦）。硫化水素が溶存酸素により酸化され(2) 手順２：密度流計算による物理特性の再現て硫黄となり、海面付近で白濁して析出するのが青潮で密度流計算は、海洋モデルを用いている。これは一般ある（⑧⑨）。長年にわたって養殖・蓄養を行っているに「プリミティブモデル」と呼ばれ、非圧縮粘性流体、内湾や漁港では底質悪化が常態化しており、青潮が発生
ブシネスク近似、静水圧近似を基礎としたモデルであり、しやすい状況にある。青潮は全国各地の内湾で見られ、ＰＯＭやＭＥＣなど様々なモデルが公開されている。こ
例えば近年の東京湾では年３回程度発生している 7)。ま5)のモデルを用いれば、流速、水位、水温、塩分についてた、後述の能取湖 でも養殖による底質悪化が進んでお空間（３次元）的、時間的な変化の解析が可能となる。
り、2007 年 9 月に初めての青潮が確認されている。
なお、密度流計算の概念とその計算方法の詳細につい以上のように、水質・底質の悪化は様々な要因が相互ては文献11)を参照されたい。に関連して生じる。水産水域環境を改善するためには、水域における栄養塩等の汚濁原因物質の挙動を把握した(3) 手順３：低次生態系モデルによる栄養塩等の再現上で、ＣＯＤ・ＤＯ等の指標に着目することで適切に水低次生態系モデルは、生態系を構成している生物(植質特性を評価する必要がある。
物プランクトンや動物プランクトンなど)と非生物(無機態栄養塩や有機物など)を構成要素とし、主に食物連鎖等に基づいて、これらの構成要素間の相互作用を計算す
3．水産水域環境改善と 水質変動予測の考え方るものである。これにより栄養塩、ＤＯ、ＣＯＤ等の３次元分布とその時空間変化を求めることが出来る。
水産水域環境の改善を図ろうとする場合、対象とする
なお、低次生態系モデルの概念とその計算方法の詳細
水域の水質悪化の特性とその要因を把握し、物理特性・については文献11)を参照されたい。
水質特性を定量的に評価する必要がある。その上で水質改善目標を定めたり、必要に応じて水質改善対策や水域現地調査に数値シミュレーションを組み合わせた水質の利用調整を行ったりする。
変動予測手法のメリットは、ある特定の日時・海域条件
水産庁による「自然調和型漁港づくり技術マニュアのもとでの数値計算結果が現地調査結果を整合性良く再8)ル」 に基づいて、水域の水質悪化状況の把握から水質現できた場合、調査地点以外の任意点（例えば現地計測不可能な場所など）での値や海水交換率等の現地計測不改善に向けての手順を示すと概ね以下のとおりとなる。
可能な項目を再現できることである。また、調査対象①対象水域の自然条件、水質悪化状況を把握する。の水域において構築された数値計算モデルは、同②対象水域に求められる水質環境を環境基準9)・水産用
10)様の水質悪化問題を抱える別の現場に応用するこ
水基準 に基づき設定する。とが可能である。
③①で得られたデータに基づき数値シミュレーションを用いて水質変動予測を行う。
④③で得られた結果が②の基準を達成しているか否かを4．風連湖における水質変動予測の事例判定する。
⑤以上①～④で得られた成果を水産水域環境改善の対策(1) 概要に反映させる。
調査対象は、図-2 に示す北海道東部の根室湾に位置
水質変動予測手法は上記の①③の部分に相当し、水産する風蓮湖である。風蓮湖は、国内第 14 位の湖水面積
Kenya Sudo , Masami Ohashi , Jun Yamamoto

N

根室湾

矢臼別川St.1 北奥部

走古丹漁港（実験場所）

St.3 流入部
風蓮川

St.6 狭窄部 （採水地点）
ADCP設置2007年

風蓮湖 St.4 湖口部

St.5 湖口外ADCP設置2006年
St.0
風蓮橋

別当賀川

5km

図-2 調査位置図（風蓮湖）(2) 現地調査
現地観測は、降雨による河川出水期を捕捉するため2006、2007 年の夏から秋にかけて風蓮湖において実施した。図-2 に示す風蓮湖において、流況、水質、底質、プランクトン、光合成速度について現地調査を行った。
風蓮川からの汚濁負荷を確実に捉え、かつ河川水が風蓮湖を経由して根室湾に出る経路を網羅する形で St.0～
St.6 の観測地点を設定した。
観測結果では、風蓮湖の湖口付近は清浄な外海水との海水交換によって水質が良好に保たれるが、湖奥部では外海水の影響が届かず、河川水による汚濁負荷によって水質悪化を引き起こしていることが分かった。図-3 には 2006 年におけるＣＯＤと塩分濃度の関係を示す。ほぼ淡水の St.3 流入部(左上点線囲い)では降雨後の 11 月に顕著に増加しているが、他の St.1･5･6（右下点線囲い）では塩分濃度と負の相関が高い。つまり、風蓮湖に

Kenya Sudo , Masami Ohashi , Jun Yamamoto

20.0
CODMn (mg/L)（57.5km2）を持つ汽水湖であり、主な流入河川は、風蓮川・別当賀川・ヤウシュベツ川である。流域は大規模な酪農地帯であり、農業由来の汚濁負荷が流入し、水質及び底質の悪化が懸念されている 6)12)。風蓮湖は環境基準 9)の類型指定水域でありＣＯＤ＝5.0mg/L と設定されているが、この基準値を超過して問題となっている 6)。
このような問題を受けて、良好な河川・沿岸環境の保持・再生と農業の持続的な発展の両立が重要な課題となっている。
風蓮川流域では酪農由来の汚濁負荷を抑制する対策（沈砂池、緩衝林帯）が検討されている。その施策を推進するためには、対象とする水域において現地観測を行い、流動・水質変動機構を解明した上で、将来の事業効果を事前に予測し、評価しておく必要がある。水産土木チームでは2006年より風蓮湖において流動・水質特性等の詳細な現地観測を行い、水質変動予測が可能な数値計算モデルを構築した。これによる陸域からの汚濁負荷の流出抑制対策の前後における水質変動予測の結果を以下に示す。

8月
9月
10月
11月

y = -0.2226x + 10.374
R2 = 0.7134

15.010.0
5.00.005101520
塩分濃度 (‰)

253035図-3 ＣＯＤと塩分濃度の関係

おける汚濁負荷の高い河川水と清浄な外海水との混合割合は、塩分濃度を用いることで概ね説明できることが分かった。
(3)水質変動予測結果
現地観測結果及び光合成速度現地試験結果に基づく数値シミュレーションによりＣＯＤ収支計算を実施した。
計算手法に関する詳細は文献11)を参照されたい。
数値シミュレーションのケースは、表-１のように①
平水時の現況を現地観測結果に基づいて再現したもの（平水時現況再現）、②平水時の現況において流量を出水時のものに変更したもの（出水時対策なし）、③出水時において栄養塩類等の負荷量を低減する対策を行ったもの（出水時対策あり）の３ケース行うこととした。対象期間は表中のとおりとし、最大流量は対象期間の観測値である。河川からの流入負荷量は、寒地土木研究所水環境保全チームが作成した全窒素・全リン・ＳＳのＬＱ式13)にこれらの流量を代入して算出し、Ｌ-Ｑ式の無い水質項目については、今回の観測結果の統計量から得られる比率を用いて配分し、それぞれの値を上流側の境界条件として与えた。ケース③では、出水時の河川流量のまま、負荷量低減率を前述の文献13) のＬ-Ｑ式に基づき0.5とした。
図-4は表層のＣＯＤの平面分布を示す。縦軸・横軸は距離（単位：ｍ）である。①と②を比較すると、ＣＯＤの濃度は河川流入部付近で高く、対策前の②では当水域の環境基準値５mg/Lを大幅に上回っていることがわかる。対策後の③では流入部のＣＯＤが低下しており、流出抑制対策の効果が現れている。
図-5は湖内各地点におけるＣＯＤの経時変化を示す。
①の平水時では、河口付近のSt.3の経時変化の傾向は、河川流量の経時変化との相関が見られた。その他の地点は時間的に一定周期の変動を示し、河川流よりも潮流の影響が強いことが伺われた。出水時における②対策なし
表-1 数値シミュレーションの計算ケース
表－7.1 計算ケース
計算ケース

対象期間

最大流量

負荷量低減率①平水時現況再現

2007.09.18～09.22

32.41m3/s

1.0

2007.09.07～09.15

132.91m3/s②出水時対策なし③出水時対策あり

1.00.5

と③対策ありを比較すると、いずれも河口付近のSt.3でＣＯＤが高い。③対策ありではＣＯＤの大幅な低下が見られ、かつ出水の最大流量時にその差が最大となっており、汚濁負荷流出抑制対策の効果が顕著に現れていることがわかる。
以上の結果より，現地観測により水質等の現況を把握し，河川流量データを条件として与えて計算することで，水質汚濁負荷物質の挙動の再現・予測を行った事例を示した。

① 平水時（H19.09.07～H19.09.15）(m)② 出水時対策なし（H19.09.07～H19.09.15）(m)①平水時：計算開始 24 時間後 9 月 19 日 11:00③ 出水時対策あり（H19.09.07～H19.09.15）

図-5 各測点におけるＣＯＤの経時変化(m)

5．能取湖における水質変動予測の事例(m)②出水時対策なし：計算開始 84 時間後 9 月 10 日 23:00(m)(1) 概要調査対象の能取湖は、図-6に示すように北海道北部のオ ホ ー ツ ク 海 に 位 置 す る 国 内 第 13 位 の 湖 水 面 積（58.5km2）を持つ海跡湖である。湖口の幅・水深が非常に小さくて静穏度が高いため、水域ではホタテ稚貝の養殖が盛んに行われ、種苗として全道に供給されている。
しかし、長年の養殖による底質悪化と閉鎖性の高さから、夏季に溶存酸素（ＤＯ）が不足して貧酸素水塊が発生しやすい状況である。2007年9月19日には，能取湖では初めての青潮が湖内南岸の卯原内地区で発生し、中間育成

能取湖

(m)
③出水時対策あり：計算開始 84 時間後 9 月 10 日 23:00

図-4 ＣＯＤ平面分布の計算結果

Kenya Sudo , Masami Ohashi , Jun Yamamoto

図-6

調査位置図(能取湖)

図-7に湖央部（st.3）における水温・ＤＯの経時変化、中のホタテ稚貝のうち、約1,200万粒の斃死被害が報告
気温経時変化、風速経時変化（風向は360°を北、180°されている。また、2008年の9月20日に能取湖北岸の定を南）を示す。水温（①上）については、8月26日以降置網で漁獲されたサケ・マスに貧酸素水塊の影響と見らに水温差が表層と中層で3℃程度、中層と底層で5℃程度れる斃死被害が生じた。となり、その状態が9月中旬まで続いた。9月19日頃より以下では、水質変動予測手法により能取湖における貧底層の水温が上昇し始め、一方で気温の低下により表層酸素水塊の挙動の解明を行った事例を示す。
水温が低下したため、9月23日以降では各層の温度差が(2) 現地調査ほぼ消滅した。ＤＯ（①下）については、8月下旬から図-6に示す能取湖の各地点において、2008年8月から水温の成層化に合わせて底層のＤＯが徐々に低下し、910月にかけての夏季の貧酸素化が生じやすい時期を狙って、流況、水質、底質等の観測を実施した。連続観測は、 月11日頃には 2mg/L 程度まで低下して貧酸素水塊の形成が認められた。底層のＤＯは9月19日頃から回復傾向
水温・塩分計、ＤＯ計、ADCPによって行い、調査地点はを示し、気温の低下（②）及び9月19日頃に観測された湖底の窪地深部における貧酸素水塊の発生とその挙動の北向きの強風（③）により鉛直混合が促進したと推察さ全体像の捕捉を意図して、湖口付近（st.1）、湖口側深れる9月23日にはほぼ平常値の8mg/L程度まで回復した。
部（st.2）、湖央部（st.3）、湖奥部（st.4）に設定し(3) 水質変動予測結果た。加えて、採水・採泥による分析、ＳＴＤ及びＤＯ計
貧酸素水塊の挙動を検証・予測するため、現地観測結による鉛直分布観測を地点ア～クにおいて実施した。ま果に基づく数値シミュレーションによりＤＯ収支計算をた、測線L1～9ではADCPを取付けた調査船を走行させ，実施した。計算方法の詳細は山本らの事例14)に倣ってい湖内流況の３次元分布を把握した。さらに，低次生態系モデル計算における生物パラメーターを取得するため、るのでそちらを参照されたい。
自記式ＤＯ計を入れた明暗筒により，現地における酸素図-8には湖央部st.3の水温・ＤＯの観測値と再現計算消費速度の連続計測を8月、9月の各３日間実施した。
値の比較を示す。計算値はほぼ観測結果と同様の変動を示し，水温変化とＤＯ収支が良好に再現されている。
表層(-1m)中層(-10m)
底層(-19m)
水温(℃)
図-9にはＤＯの計算結果の３次元分布を示す。図中の25ｓｔ．３
黒い部分が貧酸素水塊である。成層状態で底層が貧酸素20化している9月16日0:00を計算開始時刻とした。9月19日
1510
未明に観測された強い南風（図-7③）により、図-9①で8/20 8/25 8/30 9/4 9/9 9/14 9/19 9/24 9/29
は表層水が風下である北岸に移動し、その代わりに底層
表層(-1m)中層(-10m)
底層(-19m)
DO(mg/L)の貧酸素水塊が反対に風上である南岸の湖奥部付近に移ｓｔ．３
10動していることがわかる。その後，9月19日午後に風速
510m/sを超す強い北風が観測されており（図-7③）、図上下層の混合
09②ではこれに伴って南風時と逆向きの同様の現象が生
8/20 8/25 8/30 9/4 9/9 9/14 9/19 9/24 9/29
じて貧酸素水塊が北岸に移動していることがわかる。こ①湖央部 st.3 における水温(上)・ＤＯ(下)の経時変化の現象は、9月20日の6:00に北岸の定置網においてサ（気象庁網走気象データ時別値使用）
ケ・マス200～300kgが斃死したことと一致する。図-1030 気温(℃)
はサケ定置網設置箇所のＤＯの計算値を示している。強20い北風が観測された9月19日午後以降にＤＯが低下し、
10水産用水基準の6.0mg/L（海域）を下回っており、計算夜間の気温の低下→表層の冷却0
結果は貧酸素水塊による漁業被害の発生を裏付けている。
8/20 8/25 8/30 9/49/9 9/14 9/19 9/24 9/29

水温 (℃)
25②気温の経時変化(気象庁のデータ)

風速 (m/s)
16128

平均風速
強い北風

最大風速

平均風向 風向°
360

強い南風

4270
180
9000
9/16 9/17 9/18 9/19 9/20 9/21 9/22 9/23 9/24
③風速の経時変化(西網走漁協のデータ)
図-7

水温・ＤＯ・気温・風速の経時変化

Kenya Sudo , Masami Ohashi , Jun Yamamotost3 1.0m(観測値)
st3 0.5m(計算値)st3 10.0m(観測値)st3 10.5m(計算値)st3 19.0m(観測値)st3 19.5m(計算値)

20
15
10
9/16

9/17

DO (mg/L)
159/189/19

st3 1.0m(観測値)
st3 0.5m(計算値)

9/20

9/21

st3 10.0m(観測値)st3 10.5m(計算値)9/229/239/24

st3 19.0m(観測値)st3 19.5m(計算値)10509/169/17

9/18

9/199/209/219/229/239/24

図-8 水温(上)・ＤＯ(下)の計算値と観測値の比較(st.3)

サケ･マス斃死発生箇所① 南風継続後②北風継続後(被害発生)

図-9 貧酸素水塊（図中の黒い部分）の３次元分布

これ以後は、図-7①下のとおりst.3のＤＯが回復傾向となり、9月23日には平常値に回復していることも現地観測結果と一致する。
図-11はst.3におけるWedderburn数（Ｗｅ数）の経時変化を示す。Ｗｅ数は上下層の混合程度を表す指標でリチャードソン数／アスペクト比と定義され、抵抗係数ＣD（通常1.3×10-3 ）、海域の長さＬ（8km）、水深ｈ（19m）等により計算した。尹ら15)によれば、Ｗｅ数≦12が青潮発生の目安とされるが、9月19日にＷｅ数が10以下となり、ＤＯが最低値となりサケ・マス斃死が発生した9月20日にほぼ一致する。このことは水域のモニタリングを行い、気象データを用いて被害発生を事前に予測できることを示唆している。
以上の結果より、水温・塩分・ＤＯの観測により密度成層の状況を把握し、気象データを条件として与えて計算することによって、貧酸素水塊の挙動の再現・予測を行った事例を示した。

サケ定置網（表層 0.5m）

サケ定置網（底層 10.5m）

10509/16

9/17

9/18

9/19

9/20

9/21

9/22

9/239/24

図-10 サケ定置網設置箇所の DO 計算値
Wedderburn No.

We数
10000
1000
100
1010.1Δρ(St.3 表層－底層) Δρ(kg/m 3 )42
Ｗｅ数が10 を下回った時間帯

9/1

9/6

9/1109/16

9/219/2610/1（We数は気象庁網走気象データ時別値とst.3観測値を使用）

図-11

参考文献
1)北海道水産林務部総務課(2011)：平成 22 年北海道水産現勢、
概要 pp.1、北海道庁ホームページ－水産統計
2)国土交通省北海道開発局(2011)：北海道の「みなと」と「技術開発」、北海道開発局ホームページ－農業・水産
3)山本潤・牧田佳巳・山下彰司・田中仁(2007)：風蓮湖に陸域からの汚濁負荷が及ぼす影響に関する現地観測、海岸工学
論文集、第 54 巻、pp.1006-1010
4)山本潤・林田健志・峰寛明・牧田佳巳・山下彰司・田中仁(2008)：北方の閉鎖性海域に適した低次生態系モデル構築に向けた現地実験、海岸工学論文集、第 55 巻、pp.1196-1200
5)山本潤・酒向章哲・渡辺光弘・牧田佳巳・田中仁(2009)：能取湖における密度成層崩壊時の貧酸素水塊の挙動に関する現地観測、海岸工学論文集、第 56 巻、pp.1291-1295
6) 北 海 道 環 境 生 活 部 環 境 推 進 課 (2011) ： 北 海 道 環 境 白底層のDOが6mg/Lを下回った時間帯

DO (mg/L)
15等の施策を進める上で必要不可欠な手法であるといえる。
現在では、数値計算とその応用に関する研究開発が進んだため、様々な現場条件に対して水質変動予測が効果的に適用できるようになっている。しかし、その手法が完全に確立したとは言いがたく、現地観測を含めたさらなる研究を進める必要がある。
今後については、閉鎖性海域での研究を進めるとともに、開放性海域等においても水質変動機構の解明を試みる予定である。現在我々は開放性で河川流入負荷のある鵡川海域の事例に取り組んでおり、現地データの収集・解析と数値計算モデルの構築を進めている。将来的には、生態系モデルを確立し、水質変動予測と資源管理を結びつけた手法の確立、および他の海域への応用に結びつけたいと考えている。

Wedderburn 数（Ｗｅ数）の経時変化

6．おわりに
本稿では、河川流入負荷がある海域におけるＣＯＤの悪化事例や閉鎖性海域における貧酸素化の事例に対して水質変動予測手法を適用し、現地状況を良好に再現・予測できることを示した。湾内や漁港の閉鎖性水域において養殖等の漁業生産を行う場合には良好な水質の維持が不可欠であり、水質変動予測は水産水域環境の維持管理
Kenya Sudo , Masami Ohashi , Jun Yamamoto

書’11(2011 年版)、pp.129、北海道庁ホームページ
7)国土交通省関東地方整備局港湾空港部(2011)：東京湾環境情報センターホームページ－東京湾を取り巻く環境（赤潮・青潮）
8)水産庁漁港部(2000)：自然調和型漁港づくり技術マニュアル－水質浄化技術－、pp.3-26
9)環境省(2011)：環境省ホームページ、環境基準（水質）10) (社)日本水産資源保護協会(2006)：水産用水基準(2005 年版)
11)平野敏行監修（1998）：沿岸の環境圏，フジ・テクノシス
テム、pp.674-697
12)北海道環境科学研究センター(2005)：北海道の湖沼改訂版、pp.46-51
13)中村和正・山下彰司・山本潤・中山博敬・大久保天・鵜木啓
二・多田大嗣・林田寿文・渡辺光弘・加藤道生・西山章彦・斉藤勉(2010)：大規模酪農地域における水質保全策が流域の水質に与える効果、寒地土木研究所月報、第 686 号、pp.2-1814)山本潤・田中仁・佐伯信哉(2007)：野見湾の内部潮汐による海水交換と水止まりによる貧酸素水塊発生に関する研究、土木学会論文集 B、第 63 巻、1 号、pp.39-5015)尹 鍾星・中辻啓二・村岡浩爾(1993)：開放性成層水域における吹送密度流に関する研究、海岸工学論文集、第40巻、pp.241-245