水研センター研報，第34号，49−70，平成23年
Bull. Fish. Res. Agen. No. 34, 49‑70, 2011

瀬戸内海の流動と栄養塩に関するレビューと流動・低次生産モデルの開発
中山

哲嚴

＊

Review of characteristics on the current and nutrients distribution of
Seto Inland Sea, and development of the hydrodynamic model coupledwith primary production model applied to this region
Akiyoshi Nakayama＊

Abstract：Decrease in fishery resources in the semi-closed sea, such as Seto Inland Sea,
Ariake Bay and so on, has been a serious problem. The project started in 2006 for finding
how to improve or preserve their environment for restoring fishery resources in the inland
sea and how to keep or improve physical and chemical conditions closely related to fisheryresources and ecosystem. Reviews or analysis for investigating material circulation throughfood chain or ecosystem in the Seto Inland Sea were done in this project.
In this report, the current and nutrients transportation in the Seto Inland Sea as one ofimportant elements that directly influence the biological production are discussed by reviewing the results of studies on the large semi-closed sea such as the Seto Inland Sea. The resolved problems and unresolved problems are made clear. The current and nutrients transportation in the Seto Inland Sea are also discussed from point of the process of biologicalproduction, environmental capacity, resource management and increasing fishery resources.
According to the review, a hydrodynamic model coupled with primary production model
was developed and applied to the Seto Inland Sea. The results calculated by this modelshow a good agreement on the hydraulic characteristics of the Seto Inland Sea. But primary
production model included in this model should be improved by reevaluating its parameters
and its structure.
Key words；Inland sea, current, nutrients transportation, hydraulic model, primary production model
瀬戸内海や有明海などは水産資源の減少が深刻となっている。我が国有数の内湾域を豊饒の海へ回帰さ性を把握するための研究レビューやデータの解析など必要な作業を分担して行った。

せるために必要な環境とは何かを解明することを目標

本論文は，プロジェクトの課題の一つとして海域のに，漁業や生物環境に関連した栄養塩をはじめとする生物生産に直接的に影響を与える流動と栄養塩を取り

物理・化学環境などの望ましい水準を明らかにするた上げ，広域な閉鎖性水域に関するこれまでの研究成果

めの研究方向を具体化する目的で，地域連携プロジェをもとに，生物生産過程や環境容量の把握，資源管理・クト研究により，検討を平成18年に実施した。本プロ

増殖といった観点から検討したものである。さらに，

ジェクトでは，瀬戸内海を主な研究対象領域として，基本となる流動・低次生産モデルを構築し，瀬戸内海

食物連鎖を通した栄養段階間の物質循環の定量的関連全域に適用した結果を紹介する。

2010年 月30日受理（Received on August 30, 2010）
独立行政法人水産総合研究センター水産工学研究所水産土木工学部水産基盤グループ 〒314‑0408 茨城県神栖市波崎7620‑7（Leader of Fisheries Infrastructure Group，Aquaculture and Fishing Port Engineering Division，National Research Institute of Fisheries
Engineering，Fisheries Research Agency，Hasaki 7620-7，Kamisu-shi，Kashima-gun，Ibaragi-ken，Japan，zip code 314-0408）

Akiyoshi Nakayama
試料と方法
瀬戸内海という広大な閉鎖性海域に着目し，以下の点に主眼をおいて検討した。
⑴

結果及び考察流動（レビュー）瀬戸内海の地形を Fig. 1に示す。この図から，西側

瀬戸内海及び周辺外洋域を含めた流動に関する研究成果の収集と分析に豊後水道，東側に紀伊水道があり，この間を結ぶ
長大な水路とみることができる。この水路は単調で同海域における流動特性に関する研究成果を収集はなく，海峡・湾・灘によって構成されている。さらし，同海域での流動を決定づける潮汐，気象，河川流に陸域から数多くの河川が流れ込んでいる（

入に伴う全体的な流動の特性（密度特性，潮汐流・吹21）。瀬戸内海の流れは，豊後・紀伊両水道における送流・河川流等の相互干渉，
季節変動等）を把握する。

水位，流速，密度プロファイルと瀬戸内海に流れ込む⑵

級河川

瀬戸内海での長期的な栄養塩動態の実態把握及淡水の流入，風・気圧分布によって支配されており，

びデータの収集，河川・外洋からの栄養塩動態のこれらが相互に干渉して，流れが形成されると考えら

把握

れる。

同海域における栄養塩動態を長期的な視点で把握す

最初に既存の研究を検討し，瀬戸内海の流れを再るために，研究論文の収集，浅海定線調査・国交省に現・評価するために必要な重要課題を以下に列挙す

よる定点調査等多くのデータベースの収集，主要河川

る。

の流量・水質に関するデータベースの収集に努める。

⑴⑶

瀬戸内海の地形特性を適切に取り扱うことがで

豊後・紀伊両水道での水位，流れ，密度プロフ
ァイルの時空間推定（黒潮と豊後・紀伊両水道周き，流動や水質変動を比較的良く表現出来る数値
モデルの検討辺水隗の応答）⑵

多くの海峡や干潟を含んだ同海域において，⑴，⑵⑴と連動した両水道周辺の陸棚斜面水（低温・高塩分）の挙動（黒潮と豊後・紀伊両水道周辺水で抽出した課題に答え得る複数の数値モデルに関する

隗の応答）

研究成果及びソフトウエアの機能を検討し，適切なモ⑶

河川流入量の的確な把握（内海の密度流に関係）デルを選定し，瀬戸内海への適用のための改良処理を⑷

適切な拡散係数のモデル化又は設定（内海で行い，瀬戸内海全域の流動・水質等の概略検討を行う。

の混合や密度フロント周辺の流れ・物質輸送の再Fig. 1. Bathymetry of the Seto Inland Sea.

流動と栄養塩，低次生産モデル
現，強流域・緩流域の流れや成層状況を評価する
精度の良い計算が必要）⑸

では，外洋と接続する紀伊，豊後水道の流れはどうであろうか？紀伊水道，豊後水道に関する流動，水質

瀬戸内海での風外力の適切な設定（地形の影響を考慮した風速分布）に関する研究は，笠井ら（2001），高志ら（2002），藤原ら（2003），武岡（1996），兼田ら（2002），速水ら（2006）

各海域における流動構造に関する研究密度はまちま

によって系統的に行われてきた。紀伊水道では，黒潮

ちである。研究が数多く行われているのは，
豊後水道，の離接岸によって大きく流況が変化する。黒潮離岸時

紀伊水道，大阪湾，広島湾であろう。個々の海域でのに顕著に低水温水隗が底層から侵入し，その波及範囲

研究が多く，対象とする時期等による相違もある。こは，紀淡海峡を経て明石海峡にまで及ぶことがある。

こでは，瀬戸内海全体の傾向をつかみ，次にどのよう

豊後水道では急潮，「底入り潮」が直接的に宇和海をにして，流れの場を時空間的にとらえるか，あるいは

経て伊予灘南部海域に波及する。このような流れがこの流れを数値計算などで再現するに当たって検討す瀬戸内海の栄養塩分布に大きく影響すると考えられるべき課題を抽出することを主目的とした。従って，瀬（後述）。豊後水道では九州東岸への黒潮接岸時にそれ

戸内海全域を対象とした研究，外洋と内海との相互作

ぞれ低温・高栄養塩水（離棚斜面の冷水）の侵入が見用，海峡と湾での流れ構造，風による広域的な流れ等られる。どちらもキーとなるのは陸棚斜面に存在するに関する研究に絞って取り上げた。なお，流動と水質

低温・高栄養塩水隗の挙動である。低水温水隗の侵入は密接に関連しており，両者を区別することはできなに対応して，周辺水位もおそらく低下すると考えられい。本論では，流動と水質をあえて分けているが，実る。兼田（2002），藤原ら（2006）はこのことに着目際には同一の論文で検討されていることが多い。従っし，周辺水温と水位に高い相関があることを指摘してて，重複する部分もあることを記しておく。

いる。紀伊水道と豊後水道との違いは，紀伊水道に明まず，瀬戸内海全域における流れの特性について述

確な急潮がないことである。武岡（1996）は，このこべる。瀬戸内海は先に述べたように外洋と接続する紀とについて豊後水道の小潮時には水隗間に大きな密度

伊，豊後両水道に挟まれ，複数の湾・灘・海峡で構成差があり，内海に波及しやすいのに対して，紀伊水道

されている。橋本ら（1997）は，瀬戸内海を鉛直混合

では豊後水道に見られるような大きな密度差が無いた

強度の季節変動を指標としてめ，急潮となって内海へ波及しないとしている。

海域に区分し，一次生産量の地理的・季節別変動を調べて，海洋構造との関

Kobayashi and Fujiwara（2006）は陸棚内と外洋と係を考察した。船による観測データ（春，夏，秋，冬）の相互作用効果の把握及びモデル化を目的として，30を基に成層・混合状態を検討し，瀬戸内海をⅠ周年成年間の長期間にわたる数値計算による検討を行ってい層域，Ⅱ成層卓越域，Ⅲ混合卓越域，Ⅳ周年混合域のる（対象領域；備讃瀬戸より紀伊水道）。モデルは多つの海域に区分した。その結果，大阪湾奥，広島湾

層レベルモデルを用い，水平拡散係数，鉛直拡散係数，奥部をⅠ，播磨灘をⅡ，周防灘，伊予灘，広島湾，燧

大気との熱交換，河川流量データ，風データを用いて灘，備後灘，明石海峡を含む大阪湾西部をⅢ，備讃瀬

計算を実施した。計算結果は概ね観測結果を再現して

戸，安芸灘，伊予灘南部（豊予海峡）をⅣとした。こおり，水温・塩分，流動構造は，日射，河川流量，風，のように混合域，成層が卓越する海域が混在している

海水交換に依存しているとしている。また，境界条件のが瀬戸内海であるといえよう。混合域か，成層域か（外洋境界水温・塩分，河川流量），外力条件（風速）は，流動の大きさ，密度成層の度合いで変化し，各海に関して，年変動を考慮するか，しないかによって，域に流れ込む河川流量とも関係している。

どのような変動をもたらすかを検討するとともに，紀瀬戸内海を一つの水路とみると，全体的に東西どち

伊水道における水温と NO3 に高い相関があることをらに流れているのであろうか？藤原ら（2006）は，瀬用いて，外洋からの NO3 供給を計算した。その結果，

戸内海を通過する流量の長期変動に関して，通過流は

黒潮接岸時と離岸時では NO3 の分布が大きく異なっ

主に風，東西の水位差により生じるとして，経年変化ており，計算値・観測値から，NO3 に関する外洋起源を検討している。風による通過流は主に冬の季節風での割合を推定し，離岸時においては，NO3 は大阪湾西起こされるために年平均すると東向きとなる。水位差部及び播磨灘南部にまで達すること，接岸時では紀伊による通過流は長期的かつ大きな変動をする。西向き

水道南部までしか波及しないことを明らかにした。さ傾向の時期は，黒潮が大蛇行した時期に一致するとしらに，季節的な流れ・密度・NO3 の分布構造を示し，ている。このような通過流の長期変動は瀬戸内海内部

外洋，陸棚との相互作用は陸棚，内海域に大きな影響の塩分分布に大きな影響を与える。
通過流量の変動は，を与えること，内湾の物理・生物特性は黒潮など外洋

外洋域における黒潮の変動が大きく影響している。

での長期的な流動の変化に影響を受けていることを示Akiyoshi Nakayama
した。一連の検討結果から，瀬戸内海の流動や栄養塩による影響により水位の季節変動量が概ね説明できるの適切な推定手法として，黒潮離岸距離と紀伊水道のとしている。さらに駒井ら（2003）は，瀬戸内海総合

水温の相関関係，紀伊水道における水温・NO3 の相関

水質調査データを用いて，最近の20年間を対象とし，関係を利用することを提案している。

季節的に変動する瀬戸内海での流れ場について検討し田中ら（2005）は，多層モデルを用いて，瀬戸内海全域において潮汐（M潮成分）による流動を計算した。また，計算においては，風（西風m/s）を考慮した計算を行っている。計算では，単層，層の計算た。前述した計算モデルを改良し，水位・水温（黒潮流路），河川流出量の特性から分けられた

つの期間（1982‑88豊後水道水位高，89‑93紀伊水道高，河川流量大，94‑99紀伊水道水位大）に適用し，各期間の平を行い（ただし，密度一定），双方とも観測結果との均的な一致度は良かったが，位相誤差は年間の内海での流れ場（月毎）の比較を行っ

層レベルモデルのている。豊後水道から紀伊水道に向かう流れが卓越し方が少なかったとしている。瀬戸内海での数地点でのているが，内海への河川流出量や密度分布によって流断面流量を計算した結果，多層レベルモデルの方が，れ方向が変化することが明らかにされたとしている。

流量が若干小さくなることが分かった。季節風の影響

陸田ら（2003）は，複雑地形を有する瀬戸内圏を対が単層モデルでは小さくでる可能性があり，多層レベ象とした大気海洋モデルを開発し，16方位の風向によルモデルでの検討が必要であるとしている。

って形成される局地風の風況特性，吹送流の応答特

浜走ら（2001）は国交省のデータ，各地の水位デー性を検討した。瀬戸内海を対象とした大気海洋モデル

タ等を基本として，ボックスモデルを適用し，海域のは，局地気象モデル ARPS（The Advanced Regional

水塊体積時間変化率は，海域への流入，流出及び海域

Prediction System） と 海 洋 モ デ ル POM（Princeton

への河川流入量で決まるとして，海域間の海水移流量

Ocean Model，ただし密度一定）に基づきモデル構築を求めた。海域区分は伊予灘，広島湾・安芸灘，燧灘・がなされている。その結果，瀬戸内海では，海域周辺

備讃瀬戸，播磨灘・大阪湾・紀伊水道としている。移の流量の月ごとの変化から，全体として〜10月に移流が大きく異なることが明らかとなった。つまり，瀬戸

月には瀬戸内海西部，中央部で内海のように四方陸地に囲まれた閉鎖性海域の気象・量が多くなり，〜は外から内へ向かう流れが生じ，〜

月には瀬戸内次元的な陸地形の効果によって吹送流の応答特性

海象特性を把握するには，

次元局地風や瀬戸内海特海全体で内から外へ向かう流れが生じている。年平均

有の海岸・海底地形を無視することは出来ないとしての移流からは瀬戸内海から外海への流れとなり，特にいる。地形的に東西方向に開けた瀬戸内海は，その方

東部海域は年間を通して外海に向かう流れとなる。さ向の風によって大きな吹送流が発生しやすい。冬季のらに瀬戸内海東部を分割し，播磨灘，大阪湾，紀伊水

季節風により，東方向の流れが強化・卓越する傾向が道間の移流量を推定した結果から，あるといえる。

月頃から大阪湾

から紀伊水道への流れが生じ，同時に播磨灘から大阪

阿保ら（2006）は，冬季の東部瀬戸内海における無月頃から播磨灘から紀伊

機態窒素濃度に及ぼす残差流の影響を検討している。

水道に向かう流れが卓越すると播磨灘から大阪湾への密度場及び風速場から，診断モデルにより東部瀬戸内

流れは小さくなるとしている。塩分を保存量として海海の残差流を計算した。残差流は風の影響を強く受け，水交換量を求め，紀伊水道・外海間，豊後水道・伊予

季節風の強い灘間，備讃瀬戸・播磨灘間は他海域の数倍となること，

播磨灘北部では強い東向きの残差流が発生した。このまた，備讃瀬戸・播磨灘間では月に他海域の10残差流に伴い，栄養塩濃度の低い水塊が燧灘から東進倍近くになることを示した。東部海域での海水交換量

し，備讃瀬戸及び播磨灘の栄養塩が低下した。播磨灘に関しては，夏期に海水交換量は卓越する。播磨灘・北部の栄養塩濃度は，夏期に低く，秋に上昇して12月

大阪湾の海水交換は両海域の塩分逆転に対応してピーに最大となり，冬季になると低下した。この冬季の濃ク値を取っている。東部海域の海水交換量は西部海域

度低下は東向きの残差流発生後に起こり，西よりの季の伊予灘・広島湾間と同程度であるなどとしている。

節風が早期から強い年には栄養塩濃度の低下が早く，駒井ら（2001；2002）は，瀬戸内海における平年のノリ養殖不作の原因になることが指摘されている。

湾に向かう流れが生じる。

〜〜

月には東向きの通過流が発生し，密度，海面気圧と水位の季節変動量の関係について

海峡部の流れについては，小林ら（2006）は，瀬戸考察を行い，外海水位，密度，海面気圧による水位

内海中央部（備讃瀬戸と播磨灘）の潮汐フロント域に変動を推定した。外海水位を境界条件とした平面

次

おいて密度流の実測（流速・塩分プロファイル，鉛直

元数値解析モデルによって瀬戸内海の水位の季節変動

拡散係数）を行い，数値モデルと併せてフロントを横量を推定した結果，外海水位と密度分布及び海面気圧

切る流れが生じる機構について検討を行った。潮汐フ流動と栄養塩，低次生産モデル
ロント域では水平密度勾配による密度流が生じ，流れのようになっているのであろうか？海域での窒素，リの主成分はコリオリの効果によってフロントに従う方

ンの経年変化については，せとうちネット（2009）に向となる。しかし，それが陸岸に達した後は岸に沿っよると経年的な変化は少ないが，TN，TP はともにて流れるため，フロントを横切る方向となり，これに

近年では横ばい傾向で，COD，透明度は漸増傾向で対応した物質輸送が生じるとしている。これは陸岸にあった。平成15年度の海域別 COD は，大阪湾，備讃囲まれた潮汐フロント域に共通するメカニズムである

瀬戸，周防灘，播磨灘，燧灘，安芸灘，紀伊水道の順と考えられる。

であった。前述した負荷流入量と海域での濃度の順位が異なっているが，これは主にその面積及び流動環境

水質（レビュー）が海域によって異なっていることによるものと考えら

瀬戸内海の水質に関する既存の研究を検討して得られた主な結論，課題等をまとめると以下の様になる。
⑴

外洋からの栄養塩供給は無視出来ない。水温と栄養塩の相関は高い。黒潮の離・接岸が重要であ
る。

⑵

⑶

れるが，第

位は大阪湾で共通であり，負荷量も多く

かつ富栄養化が進んでいる状態となっていることが分かる。
瀬戸内海全体での栄養塩輸送機構については，藤原ら（2003）が種々のデータベースや現地調査等から検陸域からの栄養塩負荷は減少した。しかし，討している（ただし夏の場合）。外海水が紀伊・豊後

海域の栄養塩分布は必ずしもこれに対応していな

水道下層から瀬戸内海中央部に進み，備讃瀬戸で上昇い。

して，東西に分離し，上層へ数 km/ 日の流速で広が西部海域と東部海域では N/P 比が異なる。植る。海域毎にみると潮流の早い海峡部ごとに強混合域物プランクトンにとっては，西部海域では窒素制が存在する。紀伊水道側では底層の栄養塩が紀淡海峡

限，東部海域ではリン制限の傾向があると考えら

で鉛直混合を受け，一部は有光層，残りは大阪湾西部

れる。

底層を通過し明石海峡で有光層に供給される。豊後水⑷

大阪湾は，陸水負荷の寄与率が高い。

道側では，豊予海峡において有光層に供給され，伊予⑸

海峡部は強い流れのために栄養塩の水平輸送能

灘10m 層を通過し，来島海峡から再度底層を通過する。

力は高い。
⑹

この輸送ルートでは，伊予灘の有光層内で消費され，栄養塩輸送は，備讃瀬戸周辺を境界として，東西に分かれる傾向にある。

⑺

河川水に含まれる栄養塩濃度，懸濁態等の適切な推定が必要である。

来島海峡から燧灘下層への流入水は貧栄養となる。このために燧灘は全層で貧栄養となる。
柳を共通とする研究者らは，瀬戸内海での栄養塩動態（窒素・リンの現存比率，溶存態 N，P，Si の分布，

以下に，取り上げた既存の研究について述べる。

全リン・全窒素の起源と濃度変動機構，沿岸域・全域

瀬戸内海に流入するにおける TP・TN 濃度変動機構）について，一連の級河川は20程度と多い。これ

らの河川からの淡水流入による瀬戸内海への栄養塩検討を行っており，以下に述べる。

供給は大きい。年間総流量の変動が大きいが，長期な瀬戸内海で長年実施されている水質調査データ（全トレンドは見られない。李ら（1995）は，原単位法を

域に対する広域総合調査，環境省；沿岸に対する公共

用いて流入負荷量の瀬戸内海水域毎の COD，T‑N，用水域水質調査，地方公共団他）を基に，瀬戸内海各

T‑P の経年変化（昭和62年まで）を推定している。

灘・湾における水質制御指標を提案し，水域毎の特性

全体的な経年変化の傾向は，①昭和40年代の急激な増を検討した（石井・柳，2005）結果，富栄養化対策と加傾向，②それ以降の減少，③昭和50年前半から昭和

して実施される陸域からのリン・窒素負荷総量規制の60年代はほぼ横ばいの傾向となっている。最近の傾向

効果の現れやすい水域と現れにくい水域が存在するとは不明であるが，浄化施設の充実及び排出規制によるしている。陸域負荷の影響を受けやすい陸域から概ね汚濁負荷の減少等を考慮すると減少傾向であると考え〜

km の領域を「沿岸域」とし，紀伊水道から豊られる。各海域別に見ると，大阪湾への負荷流入が最後水道を含む（響灘を除く）水域全体を「全域」としも大きく，全体の37 1% を占めている。続いて播磨灘，て検討している。解析対象は鉛直混合が活発な冬季の周防灘，紀伊水道，燧灘，安芸灘，伊予灘，備讃瀬戸

データを用い，1979〜1998年の期間では，沿岸域ではの順である。そのほか，海域別に負荷量についての研TP，TN とも減少していること，TP の減少が大きい

究成果があるが，上記の経年変化の傾向は一致していこと，全域では TP に関しては減少傾向にならず，むる。

しろ微増傾向であること，全域の TN も長期トレンド

では実際の海域別での栄養塩等水質の長期的変はどとしては微増傾向であることを明らかにした。陸域かAkiyoshi Nakayama
らの負荷が減少しているのに，全域の TP，TN が増倍であった。DSi は，別府湾，大阪湾で高く，燧灘で加している理由として，
外洋起源の栄養塩流入の増加，

わずかに高く，伊予灘で最も低かった。また，海域別

底質からのリン・窒素の溶出増加の可能性を挙げていに陸域からの負荷量を推定した結果，単位容積あたりる。DO の経年変化は微増傾向であることから，底質では TN，TP 負荷量は備讃瀬戸と大阪湾で大きい。

からの溶出増の可能性は少ないのではないかと考えら

DIP，DIN を海域別分布から見ると大阪湾は対応してれた。各海域での沿岸・湾灘全域での沿岸域濃度と全いるが，備讃瀬戸は対応していない。播磨灘は備讃

域濃度の比 C-Ratio を検討した結果，紀伊水道が最も瀬戸より負荷量は小さいのに DIP，DIN は逆に大き大きく，大阪湾，備讃瀬戸の順であったが，紀伊水道

い。この理由としては，備讃瀬戸の強い潮流により，が最大という点には疑問が残るとしている。C-Ratio

負荷された栄養塩は速やかに播磨灘，燧灘に輸送されを水質制御指標として，各水域での陸域負荷量削減対ると同時に，大阪湾から播磨灘への栄養塩が輸送され策による沿岸域での水質改善効果の現れやすさについてくるためとしている。DSi 負荷量に関しては，備讃て検討した結果，陸域での負荷削減対策が効果を上げ

瀬戸，別府湾，大阪湾の順となっている。しかし，備るのは紀伊水道，大阪湾であり，播磨灘，安芸灘，広讃瀬戸では DSi 濃度は低く，潮流による水平輸送が島湾，別府湾は効果が現れにくいとしている。

大きいためであるとしている。DIN/DIP，DSi/DIP，

橋本らによる1993年（10月），1994年（，，

DSi/DIN の海域別分布から，大阪湾は，著しく N/P月）の観測データから，形態別窒素・リン・珪素の海が高く Si/N が低い。これは淀川の人為的負荷が大き

域別・季節別変化についての検討が行われた（林ら，いこと，琵琶湖での植物プランクトン増殖による珪素

2000）。瀬戸内海全体の年平均では DOP，DIP，POP

取り込みによる結果としている。別府湾は河川流量がの現存割合は 1/3 程度であり，窒素も同様であった。

相対的に多いので Si/P が大きい。各水域間での栄養

生物活動が活発なために全体的に DIN の割合が低く

塩フラックスの推定を行った結果，DIP は大阪湾からなり，DON が高くなる。特に

月燧灘〜伊予灘・周別府湾までの水平輸送が存在する。DIN に関しては，

防灘にかけて DIN が非常に低い。10月は DIN の比率

大阪湾から燧灘への水平輸送，安芸灘・別府湾から伊が高くなり，DON が低くなる。海域別では広島湾北予灘への水平輸送がある。DSi は大阪湾から備讃瀬戸

部でプランクトン含有窒素が年間を通して高く，大経由播磨灘，燧灘から備讃瀬戸・安芸灘，安芸灘・別阪湾北部ではプランクトン含有率窒素は先の海域に府湾から伊予灘への水平輸送があるとしている。

比べて低い。この理由としては，同海域の基礎生産

環境省により実施されている瀬戸内海総合水質調査

は瀬戸内海中最も高いが，陸水による栄養塩供給がで取得された TP，TN 濃度から1981〜2000年の冬季

過剰なためとしている。リンに関しては，総じて窒素

水域別平均値を使用して，瀬戸内海における陸域と外同様の傾向となっており，

洋起源の TP，TN 濃度推定を試みている（石井・柳，・

月の DIP 現存割合

が DIN に比べて大きいこと，DOP の割合が DON に2004）。大阪湾では陸起源の TP が比べて低いことが異なる点である。この傾向は瀬戸内も存在し，陸域からの影響が支配的であることが分か

割，TN が割

海西部で顕著である。これらのことから，瀬戸内海西った。これは同湾への TP，TN 負荷量が瀬戸内海総部海域での基礎生産が窒素制限の傾向にあることを示

負荷量の 1/3 と非常に多い上に，内海最大の流量の淀しているとしている。Si/DIP，DIN/DIP の海域別分

川を経由して湾内に流入するためであるとしている。

布から，大阪湾での Si/DIP は年間総じて低く，ケイ

一方，瀬戸内海西部では外洋からの寄与率が藻の生産が制限される可能性の高い海域としている。

を占めており，東西で対照的な水域特性を有すること

DIN/DIP では，西部海域で， 月ともに低い傾向にが判明した。各水域における外洋起源 TP，TN 濃度あり，窒素制限傾向であり，備讃瀬戸以東海域では変動が紀伊水道・豊後水道のどちらを通じて影響を

DIN/DIP がレッドフィールド比を上回っており，リ強く受けているかについて検討した結果，平均的には

ン制限傾向であるとしている。

TP，TN ともに備讃瀬戸付近を境に西では豊後水道，

フェリーによる DIP，DIP，DSi 観測データ（1994〜2000年）から，鉛直混合が盛んな冬季の観測結果を利用し，海域別の栄養塩特性を分析した（柳・原島，
2003）
。

割以上

東では紀伊水道を通じての影響が有意であるとしている。以上までが柳らの一連の研究である。
瀬戸内海での外洋起源の栄養塩比率に関しては，武年間平均の冬季栄養塩分布から，DIP は，岡，柳（前述）をはじめとしたいくつかの研究がある。

瀬戸内海東部で高く，西部で低かった。大阪湾は他

これらを比較するとそれぞれの研究でこの比は異なっ

海域のているが，外洋起源窒素比率は 2/3 程度であるとして〜倍であった。DIN は，別府湾を除くと

DIP と同様な傾向にあった。大阪湾は他海域の〜いる（武岡ら，2005）。解析においては，河川からの流動と栄養塩，低次生産モデル
栄養塩供給なしの場合での瀬戸内海の栄養塩濃度は豊

まかに紀伊水道と共通しており，1980年代後半は低水

後・紀伊両水道に接する栄養塩濃度と等しくなるとい

温・高栄養塩期となっている。また，大阪湾での栄養う考え方から算出した栄養塩濃度を Background 濃度

塩・水温変動を分析して，外海から流入する栄養塩量とした。この見積もりは Background 濃度をどの程度の30年間の経年変化を調べた。大阪湾西部底層には，にするかなど問題点があるものの，外洋起源が少なか

紀伊水道から外海起源の栄養塩が流入しており，このらず影響していることを示している。

流入が大阪湾内の水質に影響を与えている。また，外外洋域に接続している豊後水道，紀伊水道からの栄

海から大阪湾への栄養塩の流入量には大きな経年変化

養塩供給については，
笠井ら（2001），
高志ら（2002），があり，1980年代後半に極大があり，その後減少し，藤原ら（2003），武岡（1996），兼田ら（2002），速水

1990年代後半には極大時の約ら（2006）が系統的な研究を継続的に行い，両水道から，この変動が瀬戸内海の水質の長期変動に影響を及らの栄養塩供給は無視出来ない量であることが指摘さぼしているとしている。外洋からの栄養塩供給の多寡れている。一連の研究を総括してみると以下のようなは，陸棚斜面水の動きにより決定されており，黒潮のことが言える。

動きが大きく影響していることを示している。

紀伊水道においては，栄養塩の流出入は黒潮の離接分のになったことか

河川水に含まれる栄養物質は，無機溶存態，有機溶

岸に大きく影響を受けており，黒潮接岸期においては存態，懸濁態（POM）等である。これらは源流，上主として栄養塩の流出が主であり，黒潮が離岸してい流域，中流域，下流域で変化しているが，沿岸域ではると底層からの栄養塩の流入が卓越する。ただし，こ河口域でのこれらの量が重要になる。また，河口域はの現象は夏期に限られ，冬季は黒潮の離接岸に係わら海水と河川水が混合する領域である。河川水に含まれず，流出傾向にあるとしている。豊後水道では，以前

る栄養塩は，
流域の土地利用，
点源負荷
（工場，
畜産場，

から急潮と呼ばれる黒潮暖水の表層侵入とこのイベン

汚水処理場等），地質構造等で変化する。河川水に含ト後に底層の栄養塩が増加することが知られている。

まれる栄養塩濃度は季節的にも変化する。例えば，新急潮後に湾外からの低温水（陸棚斜面起源）侵入がそ井ら（2006）は利根川の栄養塩を年間測定したが，その原因と考えられている。加えて，武岡らが「底入りの結果，秋期・冬季には硝酸態窒素は200μmol/l と，潮」という底層からの低温・高塩分水塊が宇和海にま

春期・夏期に比べてで侵入しており，この二つの現象が外洋からの栄養塩

いる。さらに出水時における硝酸態窒素濃度を連続計

供給メカニズムであるとしている。さらにこれらの現測し，出水に伴う濃度変動はそれ程大きくないことを

象は周辺海域での動植物プランクトン組成にも大きく報告している。また，
懸濁態に含まれる窒素・リンは，影響を与えているとしている。両水道での高栄養塩水

SS 濃度と比較的相関が高いことも報告している。山隗の侵入は，いずれも黒潮の離接岸に大きく関係して本ら（2003）は涸沼川水系（農地河川），白川（阿蘇いると見られる。ただし，紀伊水道では黒潮離岸時，カルデラ集水域）について，観測を行い，溶存態・懸豊後水道では九州東岸への黒潮接岸時にそれぞれ高栄

濁態，濁度（SS 濃度），流量の関係を検討している。

養塩・低温水（陸棚斜面の冷水）の侵入が見られる。

その結果，同様に SS 濃度と窒素・リンは相関が高く，

これは，それぞれの水道周辺の地形条件，陸棚斜面起SS は流量と関係があるが，一般にその関係はヒステ

源の冷水領域，黒潮の位置関係及び黒潮の動きによっ

リシスカーブとなる。一方，溶存態窒素 NO3 はヒスて生起条件が異なってくるためと考えられる。その直テリシスではなく流量とほぼ一次の線形関係にある場

接的影響範囲は，紀伊水道側では大阪湾西部海域，豊合が多いとしている。それぞれの河川は，土地利用や後水道では伊予灘南部海域に及んでいるとしている。

点源負荷等に応じて異なっており，この点についてさ従って，間接的にはそれよりかなり広い範囲に影響をらに詳細な検討が必要になる。しかし，河川データは，与えるものと考えられる。外海から流入する栄養塩濃TN，TP，SS データに限られていることが多く，溶度は，紀伊・豊後水道ともに水温との相関が高い。

存態に関するデータは少なく，溶存態栄養塩をこれら

藤原ら（2006）は，豊後水道の水深100m 水温と細倍の濃度となることを報告してのデータから推定せざるを得ない状況にある。

島の水位がほぼ連動していることから，細島の水位の平水期あるいは渇水期には，流速が遅いうえ，ダ

長期データから100m 水温の長期変動を推定した。水ム，河口堰など流域に滞留域がある場合には，水温が温と栄養塩濃度は高い負の相関があるので，これを用植物プランクトンにとって適切であるならば増殖する

いて100m での栄養塩濃度長期変動を推定した。その結果，溶存態が減少し，懸濁態に移行するメカニズム

結果，数十年スケールの変動があり，硝酸態窒素濃度が働く。しかし，水温が低い状態では起こりにくいと

μM で変動している。変動のパターンは大考えられる。冬季に栄養塩が高くなるのは，流域でのは〜Akiyoshi Nakayama
負荷発生状況も絡んでくるが，河川水域においての一向，風速，気温，湿度，全天日射量，雲量等）を入力

次生産の活性が低くなることも大きな要因の一つではしやすいこと，④計算メッシュを小さくするために計ないかと考えられる。

算容量が大きくなることからロバスト性に優れている
こと，⑤生態系（栄養塩も含む）モデルを容易に組み

底質・ベントス等（レビュー）
水質に影響を与える重要な問題として，瀬戸内海に込めること，⑥ネスティングや可変格子形成が容易なことなどを考慮する必要がある。

おける底質分布がある。海域での生産過程や陸水から

海洋流動モデルとして，最も広く利用されているのの有機物負荷により，有機物は海底に堆積する。堆積が，プリンストン・オーシャンモデル（POM）である。

した有機物は，ベントス，細菌類等により，消費分解

POM は，以下に示すように，される。しかし，負荷量が多くなると未分解の有機物あり，拡散係数を水平方向についてはスマゴリンスキの堆積が増加し，これまで良好であった底質環境が悪ー型とし，鉛直方向については乱流輸送エネルギー方

化し，夏期の貧酸素化，底質からの栄養塩溶出増加に

程式（Mellor & Yamada，1982）から求めるようによる内部生産の増加によるさらなる貧酸素化・有機物なっている。計算の手順としては，まず，の堆積という悪循環サイクルに陥る。通常の環境であれを計算し（外部モード），それより長い時間間隔でればベントス類が生息し，適切な物質循環が行われる

傾圧的な流れを計算するようになっている（内部モーが，有機汚泥化した領域では，夏期の貧酸素化等によド）。風応力は海面摩擦係数を介して，水塊に作用する。

りベントス類の避難あるいは斃死が発生する結果，ベ

POM の計算に関する記述は付録に示す。

ントス群集のサクセションをストップ（またはリセット）させてしまう。

次元の流動モデルで次元の流POM の問題点としては，σ座標系でのスカラー量
輸送計算における数値分散の介入がある。直交座標

せとうちネット（2009）に底質の含泥率，COD，系における密度の定義位置とσ座標系における密度

マクロベントス種類・個体数の分布が公表されていの定義位置が異なることによる計算上の誤差である。

る。含泥率分布を見ると，瀬戸内海の流れを反映し，POM を ベ ー ス に 様 々 な 改 良 が 加 え ら れ，ROMS，

流れの速い海域では，含泥率が低く，閉鎖度が高い海

ECOMSED，EFDC などが開発された。その過程の域ほど高い分布となっている。COD の分布もほぼ同中で問題点である数値分散等が改善された。POM 自様である。マクロベントスの個体数・種類数に関して

身は様々なユーザーにより改良されるとともに，生態は，含泥率・COD の分布とは異なり，閉鎖性が強い系モデル，
気象モデルと結合したモデル等が開発され，海域の主要な領域では少なく，むしろ流れの強い海域

海洋分野や地球温暖化問題で利用されているが，これ周辺に多い傾向にある。個体数・種類数のみでなく，らすべてが公開されているわけではない。モデルの基その構成及びその経年変化が重要であると考えられる

本的特徴（拡散係数，座標系，干出・遡上モデルの有が，この分布を見るだけでも，海域の特性を反映した

無，土砂の移流拡散・堆積モデル有無等）の一覧表を

分布特性を示しており，各海域の生産性，負荷，水理示す（Table 1）。

環境が密接に関連していることが伺える。大阪湾や広瀬戸内海に適用するモデルは，先に述べた流動，

島湾での研究成果等から，底質からの栄養塩溶出は，水質，底質等に関する課題群に適切に対応出来るモ

季節的には夏期に大きく，
冬季に少ないと考えられる。

デルでなければならない。著者らは，以前から開放性

特に夏季における貧酸素化は，大阪湾，広島湾奥では沿岸域である鹿島灘・九十九里を対象に流動・低次生深刻である。こうした底質からの栄養塩供給は，瀬戸産モデルの開発を行ってきた。そのベースは POM で内海の低次生産や環境を評価する上では無視することある。これに低次生態系モデルを結合させている。

は出来ない。

同海域には黒潮，親潮続流が存在し，外洋の流れが直接沿岸域に波及する。このため，これら海流を適切に数値モデル
数値モデル概要

考慮することが必要である。また，利根川，那珂川という級河川があり，これらの影響（淡水流入，懸濁既に述べたが，瀬戸内海は非常に複雑な地形をして物，栄養塩等）を考慮することも必要である。前者にいること，周防灘などのように浅く，奥には干潟が広

ついては POM と沿岸・外洋観測値を用いた同化モデがっている領域もある。河川流量も多い。このことか

ルにより求められる流動，水位，塩分・水温が，既にら，流動の数値モデルを適用する場合には，①地形効

FRAJCOPE，JCOPE データ（

果を良く反映出来ること，②淡水流入に伴う密度成層ており，このデータを利用することができるようにモ③外洋の条件や気象（風の形成を良く再現出来ること，デルを構築している。後者については，河川データを日平均）で公表され流動と栄養塩，低次生産モデル
Table 1. Characteristics of hydro dynamic models

Fig. 2. Calculation domain of the Seto Inland Sea.

Akiyoshi Nakayama
用いて淡水流入，栄養塩供給等もモデルに組み込まれTable 2. A table of conditions used in the Model.

ている。POM は，他のモデルよりコード数が少なく，比較的容易に改良できることが大きな利点である。先に挙げた諸課題に対応する基本的なモデルを構築する
ために，このモデルを用いることが最も効率的と考え，
瀬戸内海に適用出来るように改良を加えた。なお，本モデルは現在でも多重ネスティング，干潟での干出・冠水等種々の改良作業が継続されている。ここでは，途中の段階の計算結果等を紹介する。なお，低次生態系モデルは KKYS モデル（Kawamiya ら，1995）をベースとし，リンのコンパートメント及び底泥の栄養塩溶出項を組み込んでいる。
以下にモデルの概要，条件等を述べる。
瀬戸内海モデルの計算領域は，瀬戸内海全域を含む
南北234km，東西421km の領域とし，豊後水道，紀伊水道，関門海峡を開境界とした（Fig. 2）。
既に述べたが，瀬戸内海の流れを駆動する主要な外間して海域ごとに設定した。外洋の流速・水温・塩分

力は潮汐及び風である。潮流は往復流であるため，正は，JCOPE 同化データ（宮澤・山形，2003）を開境

味の物質輸送は潮汐残差流・吹送流・密度流によって界に与えた。JCOPE データの空間分解能は本計算の決定される（柳，1992）。物質量の供給源は河川（陸

空間分解能より大きいので，内挿して境界値として用起源），底泥からの溶出に加えて，前述したように外いた。外洋の硝酸塩（以下，NO3）及びリン酸塩（以下，

洋からの供給が重要であることが明らかにされてきPO4）は，日本海洋データセンター（以下，JODC）

た。瀬戸内海モデルを構築する際には，上記の影響をのデータを解析して得られた水温と NO3 及び NO3 と適切に取り込むことが重要である。

PO4 の関係式（Fig. 3）を使用して，開境界の水温か

Table 2に瀬戸内海モデルの概要を示す。主要な外ら NO3 ，その NO3 から PO4 を算出して開境界に与え力である潮汐は，紀伊水道と豊後水道の開境界に天文た。ただし，40m 以浅では水温と NO3 の間に明確な潮位を与え，風は気象庁の GPV（RSM）データを使

相関関係は見られなかったため，40m 以浅の平均値用して時空間的な変化を考慮した。河川は瀬戸内海にを与えた。

流入する21の級河川を考慮し，底泥からの栄養塩溶

出量は山本ら（1998）の結果を時間的にスプライン補Fig. 3. The relationships between NO3 and water temperature and PO3.

流動と栄養塩，低次生産モデル
モデルの計算結果⑴など多数の研究がなされており，大阪湾の時計回りの流動場計算結果の検証

沖ノ瀬環流や広島湾の反時計回りの循環流が知られて流れを駆動する主要な外力である潮位の再現性にいる。計算による

年間の表層平均流速から（Fig. 5），

ついては，気象庁による観測潮位と計算結果の水位をそれら既知の循環流等を概ね表現できていることがそれぞれ調和分解して得られた分潮成分（M2，S2，

確かめられた。藤原ら（2006）によると，瀬戸内海をK1，O1，K2，N2，P1）の振幅及び遅角を比較した

通過する流量は10年スケールの変動があるが，豊後水（Fig. 4）
。その結果，相関係数はそれぞれ0 94及び0 98

道から流入して紀伊水道から流出する年が多い。計算であり，潮位の再現性は良好であると考えられた。正による月ごとの平均流速と水深から通過流量を算出し味の物質輸送を支配する循環流，恒流は柳ら（1992）

た結果，2005年は豊後水道から流入して紀伊水道からFig. 4. Comparison of field data to calculated data (amplitude and phase of tide).

Fig. 5. Calculated residual current vector (surface).

Akiyoshi Nakayama
る。計算結果は，紀伊水道の月，

月では実測値に比べて水温は高めであるが，それを除けば，観測結果を適切に再現できているものと考えられる。塩分に関しては，観測結果では，外洋に面した伊予灘及び紀伊水道以外では，量的な差はあるものの時間の経過と共に高塩分化する傾向があり，燧灘以東で特にその傾向が強い（Fig. 8左からFig. 6. Water discharge rate of the Kii channel and21 rivers in the Seto Inland Sea coast.

列目）。計算結果の塩分は，そのような高塩分化の傾向を概ね再現しているものの，全体的に観測結果より高めであった。瀬戸内海に流入する級河川の流域面積の合計は級河川の約 1/2 に相当する（せとうちネット，2009）。それゆえ，級流出する傾向がみられ，その流出量は瀬戸内海に流入河川のみを考慮した本モデルでは，流入する淡水量がする21の実際よりも少なかった可能性がある。更に，塩分の初級河川の総流量の10倍程度であった（Fig.

6）
。

期値を瀬戸内海全域の観測値が入手できた瀬戸内総合以上より，本モデルは瀬戸内海の流動場を概ね再現できているものと考えられる。
⑵の2005年冬季データ（

月〜

計算初期値が実際の月の塩分より幾分高めであった

月）から作成したため，

水温，塩分及び低次生産計算結果の検証可能性もあり，今後改善する必要がある。また，図に

水 温， 塩 分，NO3 及 び ク ロ ロ フ ィ ル a（ 以 下，は示していないが，塩分・水温分布による成層強度を

Chl.a）について，瀬戸内海総合水質調査データ（以見ると，計算結果は観測結果よりも弱く，モデルにお下，瀬戸内総合）及び各県の浅海定線データ（以下，

いて混合が過大である可能性がある。これは栄養塩輸

浅海定線）の各海域平均値と計算結果を比較検証した。

送，低次生産に直接影響するので，改善する必要があFig. 7に両データの観測地点を示す。水温に関しては，る。改善方策としては，拡散係数の再検討や鉛直分割

観測結果では，2005年の瀬戸内海における水温は数の増加などが考えられる。

頃に10℃程度の極小値をとり，に伴い水温成層が発達して

月

月以降に水温の上昇

なお，2005年当初の塩分は平年に比較して低い傾向

月に表層と底層の水温差にある。これは前年の2004年に台風等による出水が多が最も大きくなる（Fig. 8最左列）。表層水温は月に約25℃の極大値となり，底層水温は10月に極大値となかったためと考えられる。
観測結果では，NO3 はFig. 7. Location and number of observation points.

月から減少傾向で月，

流動と栄養塩，低次生産モデル

Fig. 8. Field data and calculated data.

月からから冬季に濃度が上昇する傾向を再現しているが，春月は海域ごとに傾季における NO3 の枯渇状況は再現できていない。こ向が異なり，
10月以降は増加傾向であった（Fig. 8右かれは，春季の Chl.a 極大を過小評価したためと考えら列目）。一方，Chl.a は NO3 が大きく減少する

れる。春季以降は，Chl.a の全体的な濃度は概ね再現

月には1 0mmol/L 以下まで減少している。
月は低濃度で推移し，
ら

月又は

月から

月に極大となり（春季ブルーム），次の月にできているが，一時的に高濃度となる状況は再現できは減少して安芸灘や広島湾では 1/2 以下になっているていない。春季ブルームを再現できていない原因とし（Fig. 8最右列）。 月以降は， 月， 月に表層， 月，ては低次生産モデルのパラメータ設定及び計算開始か10月には表層と底層で一時的に Chl.a が増加しているら海域がある。計算結果では，夏季に NO3 が低く秋季

しえる濃度ではなかったことが考えられる。今後は，ヶ月後であるため，各物質が春季ブルームを再現Akiyoshi Nakayama
参考文献

複数年にわたる検証などを通じて春季ブルームや一時的な Chl.a 濃度の上昇を支配する要因を解明し，モデ
ルの精度を向上させる必要がある。
なお，図には示していないが，紀伊，豊後両水道の底層から栄養塩が瀬戸内海に輸送されており，その強さは外洋の条件（境界条件）により変動していること

阿保勝之，杜多哲，藤原建紀，2006：冬季の東部瀬戸内海における残差流と栄養塩環境．海岸工学論
文集，53，1096‑1100．
新 井 雅 之， 中 山 哲 嚴， 足 立 久 美 子， 齊 藤肇， 奥が計算結果に現れている。詳細な検討は今後行う予定

西であるが，前述した外洋からの栄養塩供給機構を反映

岸での一次生産に及ぼす利根川・那珂川の影響に

武，八木 宏，2006：鹿島灘・九十九里浜沿していると考えられる。

ついて．海岸工学論文集，53，1101‑1105．
藤原建紀，小林志保，高志利宣，2003：瀬戸内海の窒まとめ素・リンの輸送と起源の現地観測．海岸工学論文
集，50，951‑955．

以上より，これまでの調査研究成果レビューを行藤原建紀，小林志保，國井麻妃，宇野奈津子，2006：

い，諸課題の解決に向けて必要な基本的な流動・低次

瀬戸内海に存在するリン・窒素量の経年変動．沿生態モデルを，流動モデルとして POM，低次生態系

岸海洋研究ノート，43⑵，129‑136．

モデルとして KKYS モデルをベースに開発し，瀬戸

浜走幸育，日比野忠史，駒井克昭，松岡純作，2001：

内海に適用した。低次生態系モデルの諸パラメータは瀬戸内海における湾・灘水塊の移流・交換量の季

基本的に KKYS のデフォルト値を用いており，再現性に課題を多く残している。瀬戸内海における基礎生
産調査データ等による生態系モデルのパラメータの見直し，生態系モデルの構造，底泥からの溶出過程等を詳細に検討する必要がある。さらには流動モデルにお節変動特性．海岸工学論文集，48，1036‑1040．
橋本俊也，山本民次，多田邦尚，松田治，永末寿宏，1997：瀬戸内海の一次生産と海洋構造．沿岸海洋研究ノート，35⑴，109‑114．
速水祐一，兼田淳史，小濱剛，中野伸一，武岡英隆，

いても，既に述べたように混合が過大気味のために成2006：豊後水道における外洋起源栄養塩の供給機

層強度が弱めの計算結果を与えているなどの問題があ構とその生態系への影響．沿岸海洋研究ノート，る。これらの点に関しては，現在，鋭意改良作業を継43⑵，143‑149．

続している。また，干潟上の流れなどを精度良く再現

石井大輔，柳哲雄，2004：瀬戸内海各灘・湾におけするために一部領域をより解像度を上げる多重ネステる全リン・全窒素の起源と濃度変動機構．海の研

ィングや干出・冠水を表現出来る移動境界などの検討

究，13，389‑401．

も行っている。本論文では基本モデルの紹介にとどめ石井大輔，柳哲雄，2005：瀬戸内海の沿岸域と全ているが，上記の改良・開発に関しては成果が出次第

域における TP・TN 濃度変動．海の研究，14，報告する予定である。

35‑45．

今後，本モデルとデータベース，現地観測等を総合

石井大輔，柳哲雄，2006：瀬戸内海に存在する太平的に利用することにより，水産有用種の産卵場から幼洋起源のリン・窒素．沿岸海洋研究ノート，
43⑵，

生分散・着底過程の検討，着底場における生産環境評

119‑127．

価，漁場環境の改善効果評価，中長期の瀬戸内海環境

笠井亮秀，藤原建紀，多田光男，2001：紀伊水道の変動予測，準リアルタイム流動・低次生産情報把握な海洋構造と栄養塩輸送．海岸工学論文集，48，

どについて検討することを考えている。

436‑440．

なお，浅海定線データ等の提供等に関しては，瀬戸兼田淳史，乗松桂輔，渡辺浩三，小泉喜嗣，武岡英隆，

内海沿岸の各県の水産試験場並びに県職員の方々に協

2002：黒潮の離接岸が豊後水道の水温に与える影力頂いた。

級河川の流量，水質データの提供は該当

響．沿岸海洋研究ノート，39⑵，181‑188．

する21の級河川の国交省下流事務所の方々に協力頂

Kawamiya, M., M. Kishi, Y. Yamanaka, and N.

いた。それぞれ関係機関，氏名を明記すべきところでSuginohara, 1995：An Ecological-Physical

すが，あまりに多いので，省略させていただきました。

Coupled Model Applied to Station Papa.

ここに謝意を表します。

51, 635‑664.
Kobayashi, K. and T. Fujiwara, 2006：Physical andbiochemical responses to decadal variability of
shelf water intrusion in the Seto Inland Sea,

流動と栄養塩，低次生産モデル
Japan. Workshop on Coastal Observatories

よびリン溶出量の見積もり．海の研究，- Best practice in the synthesis of long-term

151‑158．

observations and models, Proudman

山本浩一・二村貴幸・坂野

Oceanographic Laboratory (Liverpool, UK).

⑶，章・日下部隆昭・末次忠

司・横山勝英，2003：濁度計による懸濁態栄養塩

小林志保，藤原建紀，高杉由夫，2006：陸岸に囲まれた潮汐フロント域における物質輸送機構．海岸工

負荷推定に関する研究．河川技術論文集，
柳哲雄，1992：瀬戸内海の流動とそれに関する諸問

学論文集，53，921‑925．

題．海の研究，

駒井克昭，日比野忠史，清水勝義，2002：瀬戸内海に

，

柳

哲雄，原島

⑵，9‑19．

省，2003：瀬戸内海における溶存態

おける水位の季節変動量の推定．
海岸工学論文集，

無機リン・窒素・珪素分布の特徴とその要因．海49，381‑385．

の研究，12，565‑572．

駒 井 克 昭， 竹 内 健 太 郎， 日 比 野 忠 史， 松 本 英 雄，
2003：瀬戸内海における湾・灘間での海水交流量＊以下は，本文中で引用していないが，関係論文としの推定およびその長期変動に関する研究．50，
921‑925．

て掲載する。
Blumberg，A. F. and G. L. Mellor, 1987：A

駒井克昭，日比野忠史，仁子幸子，浜走幸育，松岡

description of a three- dimensional coastal ocean

純作，2001：瀬戸内海における栄養塩と一次生

circulation model, in Three-Dimensional Coastal

産力の分布・変動特性．海岸工学論文集，48，

ocean Models (ed. by N. Heaps), 1‑16.

1106‑1110．

呉碵津，松山幸彦，山本民次，中嶋昌紀，高辻英之，Mellor, G. L., and T. Yamada, 1982: Development

藤沢邦康，2005：近年の瀬戸内海における有害・

of a turbulence closure model for geophysical

有毒渦鞭毛藻の分布拡大とその原因−溶存態有機

20,

リンの生態学的重要性−．沿岸海洋研究ノート，

fluid problems.

43⑴，85‑95．

851‑875.
宮澤泰正，山形俊男，2003：JCOPE 海洋変動予測シ

ける窒素・リンの存在量とその長期変動．海と空，

ステム．月刊海洋，12，881‑886．
李弘，1995：瀬80⑵，75‑78．

戸内海への栄養塩類の流入負荷解析と物質収支に

日比野忠史，浅井寅鉄，浮田正夫，関根雅彦，中西関する研究．海岸工学論文集，42，1106‑1110．
陸田秀実，市位嘉崇，秋山佳明，土井康明，2003：局
地気象モデルを用いた瀬戸内圏の風況解析と吹送林

速水祐一，碓井澄子，武岡英隆，2004：瀬戸内海にお正，細川恭史，2000：瀬戸内海における湾・灘での水質の季節変動特性．海岸工学
論文集，47，1036‑1040．
眞鍋武彦，反田寛，堀豊，長井敏，中村行延，流の応答特性．海岸工学論文集，50，436‑440．

1994：播磨灘の漁場環境と植物プランクトンの変

美鶴，柳動─20年間のモニタリングの成果─．沿岸海洋研哲雄，橋本俊也，2000：瀬戸内海における窒素・リン現存量比率．海の研究，⑵，

せとうちネット，2009：瀬戸内海の環境情報─自然環
境に関する情報─．http://www.seto.or.jp/seto/
kankyojoho/sizenkankyo/top-data/hyou1-23.htm
参照2009/5/17.
高志利宣，藤原建紀，住友寿明，竹内淳一，2002：外洋から紀伊水道への窒素・リンの輸送．海岸工学
論文集，49，1076‑1080．
武岡英隆，1996：沿岸海洋と外洋の相互作用．沿岸海洋研究ノート，34⑴，3‑13．
田中陽二，河原能久，末永慶寛，2005：瀬戸内海全域
における海水流動特性の数値実験．水圏環境工学
論文集，

究ノート，31⑵，169‑181．
丸谷尊彦，韓83‑89．

⑴，9‑14.

山本民次，松田治，橋本俊也，妹背秀和，北村智顕，1998：瀬戸内海底泥からの溶存無機態窒素お

銅珍，中辻啓二，2000：大阪湾における水質の長期予測計算．海岸工学論文集，47，
1051‑1055．
中村充博，小林弘二朗，村上和男，2005：瀬戸内海の水質環境に及ぼす黒潮流路と外海の海況変動の影
響．海岸工学論文集，52，361 365．
西田修三，金漢九，高地慶，入江政安，中辻啓二，

2006：紀淡海峡における水質変動特性と栄養塩輸送．海岸工学論文集，53，996‑1000．
陸田秀実，土井康明，稲山丈晴，2002：瀬戸内海における栄養塩構成比の季節変動と赤潮発生の関係．
海岸工学論文集，49，1081‑1085．
陸田秀実，土井康明，河村和憲，亀山直記，2001：広島湾奥部の赤潮発生に及ぼす陸域起源の降雨量の影響．海岸工学論文集，48，1071‑1075．

Akiyoshi Nakayama
多田邦尚，和田彩香，一見和彦，橋本俊也，2006：海砂利採取船からの高濁度排水中の微粒子の挙動─
微粒子の特性と沈降速度─．
沿岸海洋研究ノート，
43⑵，157‑162．
武岡英隆，速水祐一，兼田淳史，2005：瀬戸内海における環境の長期変動．沿岸海洋研究ノート，
43⑴，45‑50．
武岡英隆，速水祐一，兼田淳史，松下太郎，紀本岳志，
渡辺浩三，藤川淳一，2001：瀬戸内海の栄養塩環
境の長期自動モニタリング．
沿岸海洋研究ノート，
38⑵，91‑97．
柳哲雄，2001：海洋環境モニタリングと情報ネット
ワーク─「せとうちネット」について─．沿岸海洋研究ノート，38⑵，99‑102．

流動と栄養塩，低次生産モデル

付録

Akiyoshi Nakayama流動と栄養塩，低次生産モデル

Akiyoshi Nakayama流動と栄養塩，低次生産モデル

Akiyoshi Nakayama