土木学会論文集B2（海岸工学），Vol. 73, No. 2, I_1315─I_1320, 2017.

琉球諸島周辺海域における生態系ネットワーク
形成に対する黒潮の影響について	 
小谷 瑳千花1・内山 雄介2・鹿島 基彦3・上平 雄基4,1・御手洗 哲司5
1学生会員	  神戸大学大学院 工学研究科市民工学専攻（〒657-8501 神戸市灘区六甲台町1-1）

E-mail:169t107t@stu.kobe-u.ac.jp
2正会員	  神戸大学教授 大学院工学研究科市民工学専攻（〒657-8501 神戸市灘区六甲台町1-1）
3 非会員	  神戸学院大学准教授 人文学部人文学科（〒650-8586 神戸市中央区港島 1-1-3）
4学生会員	  日本原子力研究開発機構 原子力基礎工学研究センター 環境動態研究グループ（〒319-1195 茨城県那珂市東海村白方白根2-4）
5沖縄科学技術大学院大学准教授 海洋生態物理学ユニット（〒904-0495 沖縄県国頭郡恩納村字谷茶1919-1）

JCOPE2-ROMS ダウンスケーリングシステムによる 2 段ネスト琉球諸島海域流動モデリング及びサンゴ卵・浮遊幼生を想定した 3 次元 Lagrange 粒子追跡数値実験を行い，コネクティビティの評価と黒潮に捕捉されて長距離輸送される粒子（LDP）に着目した広域分散解析を実施した．粒子の多くは放卵海域内及びその周辺に漂着するが，黒潮により北東方向へ輸送される LDP は放卵粒子総数の約 1 割に達し，サンゴ生態系の島間リンク形成に大きく寄与することを示した．石西礁湖北部海域を放出源とした粒子は常にLDP に寄与するが，南部海域起源の LDP 発生確率は放卵時の経度方向流速に強く影響されていた．また，黒潮により北東方向へ輸送された LDP は，琉球諸島と黒潮間に発達するメソスケール渦に伴う黒潮反流により黒潮を離脱し沖縄本島へ接近することが示唆された．

Key Words : Kuroshio, Ryukyu Islands, coral spawn and larvae, connectivity, ROMS

1.	  はじめに	 	  琉球諸島海域では，北進する黒潮に伴う低緯度域からの暖水輸送のため，比較的高緯度にありながら豊かなサ	  サンゴは海洋生物への生息場所の提供を通じて生物多ンゴ生態系が維持されている．黒潮から横断方向に 100様性維持に貢献するとともに，その造礁機能によって天km 以上離れた琉球諸島沿岸部への暖水波及メカニズム

然の防波堤として海岸への波浪の作用を軽減し，津波やとして，黒潮と沖縄本島間に発生する負のサブメソスケ台風の被害を抑制する．また，共生する褐虫藻の光合成ール渦に伴う eddy heat flux3)，黒潮−沖縄本島に間欠的にによって海水中の CO2 濃度のバランスを保つとともに， 生じる負のメソスケール渦に伴う黒潮反流 4), 5) などが知観光資源として地域経済に大きく寄与する側面も持ってられている．このうち黒潮反流は，春季に強化され冬季いる．しかし，温暖化などの影響によって世界規模でサに弱化するように季節変動している 4)．したがって，黒ンゴの白化や衰退が生じ，深刻な環境問題と認識されて潮はそれに付随する時空間的に変動する広域流動を通じ久しい．1998 年には大規模な海水温の上昇によって世て，熱だけでなく海洋生物の卵・浮遊幼生，栄養塩の輸界各地でサンゴの白化・衰退が発生したが，琉球諸島海送を惹起し，生態系ネットワークの形成に大きな影響を域（図-1）も例にもれなかった．1998 年以降，本海域及ぼしていると考えられ，本海域内において黒潮を中心のサンゴ被度は大きく低下して被度 30 %程度が続いてとした海洋構造を把握することが，サンゴ生態系保全におり，2004 年に一時的に約 40%まで回復したものの，向けて重要となる．

根本的な回復にはほど遠い 1)．さらに，2016 年には夏季	  本研究では，琉球諸島海域における広域サンゴ生態系高水温による大規模な白化現象が再び発生し，とりわけネットワーク構造を精査するため，JCOPE2−ROMS ダウ

同海域南部に位置する石垣島，西表島および石西礁湖周

ンスケーリングモデル 3), 4), 6)をベースに，2 段ネスティン

2)

辺での白化率は全体の約 7 割に達し ，現在に至るまで

グによる 2 種類の流動モデリングを行い，サンゴ卵・浮本海域でのサンゴ生態系保全は依然として喫緊の課題の遊幼生を模擬した Lagrange 中立粒子の 3 次元追跡計算をままである．

実施した．石西礁湖周辺海域におけるサンゴの産卵は 5I_1315

トワーク構造について解析を行う．次に，2006 年と2007 年に実施した GPS 表層漂流ブイ追跡観測実験結果と，2012〜2015 年を対象とした Lagrange 粒子追跡実験結果から長距離移動する粒子に着目し，黒潮に捕捉されて長距離輸送される粒子の発生確率とその分散過程について解析を行う．

2.	  海洋モデル	  本研究では，JCOPE2 8) 海洋再解析データ（水平解像度約1/12 度）を最外側境界条件及び初期条件として与え，
領域海洋循環モデルROMSをベースとした2段階の1-wayoffline 海洋ダウンスケーリングにより2種類のROMS-L2

図-1	  JCOPE2 に内包された 2 段階ネスティング ROMS モデ
ル領域（L1：水平解像度約 3	  km，L2v1：同 600	  m，L2v2：同 1	  k	  m）．カラーは水深（m）	  	 モデルを構築した．まず，著者ら3), 4), 6) による日本南西の東シナ海及び北太平洋の一部を内包する中解像度のROMS-L1モデル（同3 km）1日間平均値出力を境界条件としてダウンスケーリングすることで，琉球諸島海域

表-1	  ROMS-L2v2 の計算条件
計算期間
格子数
外力（海上風）海面フラックス
海表面水温・塩分
河川
海底地形
潮汐
境界条件	  を対象としたROMS-L2v1モデル（同600 m）を構築した

2010/12/27 - 2015/11/2
1280×1120×32 層，水平解像度 1 km
JMA GPV-MSM
COADS（月平均気候値）
JCOPE2（20 日平均値）
月平均気候値（Dai and Trenberth ）SRTM30TPXO7.0
ROMS-L1（日平均値）（図-1）．L2v1モデルの海上風応力には，気象庁GPVGSMの再解析値を用いた．	  さらに，より広域をカバーするために琉球諸島を中心としたROMS-L2v2モデル（同1 km）を構築し，L1モデルからのダウンスケーリングを別途行った．このL2v2モデル（図-1，表-1）では，海上風応力にGSMより高解像度の気象庁GPV-MSM再解析値を用い，さらにTPXO7.0による主要10分潮調和定数を用いた潮汐を開境界で与えた．また，Dai and Trenberth (2009) 9)	  の月平均気候値を用いて揚子江などの主要河川を考慮した．両L2モデルではL1モデルと同様に，気候値を用いることにより生じる海面フラックス誤差に伴うドリフトを補正するために，JCOPE2再解析データの20日間平均値SSTとSSS値を用いて緩和補正した．	  スピンアップ期間を含む計算期間は，L2v1モデルは2010年12月27日から2013年8月12日（解析期間は2013年5月〜6月），L2v2モデルは2010年12月27日から2015年11月2日の約5年間（解析期間は各年5月〜6月）である．

3.	  コネクティビティ解析	  琉球諸島海域内のコネクティビティを定量化するため，	  ROMS-L2v1モデルの再解析結果を用いて3次元Lagrange

図-2	  琉球諸島周辺海域拡大図と水深分布．数字はコネクティビティ解析のための Lagrange 粒子のソース及びシンクパッチが位置するエリア番号（図 -3 参照）．黒枠：左上 inset に示す拡大図領域，赤枠：

中立粒子追跡計算を行い，移流時間毎のLagrangian PDF（粒子の変位に関する確率密度関数）を求めた．ここでは主な島の沖合い 3 kmまでの海域に計 19,049 個の月の満月大潮の翌日の夜間に行われるため，春季を粒子追跡計算の対象期間とした．まず，海域中の島を中心に18 のエリアに分割し（図-2），Lagragian PDF 7) に基づく海域内コネクティビティを定量化し，サンゴ生態系ネッ

Lagrange 粒子の放出源（放卵源）を等間隔に定義し，さらにこれらを島ごとに18のエリアに分類して整理した（図-2）．2013年5月1日から6月30日まで12時間おきにLagrange 粒子を表層（水深 1 m）より放流し，移流時間I_1316

図-4	  石西礁湖からリリースした GPS 漂流ブイ追跡観測
結果．黒線：各漂流ブイの軌跡，丸印：漂流ブイのリリース位置（赤），最終漂着位置（青），カラー：L2v2 による Lagrange 粒子追跡計算結果の移流時
間 21 日間の PDF 時間積分値（4 年平均値），左上
inset：点線領域の拡大図．	  	  	  図-3	  	   ROMS-­‐L2v1 モデルを用いたコネクティビティ解析
結果（移流時間 30 日）．縦軸：Lagrange 粒子のシンクパッチ位置のエリア番号，横軸：ソースパッチ位
置のエリア番号．エリア番号と位置は図-2 に示すとおりである．	  

が30日を経過するか，粒子が着岸（陸グリッドに到達）あるいは計算領域外へ流出するまで3次元追跡計算を行

4.	  GPS漂流ブイ観測との比較った．計算期間の2ヶ月間にリリースした総粒子数は約	  琉球諸島南西部の石西礁湖と沖縄本島間の長距離コネ2,286,000 個である．粒子追跡計算方法およびLagragian

クティビティの形成機構を把握するため，ROMS-L2v2モ

PDF, Connectivity matrixの算出方法などについては，Mitarai

デルの出力結果を用いた3次元Lagrange 粒子追跡計算実ら7)，内山ら10) を参照されたい．

験を行った．ここでは，石垣島および西表島の沿岸，石	  ROMS-L2v1出力＋粒子追跡から求めたConnectivity ma-

西礁湖海域に半径 1 kmの円形パッチを合計145個設置し，trixを図-3に示す．全18エリアにおいて，移流時間30日間の粒子はリリースされた各エリア内に留まる傾向が強

各パッチ内の水深下3 m地点に約35 m間隔で放卵源を定が，これは北東への輸送が卓越すること，すなわち海域内を南西から北東方向へ通過する黒潮によって南西エリ4カ年の春季（すなわち4ケース）を対象に粒子追跡計算を行った．各粒子は移流時間が21日を経過するか，着岸アの放卵域から北東エリアの島近傍海域へ粒子が輸送さあるいは計算領域外に流出するまで追跡した．すなわち，れたものと解釈される．例えば，サンゴが豊富に生息す（放卵14日＋追跡21日）× 4ケースの粒子追跡計算を行る海域である南西部の石垣島，西表島，およびそれらのった．各ケースにおける総放卵粒子数は約4,006,000 個で間に存在する石西礁湖（図-3，エリア15-17）から放卵ある．

された粒子が沖縄本島（図-3，エリア8）へ漂着するコ	  モデルの再現性を検証するために，石西礁湖においネクティビティは 10 -3 であることから，石西礁湖周辺でて行った GPS 漂流ブイを用いた観測結果と，2012 年か発生したサンゴの卵や幼生が沖縄本島へ直接輸送される

ら 2015 年を対象とした 4 ケースの粒子追跡計算結果か可能性は無視し得ないと解釈される．つまり，琉球諸島ら求めた Lagrangian PDF の移流時間 21 日間積分値の 4 ケ海域でのサンゴ卵・浮遊幼生の輸送に対しては，黒潮をース分の平均値との比較を行った（図-4）．観測は，中心とした海域の流動特性の影響が強く，放卵海域外へサンゴの産卵時期に相当する 2006 年 6 月 12 日，13 日おの流出はもとより，遠方のエリアまで漂流する確率も高よび 2017 年 5 月 31 日の深夜に合計 6 機の GPS 漂流ブイ

く（最長輸送距離は400 km以上），長距離移動するサンを石西礁湖から放流し，30 分おきに GPS 位置情報を取ゴの卵や幼生の分散過程が島間のリンク形成に大きく寄得した．漂流ブイの軌跡から，6 機のうち 4 機は，石垣与することが示唆される．したがって，琉球諸島海域全島周辺を時計回りに浮遊して石西礁湖周辺に漂着，ある義した．サンゴの産卵を模擬するために5月の満月大潮
いが，遠方エリアとのリンク構造も多く見られる．また， の翌日を中心とした前後1週間（計14日間）を放卵期間
傾き45度ラインよりも上部で比較的高い値を示しているとし，ROMS-L2v2モデル計算期間の2012年から2015年の体のサンゴ生態系コネクティビティを解明するためには， いは南方向へ輸送されたことがわかる．一方，2006 年 6
長距離移動する粒子に着目することが重要であると理解
月 12 日に放流されたブイは石西礁湖北部から黒潮に捕される．

捉されたのち北東方向へ輸送され，21 日後には九州南

I_1317

図-5	  年毎の Lagrangian	  PDF 時間積分値	  ．粒子放出期間は，	  (a)	  2012 年 4 月 30 日から 5 月 13 日，(b)	  2013 年 5 月 19 日から 6 月1 日，(c)	   2014 年 5 月 9 日から 5 月 22 日，(d)	   2015 年 4 月 28 日から 5 月 11 日である．黒線は ROMS-­‐L2v2 モデルの領域を示し，グレー線は流速 0.5	  m/s おきの表層流速絶対値のコンターである．	  

岸まで到達した．

度線を越えた粒子（Long distance particle	  以下，LDP と呼	  ROMS-L2v2 モデルを用いた粒子追跡実験の結果から称）を黒潮に捕捉された粒子と考えて，全粒子から求めた Lagrangian PDF の移流時間積分値の 4 年平均値をLDP を抽出した．各ケースのリリース期間である 14 日

見ると，石西礁湖からリリースされた粒子は，石垣島・間に放出された全粒子に対する LDP 発生の平均確率は，西表島及び石西礁湖周辺を通過する確率が非常に高く，2012 年は 16.8％，2013 年は 5.7％，2014 年は 8.5％，2015

観測結果の漂流ブイの軌跡とも概ね一致している．また， 年は 17.6％であり，年によって多少差があるものの，平
確率的には 2 桁程度低いものの，黒潮流軸に沿って粒子均的に全放卵粒子の約１割程度は LDP として黒潮に捕が通過する確率も無視できないことが示されている．以捉されて北上することがわかった．

上のことから，粒子追跡モデルは統計的に漂流ブイ追跡	  次に，LDP の発生確率と粒子の放出位置との関係を観測結果と良好に一致しており，十分な再現性があると調べるために，各 LDP を放流パッチ毎に整理し，LDP判断される．

発生確率の空間分布を求めた（図-6）．石西礁湖北側から放出された粒子は，南側から放出された粒子に比べ5.	  黒潮による長距離粒子輸送て LDP 発生確率が 4 ケースともに高い．これは，石西礁湖西部および北部を通過する黒潮に伴う中規模渦によ(1)	  粒子追跡計算結果	  2012 年から 2015 年までのケースごとの Lagragian PDF
移流時間積分値を図-5 に示す．PDF 積分値は，21 日間の移流時間内に粒子が一度でもそこを通過する累積確率密度を表しており，値が大きいほどその地点を通過する確率が高いことを意味する．各年において放出源である
石西礁湖周辺の値が高く，放出源から離れるほど値が小さくなることから，石西礁湖を中心に琉球諸島海域内で粒子が広域的に分散している様子が確認できる．全ケー
スに共通して PDF 積分値は石西礁湖から北東方向へ偏在するように分布しており，移流時間 21 日間に黒潮の影響を受けて黒潮流軸に沿った北東方向への長距離輸送が生じていたことが分かる．同時に年差も見られ，例えば石西礁湖付近では 2012 年と 2014 年は全体的に南側境界に向かって輸送される粒子が多いのに対して，2013年および 2015 年は石西礁湖東部で PDF 積分値が高い．
つまり，石西礁湖から放出された粒子の輸送傾向には各年で有意な差が見られることから，粒子はリリース時や浮遊期間における周辺の海流の影響を鋭敏に受けていることが確認される．

る分散作用によって，北側パッチ群からの放流後，粒子は速やかに黒潮に捕捉されやすいことを示している．一方で，石西礁湖南側のパッチ群から放出された粒子のLDP 発生確率は，2013 年と 2014 年は非常に低く，2012年と 2015 年では比較的高い．このケース間の差異から，リリース時の石西礁湖南部の流動環境が LDP 発生確率に影響を及ぼしていることが予想される．	  そこで，放卵時の放卵源周辺の表層流動と，LDP 発生確率の時間変化を 4 ケースについて比較した（図−7）．石西礁湖南部（図-2，赤枠）のリリース時の空間平均流速（北および東方向）成分の日平均値はリリース日によって大きく変動しており，全粒子数に対するLDP 発生確率と流速（特に北向き成分）との関係は必ずしも明白ではない．しかしながら，LDP 発生確率がほぼ 0％となる 2013 年のリリース期間後半と 2015 年のリリース期間前半では，明らかに東向き流速が発達しており，放流された粒子が直ちに石西礁湖東側海域へ輸送されることで，黒潮に捕捉される LDP 数が大幅に低下したものと解釈される．この東向き粒子輸送は Lagragian
PDF 積分値分布（図-5）にも見られ，該当する 2 ヶ年(2)	  黒潮による輸送

では石西礁湖東側海域で PDF 積分値が明らかに高くな	  石西礁湖から放出された粒子のうち，黒潮の影響を受っている．つまり，石西礁湖北側から放出された粒子はけて北東方向へ長距離輸送される輸送パターンについて常に黒潮に捕捉されて LDP になりやすいのに対し，石着目する．ここでは，放卵源から十分に離れた北緯 26.5

西礁湖南側から放出された粒子は，放出直後の経度方向

I_1318

図-6 	  リリースパッチ毎の LDP の発生確率（パッチ内全放流粒子に対する LDP の割合）．丸印は西表島・石垣島および石西礁湖周辺に設置した合計 145 個のパッチ位置を表す．(a) 2012 年，(b) 2013 年，(c) 2014 年，(d) 2015 年．	  
図-7	  リリース日毎（横軸）に整理した LDP 粒子の全粒子数に対する発生確率（％，青線）とリリース時の時間前後 1 時間平均表層流速の石西礁湖南側海域（図-2 赤枠）に対する空間平均値の東方向成分(m/s，緑線)，北方向成分	  (m/s，赤線)の時系列．横軸はリリース日を示す．(a)	  2012 年，(b)	  2013 年，(c)	  2014 年，(d)	  2015 年．	  

図-8 黒潮反流にトラップされ沖縄本島に接近した粒子に対する Lagragian PDF 移流時間積分値．マゼンダ円形
枠を通過した粒子を黒潮反流捕捉粒子と仮定して抽出した．グレー線：各年計算期間平均表層流速の 0.5
m/s おきのコンター，青線：黒潮反流通過流量換算のための検査線．(a) 2012 年，(b) 2013 年， (c) 2014
年，(d) 2015 年．

流の影響を強く受け，特に東向き流れが発達するときはLDP 発生確率が著しく低下することがわかった．

表-2	  各ケースにおける検査線断面（図-8 青線）での黒潮反流流量（Sv）と検査海域（図-8 マゼンタ枠）を通過した粒子の全放流粒子数に対する割合（%）の比
較．

6.	  黒潮反流の影響	  著者ら4) が示したように，サンゴ産卵期である春季は，沖縄本島と黒潮流軸間において負のメソスケール渦に伴う黒潮反流が間欠的に発達する期間であり，LDPはメソ黒潮反流 (Sv)
粒子割合 (%)

スケール渦に伴う水平分散作用と，黒潮反流による南下流の影響を同時に受けて，沖縄本島方向へ輸送されることが考えられる．そこで，沖縄本島と黒潮間に半径502012 年2013 年

2014 年2015 年

3.94
0.70

6.24
0.03

7.02
1.31

17.212.89

過する流速を正と定義し，海表面から水深 600 mまで鉛直積分して通過流量の検査線方向分布を求め，さらに沖kmの検査領域（図-8，マゼンダ枠）を設け，そこを通過した粒子を黒潮反流に捕捉された粒子と仮定して，該

縄本島−黒潮間に発生する負の通過流量を検査線方向に当する粒子のみを抽出し，対応するLagragian PDF 移流時

いて，黒潮反流通過流量と黒潮反流に輸送された粒子数間積分値をケース（年）ごとに求めた（図-8）．2012年，2014年および2015年では，黒潮に捕捉されて北東方向へ輸送された粒子が沖縄本島の北側海域でメソスケール渦にトラップされ，黒潮反流によって南向きに輸送されていることがわかる．そこで，黒潮反流の通過流量のケースごとの35日間平均値を求め，黒潮反流に輸送された粒子数との比較を行った（表-2）．ここでは沖縄本島及びその西海域上に検査線（図-8(b)，青線）を設け，その断面内における検査線直行方向（概ね北東方向）へ通積分することで黒潮反流流量を評価した． 2013年を除には線形関係が見られ，反流量が大きいと粒子の割合も大きくなる傾向にある．したがって，黒潮反流の発達は沖縄本島へ接近する粒子数の増加に寄与していることが示唆される．一方，例外的に反流量が 6.2 Sv と比較的高いにも関わらず，沖縄本島への接近粒子が少ない2013年
は，そもそもLDPの発生割合が他ケースに比べ少なかっ
たこと，また，黒潮に捕捉されて北上してもメソスケール渦の分散作用をあまり受けず，反流に到達しなかったことなどが原因であると考えている．	 

I_1319

タリングサイト 1000 サンゴ礁調査報告書 , 151 p., 2017.
環境省那覇自然環境事務所：西表石垣国立公園 石西礁湖のサンゴ白化現象の調査結果について, 2017．	  琉球諸島海域内における広域的なサンゴ生態系ネット
3) Kamidaira, Y., Uchiyama, Y. and Mitarai, S.: Eddy-induced transport of
ワーク構造と，それに及ぼす黒潮の影響に着目した数値
the Kuroshio warm water around the Ryukyu Islands in the East ChinaSea, Cont. Shelf Res, 2016 (available online).
解析を行った．まず，3次元 Lagrange 粒子追跡から島間4)
内山雄介,
小谷瑳千花, 山西琢文, 上平雄基，御手洗哲司：のコネクティビティを定量的に評価し，粒子は放出され黒潮暖水波及に伴う沖縄本島周辺海域における非対称海たエリア内に漂着する確率が高い一方で，黒潮に捕捉さ洋 構 造 の 形 成 機 構 . 土 木 学 会 論 文 集 B2 ( 海 岸 工れて北東方向へ長距離輸送される確率も高いことが示さ学), Vol.72,No2, I_481-I_486,2016.
5) Qui, B. and N. Imasato: A numerical study on the formation of the Kuroれた．粒子輸送パターンは，石西礁湖から放流したGPS
shio Counter Current and the Kuroshio Branch Current in the East China
漂流ブイ追跡観測結果と統計的に整合的であった．次に，Sea, Cont. Shelf Res., Vol. 10, No. 2, pp. 165-184, 1990.
石西礁湖起源粒子を例に，広域ネットワーク形成に寄与
6) 内山雄介，上平雄基，御手洗哲司，坂上太一郎：琉球諸する長距離輸送粒子（LDP）を抽出し，解析を行った．
島周辺海域におけるサブメソスケール乱流混合による黒潮波及効果に関する研究，土木学会論文集 B2（海岸工石西礁湖北側から放出された粒子は常に黒潮に捕捉され学），Vol.69, No.2, I_471-I_475, 2013．
やすいものの，南側からの放流粒子は放流時に東向き流
7) Mitarai, S., D. Siegel, J. Watson, C. Dong, and J. C. McWiliams: Quanti速が卓越するとLDPが大きく低減することが示された．
fying connectivity in the coastal ocean with application to the Southern
さらに，黒潮に捕捉されて北東方向へ長距離輸送される
California Bight, J. Ge- ophys. Res., Vol. 114, C10026, doi:10.
1029/2008JC005166, 2009.
粒子の一部は，琉球諸島と黒潮流軸間に間欠的に発生す8)
Miyazawa, Y., R. Zhang, X. Guo, H. Tamura, D. Ambe, J-S. Lee, A.
る黒潮反流に再捕捉されて黒潮から離脱し，沖縄本島方Okuno, H. Yoshinari, T. Setou and K. Komatsu: Water Mass Variability
向へ輸送されることが示唆された．
in the Western North Pacific Detected in 15-year Eddy Resolving Ocean
Reanalysis, J. Oceanogr., Vol.65, pp. 737-756, 2009.
9) Dai, A., Qian, T., Trenberth, K.E. and Milliman, J.D.: Changes in conti謝辞：本研究は科研費補助金および基金（15H04049，nental freshwater discharge from 1948-2004, J. Climate, Vol. 22, pp.
15KK0207，17K00653）の援助を受けた．また，沖縄科
2773-2791, 2009.
学技術大学院大学の計算機 SANGO を一部使用した．
10) 内山雄介,小硲大地,西井達也,御手洗哲司:コネ クティビテ
ィを用いた瀬戸内海全域における海洋生態系ネットワー
ク構造の解析．土木学会論文 B2 (海岸工学)，Vol. 70, No. 2,
参考文献
I_1076−I_1080, 2014 .	  	  	  	  	  	  	  	  	  	  	  	  	  	  	  	  	  (2017. 3. 15 受付)
1) 環境省自然環境局 生物多様性センター：平成 28 年度 モニ7.	  おわりに2)

INFLUENCE OF THE KUROSHIO ON DISPERSAL OF CORAL SPAWN AND
LARVAE AROUND RYUKYU ISLANDS
Sachika ODANI, Yusuke UCHIYAMA, Motohiko KASHIMA,Yuki KAMIDAIRA and Satoshi MITARAIRyukyu Islands in the East China Sea are in a subtropical climate, hosting desirable environment for abundant coral ecosystem. It is anticipated that the Kuroshio and associated eddy mixing promote the transport and dispersal of coral spawn and larvae across the islands, whereas thearea has suffered from coral bleaching in the recent decades. Toward optimal preservation andprotection of the coral habitats around Ryukyu Island, we conduct a double nested highresolution synoptic ocean modeling using ROMS with grid spacing down to 600 m coupled with
an offline Lagrangian particle tracking model. The modeled dispersal of coral spawn and larvae
released from 18 major islands and lagoons are examined to quantify connectivity using Lagrangian probability density functions (PDFs) of the Lagrangian particles among Ryukyu Islands. The PDF analyses suggest that the particles mostly remain near the released areas, while
about 10% of the particles are entrained in the Kuroshio to travel long distance, contributing tointer-island connectivity. Some of these particles are further transported northeastward, and thenpartially trapped by the Kuroshio Counter Current formed between the Kuroshio and the islands
to approach Okinawa Main Island. We reveal that eastward current occurred in the south theSekisei Lagoon shortly after the releases significantly suppresses the entrainment and the resultant long-distance transport.

I_1320