土木学会論文集B2（海岸工学），Vol. 74, No. 2, I_667─I_672, 2018.

半閉鎖性内湾における海水交換特性の季節変化とその形成メカニズムについて
岡田 信瑛1・内山 雄介2・多田 拓晃3・遠藤 颯4・馬場 康之5・
水谷 英朗6・久保 輝広7・森 信人8・渡部 靖憲9・大塚 淳一10・
山田 朋人9・猿渡 亜由未11・二宮 順一12
1 学生会員神戸大学大学院 工学研究科市民工学専攻（〒657-8501 神戸市灘区六甲台町1-1）
E-mail: eybuieybui@gmail.com
2 正会員 神戸大学教授 大学院工学研究科市民工学専攻（〒657-8501 神戸市灘区六甲台町1-1）
3 正会員 日本工営株式会社 名古屋支店技術部（〒461-0005 愛知県名古屋市東区東桜2丁目17-14）
4 非会員 伊藤忠商事株式会社 住生活カンパニー生活資材第三部（〒107-8077 東京都港区北青山2-5-1）
5 正会員 京都大学准教授 防災研究所（〒649-2201 和歌山県西牟婁郡白浜町堅田2500-106）
6 正会員 株式会社ハイドロ総合技術研究所（〒530-6123 大阪市北区中之島 3-3-23 中之島ダイビル 26F）
7 正会員 京都大学技術職員 防災研究所（〒649-2201 和歌山県西牟婁郡白浜町堅田2500-106）
8 正会員 京都大学准教授 防災研究所（〒611-0011 京都府宇治市五ケ庄）
9正会員 北海道大学准教授 大学院工学研究院（〒060-8628 札幌市北区北13条西8丁目）
10正会員 （国研）寒地土木研究所 寒地水圏研究グループ（〒062-8602 札幌市豊平区平岸1条3-1-34）
11 正会員 北海道大学助教 大学院工学研究院（〒060-8628 札幌市北区北13条西8丁目）
12 正会員 金沢大学助教 理工研究域 環境デザイン学系（〒920-1192 金沢市角間町）

半閉鎖性内湾である和歌山県田辺湾湾口部において，2016年9月から12月にかけて田辺中島高潮観測塔
計測システムおよび自記式ADCPを用いた長期連続観測を行い，海水交換特性に関する解析を行った．観測結果より，長期的には観測の前半期にあたる成層期は湾内から湾外への流出傾向，後半期の成層崩壊期は流入傾向にあることがわかった．流出期には台風16号が田辺湾を直撃し，瞬間的に強い湾外流出イベントが惹起された．流入期・流出期ともに流動のスペクトルでは半日周期成分が卓越したのに比べ，風速・水位とのコヒーレンスでは半日周期成分は小さくなり，亜表層もしくは混合層以深における傾圧的な流速分布により流動方向が逆転したことが示唆された．また，流出期には黒潮流路変動やメソスケール還流の影響に伴うサブタイダル周期での水温上昇イベントが複数回発生し，海水交換に強い影響を与えていた．

Key Words : semi-enclosed bay, water exchange, ADCP, Kuroshio and subtidal fluctuations

1. はじめに

中島高潮観測塔（以下，観測塔）を用いて，毎年夏に台風接近時の気象・海象データ取得を目的とした現地観測を断続的に実施し，これまでに2009年18号3)，2013年17号，和歌山県・田辺湾は紀伊半島南西端に位置する半開放系の湾である（図-1）．本海域には台風が頻繁に接近し， 18号，26号1), 4), ，2014年18号，19号5)などの台風接近イベ
特に秋季に散発的な気象擾乱の影響を強く受ける一方で， ント時の海象情報を得ることに成功した．特に，観測結
地形的には半閉鎖性であるという極めてユニークな条件

果と高解像ダウンスケーリング海洋モデルによる数値解下にある1)．また，紀伊半島南端沖を東進する黒潮に伴析を組み合わせることにより，湾口部における流動構造

う外洋の影響を間接的に受け，豊かな生態系が維持されと湾内外の海水交換特性に着目した詳細な解析を行ったている．そのため，田辺湾南部や奥部の静穏域は古来よ1)り養殖漁場として利用されてきたが，その半面，有機物至る約3ヶ月間の長期連続係留観測を行い，欠損の少な．それに加え，2015年には，成層期から成層崩壊期に負荷による赤潮や貧酸素水塊の形成がしばしば生じる2)． い良質なADCPデータを取得することに成功し，数値海したがって，汚濁負荷低減策と並行して，湾の海水交換洋モデルに基づく推算結果1)の整合性を確認した．しか特性を定量的に理解し，水質管理計画に応用する技術のしながら，台風接近イベント時を含むADCPデータを取構築が求められている．これに対して筆者らのグループ

得するには至っていなかった6)．本研究では， 2016年夏

は，田辺湾湾口部に位置する京都大学防災研究所・田辺季〜冬季にかけて，田辺湾湾口部に位置する観測塔に常I_667

図–1 観測海域（和歌山県田辺湾，黒枠）．カラー：水深．図中左上のインセットは，田辺湾周辺（黒枠部）の拡大図であり，インセット中の赤三角印は京都大学防災研究所・田辺中島高潮観測塔の位置を示す．

図–2 田辺中島高潮観測塔周辺における計測器の設置概況．

台風12号，台風16号の気象擾乱の影響を含んだ連続的で分と，鉛直方向上向き），流向（北向きから反時計周り）

ノイズレスなADCPデータの取得に成功した．本論文での時系列である．12月の気圧・風速データは欠測となっは，台風イベントを含む夏季〜冬季における田辺湾湾口ている．まず，ADCPによる東向き流速および流向から，部における海水交換特性に着目し，長期連続観測データ

観測期間前半の9月から10月下旬にかけては西向きの流解析を通じて湾口フラックスの時間変動・季節変動特性

速成分が強く，湾内から湾外への流出傾向にある（以下，を定量化し，その形成メカニズムを明らかにする．さら流出期と呼称）一方で，11月以降は東向きの流速成分がに海上風速に代表される外力と海水交換の周波数特性と強く，流入傾向にあった（以下，流入期と呼称）と判断の関係について着目した解析を行ったので，その結果をできる．水温の鉛直分布に着目すると，期間を通じて全報告する．

般的に徐々に低下しており，夏季（成層期）から秋季

系列変化を図-3に示す．上から順に，海面気圧，平均水設された波高・周期，潮位，海水温，海上風計測に加え， 位，海上風速ベクトル，海水温の鉛直分布，塩分（水深長期連続ADCP観測（2016年9/1〜12/7，図−1）を実施し， 30 m地点），流速鉛直分布（東向き，北向きの水平2成（成層崩壊期・非成層期）への季節的な長期変化が卓越

2. 現地観測している．成層期はおおむね流出期と対応しており，成層崩壊期は流入期と対応しているものと考えられる．前(1) 観測の概要

半の成層期において，台風12号，台風16号が接近した9田辺湾湾口に常設された京都大学防災研究所白浜海象

月8日，9月20日に気圧の大幅な低下と，水深30m付近に所田辺中島高潮観測塔（図-1，観測塔）およびその周辺おける水温の一時的な低下が見られる．その後，流出か海域において，2016年9月1日から2016年12月7日まで長ら流入へ転じる期間中である10月21日以降の2–3日間で期連続係留観測を実施した．観測塔における波高，潮位， は，水深10m以深において水温の急激な上昇が発生し，海水温，海上気象計測に加えて，ADCP（Acoustic Doppler
同時に塩分も上昇した．また，後半の非成層期には表層Current Profiler, 超音波流向流速計，Teledyne RDI社製付近で複数回の水温の上昇が見られる．これらの水温上
Workhorse Sentinel 1200 kHz）による流速分布計測，水深30昇イベントについては，3. (4)で検討する．後半期においmまでの海水温鉛直分布計測，水深 30 m地点での塩分計ては，表層の水温変動は見られるが，水深30m地点の塩測を行った（図-2）．湾口部の水深は平均 30 m弱である

分は時間的にほぼ一定であり，同様に水深30m付近の水が，観測塔は天端水深約10 mの小規模な孤立シルの上に温変動も前半に比べると非常に小さい．また，前半は風敷設されており，観測塔の海面上では，風向・風速，気向の反転がみられるが，後半は北向き風が卓越していた．

圧，波高などが観測されている．ADCPは水深 10 mの位以上のことから，本観測期間における田辺湾湾口部の置に上向きに設置し，水深10 mまでの流速分布を計測し海洋表層力学特性は，成層期から非成層期への季節的なた．なお，水深10 mまでの表層付近については，水温計表層混合層の発達と，卓越風向変化に連動した流動特性を集中的に配置した．連続観測の詳細については，馬場の変化で特徴づけられる．加えて，風に代表されるイベ7)

ら を参照されたい．

ント的な水温塩分低下と風速強化が重畳していたと考察される．そこで以降では，周波数スペクトル解析によっ(2)

観測結果て流速の周期数特性および外力との相関を評価するとと約3ヶ月分の連続観測データのうち，主要な項目の時もに，湾口通過フラックスを算定し，流出入量変動を詳

I_668

図–3 連続観測結果の時系列．上から，田辺湾口観測塔における海面気圧，平均水位，海上風速ベクトル，水温鉛直分布，塩分 (水深 30 m 地点)，流速鉛直分布 (上から東向き，北向き，鉛直方向上向き)，水平流向 (北向きから反時計周り)の時系列（2016 年 9 月 1 日〜12 月 7 日）．気圧・海上風速は観測塔頂部付近の気圧・風向風速計，水温は水深 30 m まで設置した海水温計，塩分は水深 30m 地点の塩分水温計，他は ADCP（設置水深 10 m）による観測値を使用．

細に評価する．さらに，通過フラックスを周波数分解す

外力変動の影響を含んでいると考えられる．同様の傾向

ることで，台風の間欠性を反映した乱流成分，内部潮汐は北向き流速に対しても確認された．鉛直流速についてに影響された潮流成分などの各周波数成分の寄与率，歪も定性的に同様の傾向が見られるが，表層付近よりも水度を求め，湾口流動特性について詳細な解析を行う．

深5m以深において顕著に見られ，内部界面変動（内部波）の影響をより強く反映していることが分かる．

3. 解析と考察(1)(2) 湾口通過フラックス
田辺湾の海水交換特性を評価するために，ADCPによ

流速変動特性に関するスペクトル解析

観測期間を10月20日より前の流出期（成層期）と，そる流速の湾口ラインを直交する成分を取り出し，海表面れ以降の流入期（非成層期）に分割し，それぞれについから水深10mまで鉛直積分することにより，単位幅表層て東向き・鉛直上向き流速の周波数スペクトル密度鉛直体積フラックスQとして評価した（図–5，流入を正とす分布を求めた（図–4）．東向き流速については，流出期， る）．なお，田辺湾の湾口ラインは概ね南北方向に約長周期帯において，全水深にわたってパワーが大きい．

3.5 km，水深は約30 mであり，一点の流速データで海水交換を代表させることはできないが，観測塔位置は水平両期間を比較すると，流入期では日周期帯のパワーが表地形の影響を受けづらい湾口ライン中央部付近であるこ層付近に集中し，周期3日以上の長周期成分が特に強化と，ADCPが計測する混合層内外の表層流は湾口流動のされている．前半の流出期においては，日周期帯の両側大部分を占めることから，Qは湾口フラックスを特徴づの帯域を含めて全体的にパワーが強く，様々な周波数のける量と考えることができる．

流入期ともに半日周期，日周期，および2日程度以上のI_669

流出期の期間平均Q は−0.036 m2/s，流入期は+0.104 m2/sと見積もられた．期間最大値は9月20日午後の台風16号接近時の−1.411 m2/sであり，瞬間的であるが非常に強い流出が惹起されている．さらにQを周波数分解し，各成分のRMS，およびRMS寄与率を求めたところ，最も寄与の高い長周期成分は流出期に32.9％，流入期に42.7%であった（表–1）．台風に伴う間欠的な外乱の寄与を強く反映する乱流成分（周期8時間以下）は，流出期に18.8%，流入期で14.7%であり，流出に対する気象擾乱の影響は有意であったと考えられる．一方で，台風12号接近時については，乱流成分の変動は確認できなかった．
周波数分解後の湾口通過フラックスの時系列分布を見ると（図–5下），流入期である10月下旬以降において，期間を通してsubtidal成分に明確な流入傾向が現れていることが確認できる．そこで，それぞれの成分に対して歪度を求めたところ，流入期のsubtidal成分は他の成分と比べて最も強い正の歪度（+0.6048）を持っていた (表–1)．
流出期のsubtidal歪度は+0.1154と流入期の1/4程度にとどまっており，湾口フラックス形成に対してsubtidal成分が極めて重要であることが理解される．つまり，日周期や潮汐周期で生じる連続的な外力というよりは，散発的に発図–4 周波数スペクトル鉛直分布．上 2 枚：東向き流速，下2 枚：鉛直上向き流速．それぞれ，上：流出期（観測前半の成層期），下：流入期（後半の非成層期）．
縦線は左から順に，3 日, 2 日, O1 （~25 時間）， K1（~24 時間），M2（~12 時間）周期に対応する．

生する現象や，長い時間スケールでゆっくりと変動する現象（外洋影響）の影響が強いことが示唆される．
(3) 外力との関係
各期間における湾口通過フラックスQの変動の形成要れる．流出期では，北向き風速は4~5日周期，東向き風因を定量的に評価するため，外力である北向き風速，東速は2~3日周期でQとの相関が高いことから，北風に伴向き風速，水位と通過フラックスQのコヒーレンスをそう東向きエクマン輸送と，局所的な東風による吹送流のれぞれ求めた（図−6）．流出期の短周期・乱流成分で風影響が両方とも強かったことが示されている．実際に，速と強い相関を有することから，台風の影響が再確認さ

流入期においては風向は北寄りの傾向が強く（図–1），

図–5 観測点における湾口通過フラックスの時系列．海表面から水深 10m までの鉛直積分値（単位幅，m2/s）．流入が正．
（上）1 時間平均値．（下）周波数分解された湾口通過フラックス．下段では，周期 28 時間以上の subtidal 成分（黒），18~28 時間の日周成分（マゼンタ），8~18 時間の半日周成分（緑），8 時間未満の乱流成分（青）に分解した．

I_670

表-1 湾口通過フラックスの各周期成分における流量変動の RMS，RMS 寄与率，歪度．
流出期

流入期

RMS

RMS 寄与（％）

歪度

RMSRMS 寄与率（％）

歪度

subtidal
日周期0.187432.87

0.15540.192042.66

0.60480.071312.50

0.00450.059013.110.0117

半日周期

0.204335.83-0.0268

0.133029.55-0.0391乱流

0.107218.80

0.05200.066114.69-0.0379

図–6 観測塔における各外力と湾口フラックス Q とのコヒーレンス．横軸は cpd（cycle per day）．左から順に，北向き風速とQ，東向き流速と Q，水位と Q．黒：流出期（観測前半），マゼンタ：流入期（後半）である．

れは，Qが水深10mまでの流速の鉛直積分であるため，表層から下方に向かうにつれて流入出の方向が反転するような傾圧流の影響を反映した結果，半日周期の流量が鉛直方向に相殺され，単独の外力とQとの直接的なコヒーレンスが小さくなったものと推察される．
(4) 外洋影響
水温鉛直分布からは，流出期から流入期へ遷移する時期である10月21日頃に，水深30m付近の水温と塩分がと図–7 関東・東海海況速報（三重県水産試験場）による海もに上昇した．当時，黒潮は潮岬において接岸流路を取表面水温分布．2016 年 11 月 21 日．っていたため，高温高塩分である黒潮水塊の影響を反映
エクマン輸送による東向き輸送（流入）が卓越していたしたことものと推察される．その後の塩分変動は小さかと考えられる．流入期の風とQの相関の大小は，2~5日ったが，海面冷却の進展に伴い水温は徐々に低下した．

周期においては東向き，北向きともに流出期と反対の関

水深30m付近での弱い水温変動傾向から，流動もまた弱係にあり，外力への応答特性の差がQの差を生じたと考かったと考えられる．一方，図–5より，黒潮水塊の流入えられる．なお，日周成分の風との相関も高く，海陸風が生じたと思われる同時期の表層付近は流出傾向にあり，モードに代表される気象場の日内変動の影響がQ形成に表層と亜表層での流動特性に差異が生じたことが示唆さ寄与していたことが確認できる．この日周期の風とQのれる．また，11月中に複数回水温の上昇が確認されるが，相関は非成層期である流入期のほうが若干高い．流出期表層付近での現象であることから，黒潮水塊の流入の直である成層期には強い海陸風モードが表れるはずである

接影響を反映したものではないと考えられる．

が，乱流成分や数日周期成分の方がより強かったため，

内山ら1) は，田辺湾の長周期の海水交換と湾内環境は，相対的に小さく評価されたものと推察される．また，水紀伊水道と黒潮域に挟まれた海域に発達するメソスケー

位とQの関係を見ると（図–6右），日周期と長周期でのル還流の影響を強く受けていることを報告している．

相関が高いことから，日周期の内部潮汐と，成層や流動2013年に発生した流入イベント時においては，紀伊半島場の季節消長・黒潮流路変動等による長周期の外洋条件

南西部において反時計周りの中規模渦が生じ，その補償変化の影響が重畳して表れたことが分かる．

流として紀伊半島側沿岸部に負の渦度が生じ，紀伊半島一方，流入期・流出期ともに半日周期成分の流速は全西岸に沿って強い南下流が発達したと述べている．今回

水深にわたってパワーが大きく（図–3），順圧的な潮流の観測結果より，subtidal成分による流入（図–5）と，南の影響が現れたと推察される．しかしながら，その鉛直向き流速の発達（図–4）は同時期に発生しており，2013

積分値であるQについては風速，水位ともに半日周期で年と同様の現象が発生した可能性が示唆される．11月21の相関が低く，流入期の方がその傾向が顕著である．こ日から22日にかけて起こった流入に着目すると，田辺湾

I_671

の北西部における表層水温分布が23~23.5℃となってお参考文献

り（図–7），今回の観測結果と整合的であるといえる．

1)

4. おわりに本研究では，田辺中島観測塔および ADCP を用いた

2)

約 3 ヶ月間の連続観測を行い，9 月から 10 月下旬の湾内から湾外への流出期と，10 月下旬から 12 月上旬にかけての流入期に分割した解析を行った．台風による気象擾

3)

乱，黒潮流路，気象場の季節変動に伴う長周期外力変動の影響を最も強く受け，ついで内部潮汐に伴う短周期変4)動が寄与していた．主要な結論は以下の通りである．
(1) 流出期の台風 16 号接近時には瞬間流出量が 1.411 m2/sとなり，期間最大値を記録した．流入期は期間を通してsubtidal 成分による安定した流入が続いていたのに対し，

5)

流出期には，台風などの突発的なイベントの影響が相対的に強く現れた．(2) 水深 10m までの東向き・北向き流速の周波数スペクトルには半日周期のパワーが強く表れていたが，これを鉛直積分した湾口通過フラックスでは， 6)半日周期の外力との相関が小さかった．これは観測地点における流速分布の傾圧性による流向の反転により表層と亜表層の流速が相殺された結果であると推察された．
(3) 流入期へ切り替わる時期，および流入期全体にわたり暖水の流入が見られた．これは主に外洋影響，特に黒7)

潮が接岸流路を取った結果，紀伊水道スケールでのメソスケール還流の影響を受けたものと考察された．
謝辞：本研究は科学研究費（15H04049，15KK0207，

内山雄介，松川大佑，神吉亮祐，馬場康之，森信人，
水谷英朗，渡部靖憲，大塚淳一，山田朋人，猿渡亜由未，二宮順一: 紀伊半島田辺湾湾口部における海水交換特性に関する研究，土木学会論文集 B2 (海岸工学)，Vol.70，No.2，pp.I_446-I_450, 2014.
上出貴士: 和歌山県沿岸における赤潮 – 1981 年から2000 年の赤潮発生状況 - ，和歌山県農林水技セ研報，Vol.5, pp.107-116. 2003.
森 信人，鈴木崇之，木原直人：海洋表層鉛直混合におよぼす風応力と波浪の影響，土木学会論文集 B2（海岸工学），Vol.66, No.1, pp. 311-315, 2010.
馬場康之，水谷英朗，久保輝広，内山雄介，森信人，渡部靖憲，大塚淳一，山田朋人，猿渡亜由未，二宮
順一: 台風通過に伴う田辺湾湾口部における海水温変化について，土木学会論文集 B2 (海岸工学) ，Vol.70，I_476-I_480. 2014.
久保慎也，二宮順一，森 信人，馬場康之，水谷英朗，
久保輝広，内山雄介，渡部靖憲，大塚淳一，猿渡亜由未：現地観測に基づく水温躍層の破壊と鉛直混合に関する研究，土木学会論文集 B2（海岸工学），
Vol.72, No.2, pp.I_487-492, 2016.
内山雄介，吉木智軌，多田拓晃，馬場康之，水谷英朗，久保輝広，森信人，猿渡亜由未，大塚淳一，二宮順一，渡部靖憲，山田朋人: 半閉鎖性海域湾口部における成層期および成層崩壊期の流動構造と海水交換について，土木学会論文集 B2（海岸工学）, Vol.
72, No. 2, pp. I_493-I_498，2016．
馬場康之，水谷英朗，久保輝広，内山雄介，森信人，渡部靖憲，山田朋人，猿渡亜由未，大塚淳一，二宮
順一: 田辺中島高潮観測塔における 2016 年夏季観測結果，京都大学防災研究所年報, 第 60 号 B, No.2,pp.664-671, 2017.
(2018. 3. 15 受付)

18H03798，代表：内山雄介），京都大学防災研究所一般
共同研究費の援助を受けた．

SEASONAL VARIABILITY OF WATER EXCHANGE AT THE ENTRNCE OF
A SEMI-ENCLOSE BAY AND ITS UNDERLYING MECHANISMS
Nobue OKADA, Yusuke UCHIYAMA, Hiroaki TADA, So ENDO, Yasuyuki BABA,
Hideaki MIZUTAMI, Teruhiro KUBO, Nobuhito MORI, Yasunori WATANABE,Junichi OTSUKA, Tomohito YAMADA, Ayumi SARUWATARI, and Junichi NINOMIYA
An extensive field campaign was carried out to examine water exchange at the entrance of the semienclosed Tanabe Bay, Wakayama, from September until December 2016 using the observing system of
the Tanabe-Nakashima observatory and an ADCP. The first half of the observation period mostly understratified condition is characterized by a gross outflow from the bay, while the second half under unstratified condition is dominated by inflow into the bay. Typhoon 201616 attacked the bay during the first half,
which excited intensive instantaneous outflow. Although the semi-diurnal component prevails in the upper ocean velocity during both periods, the coherent analysis indicates that the external forcing such as
wind and sea surface height is less correlated at this frequency. This suggests that baroclinic velocity profiles at depth would yield inversion of flow direction that cancels the vertically-integrated volume flux,
which is viewed as a proxy of the water exchange. During the second half, the subtidal low-frequency
component is responsible for the exchange, which may be caused by meridional shift of the Kuroshio pathand associated mesoscale circulations formed in the Kii Channal that leads to multiple warming events.

I_672