土木学会論文集B2（海岸工学），Vol. 78, No. 2, I_1009─I_1014, 2022.

領域海洋モデルを用いた洋上風力発電施設の適地選定システム開発に関する基礎的研究
玉井 駿壮1・内山 雄介2・Xu Zhang3・増永 英治4・辻 一洋5・山崎 秀勝6, 7
1 学生会員

神戸大学 大学院農学研究科 食料共生システム学専攻（〒657-8501 神戸市灘区六甲台町1-1）
神戸大学教授 大学院工学研究科市民工学専攻（〒657-8501 神戸市灘区六甲台町1-1）
E-mail: uchiyama@harbor.kobe-u.ac.jp (Corresponding Author)
3
Department of Atmospheric and Oceanic Sciences, Peking University（Beijing 100871, P. R. China）
4 正会員 茨城大学講師 地球・地域環境共創機構（〒316-8511 茨城県日立市中成沢町4-12-1）5(株)アルファ水工コンサルタンツ 技術開発センター（〒063-0829 札幌市西区発寒9条14丁目）
6
東京海洋大学名誉教授 大学院海洋科学技術研究科（〒108-0075 東京都港区港南4-5-7）7(株)アルファ水工コンサルタンツ 東京本社技術開発センター（〒104-0045 東京都中央区築地3-9-9）
2 正会員

カーボンニュートラル実現に向けて洋上風力発電施設の拡充が社会的課題になっているが，その敷設を巡る日本特有の問題として，建設に際して漁業権者の同意を得る必要性が挙げられる．合意形成促進のために施設設置によって強化される鉛直混合による海洋表層への栄養塩供給とそれに伴う「漁礁効果」を定量的に評価する総合的なシステムの開発を最終目的に，システムのコア部分を構成するROMS-NPZD広域流動・生態系予測モデルの概要を紹介し，洋上風力発電施設の適地評価に向けた海域ごとの空間評価方法について検討した結果を報告する．

Key Words : offshore wind farm, suitable site selection, biological productivity, ROMS-NPZD model

1. はじめに
近年，地球温暖化抑止のための低炭素電力の主力として，洋上風力発電などの再生可能な海洋エネルギー資源の利用がめざましく進展している．特に北海やバルト海では 2000 年代に入って急速に洋上風力発電ファームの設置海域が拡大し，2018年末までに欧州各国に設置された累計のタービン数は 4,543 基，累計出力 18,499 MW に達する 1)．洋上風力発電施設（offshore wind power generation
facility, or offshore wind farm：以下では OWF と略記）の設されているが 1)，事前に定量的な評価を行われた事例はない．そこで我々は，漁業協調および OWF 誘致に向けた合意形成促進のために，施設設置によって強化される鉛直混合による海洋表層への栄養塩供給とそれに伴う「漁礁効果」を定量的に評価する総合的なシステムを開発すべく，新エネルギー・産業技術総合開発機構（NEDO）の支援を受けて各種基礎研究を実施している．
開発中の評価システムは大きく分けて 3 つのサブシステムから構成されている．それらは，日本沿岸全体を解析対象とし，海洋環境科学的に見た適地選定機能と置にあたっては，施工性・コスト・管理運用・環境負荷などが課題となっているが，日本特有の問題として漁業OWF 周りの流動再現用の超高解像度モデルへの境界条件提供を行う「広域モデル」，OWF を解像して近傍場権者の同意を得る必要性が挙げられる．OWF の設置に適した我が国の沿岸 5 km 以内の海域のほぼ全ては共同の流動・生態系環境を解析する狭領域向けの「超高解像度モデル」，予測された生態系に対するパラメータ漁業権などの漁業権が付与されている区域である．そのため，発電事業者と漁業権者が双方ともに潤う Win-Win（Chl. a など）と実際の水産有用種を主体とした魚類の彙集効果を計測する「モニタリングシステム」である．

方式（メリット共有方式）が OWF 事業における漁業協調の基本的考え方である．漁業協調施策としては OWF

本稿はこのうち「広域モデル」の開発と運用の方向性について報告するものである．広域モデルのベースとなの多目的利用，例えば海洋牧場の整備など幾つかの方策があるが，最も達成可能であると考えられているのは
OWF 敷設に伴う魚礁効果による漁獲量の増加である．
実際に，OWF 敷設後の調査によって魚礁効果が報告

るのは 3 次元海洋循環モデル ROMS を用いた JCOPE 2ROMS ダウンスケーリングシステムに窒素ベースの物質循環・生態系モデルである NPZD（Nutrient Phytoplankton
Zooplankton and Detritus）モデルをカップリングさせた

I_1009

表-1 ROMS-NPZDモデル計算条件一覧．
計算期間2015年 1月 1日〜2015年 12月 31日

格子数784× 864× 鉛直 20層

風応力

気象庁 JMA-GSM （日平均値）

T-S nudging
海面フラックス

JCOPE2 再解析値(10 日平均気候値)

海底地形
境界条件

SRTM30_PLUS + J-EGG 500
JCOPE2再解析値（日平均値）

水平解像度

3.0 kmCOADS（月平均気候値）

析期間は2015年の1年間とした．

図-1 JCOPE2-ROMS ダウンスケーリングシステムによる
ROMS-NPZD モデル領域（JCOPE 2 から直接ネスティングさせた ROMS-L1 モデル領域：黒枠）．カラーは標高および水深 (m) である．

ROMS-NPZD モデル 2) である．すなわち本研究は，OWF設置に向けた合意形成促進のために，施設設置によって強化される鉛直混合による海洋表層への栄養塩供給とそれに伴う「漁礁効果」を定量的に評価する総合的なシステムの開発を最終目的に，システムのコア部分を構成する広域流動・生態系予測モデルの概要を紹介し，洋上風力発電施設の適地評価に向けた海域ごとのモデル再現性の空間評価方法について検討したので，その結果を報告するものである．

2. 研究方法(1) JCOPE2-ROMS流動モデル
本研究では，黒潮大蛇行等の間欠的な個別イベントを包含し，時々刻々変化する海洋流動やそれに伴う生態系（一次生産）の変動をできるだけ精緻に表現することを目指し，流動・生態系モデルともに可能な限り現実的なデータで駆動させる総観値（synoptic）解析を行う．流動モデルでは，JCOPE2海洋再解析データ（水平解像度約10 km）の日平均値を初期条件・境界条件として，領域海洋循環モデルROMSを用いた1-way offlineネスティン
グにより，解像度3 kmのROMS-L1領域へとダウンスケーリングを行う（図-1，表-1）．海上風応力には気象庁
GPV-GSMの日平均値を，その他の海面フラックスにはCOADS05の月平均気候値を与える．さらに，黒潮流路等メソスケールでの再現性を向上させるためにTS Nudging（γ = 1/20 day-1； Uchiyamaら3）を参照）を適用する．なお，同一グリッドを用いた試行的な気候値（climatology）計算結果から，本モデルのプロトタイプが対象海域の3次元海洋構造を良好に再現可能であることを確認済みである2)．ROMS-NPZDモデルでは，助走計算期間を除く解(2) NPZD生態系モデル
Fashamら4)による窒素ベースの低次生産NPZDモデルを用い，Gruberらの設定5) を参考に，北太平洋海域向けに各種パラメータの調整や計算条件の設定を行った．本NPZDモデルでは，栄養塩（N：硝酸塩NO3-，アンモニウム塩NH4+），植物プランクトン（P），動物プランクトン（Z），デトライタス（D）の5要素を考慮する．初期・境界条件にはWorld Ocean Atlas 2015（WOA）による硝酸塩（NO3-）の3次元分布の季節平均気候値（深層は年平均気候値），MODIS/Aquaによる海表面Chl.aデータの月平均値に対してMorel & Berthon6)による鉛直分布モデ
ルを適用して得られる3次元Chl.a分布，Chl.aに比例させた動物プランクトン濃度初期分布を与え，初期デトライ
タス，アンモニウム塩はゼロとした総観値計算を行う．
なお，低解像度のWOAでは沿岸域でのNO3-分布を十分に表現できないため，水深500 m以浅の海域に対しては水深500 mでのNO3-鉛直分布を岸に向かって水平方向に線形外挿した値を初期分布として与えた．初期条件に含まれる不確実性は1年間の助走計算中に解消され，2年目以降では各生態系変数は初期条件によらずに内因的な物質収支にしたがって変動することを確認している2)．
(3) 数値モデル再現性の評価指標
ROMS-NPZDモデルを含む数値モデルが現実の海洋状態をどの程度再現できているのかを客観的に評価する方法について検討する．モデルの評価法には実に様々なアプローチがあるが，本研究では，平面的に得られる人工衛星データとの比較を念頭に，二乗平均平方根誤差（RMSE）と，Willmottの適合度指標7)（skillと呼称）を用いることとする．両指標の定義は以下の通りである．

I_1010

1𝑅𝑀𝑆𝐸＝% ((𝑋!"#$% − 𝑋"&' )(𝑁𝑠𝑘𝑖𝑙𝑙＝1 −∑|𝑋!"#$% − 𝑋"&' |(∑(|𝑋!"#$% − 333333
𝑋"&' | + |𝑋"&' − 333333
𝑋"&' |)((1)(2)

(a) SST: MGDSST (obs.)(b) |us|: AVISO (obs.)(c) Chl.a: MODIS/Aqua (obs.)

[oC]

[m/s]

[µg/L](d) SST: ROMS-NPZD (model)(e) |us|: ROMS-NPZD (model)(f) Chl.a: ROMS-NPZD (model)

[oC]

[m/s]

[µg/L]

図-2 広域 ROMS-NPZD モデルの再現性確認結果の一例（2015 年 1 年分の平均値）．a, d は海表面水温 SST（℃），b, e は表面流速絶対値 |us|（m/s），c, f は海表面 Chl.a 濃度（µg/L）．上段は対応する衛星観測データであり，a：気象庁 MGDSST による
海表面温度，b：AVISO 海面高度計データによる表層地衡流速，c：MODIS/Aqua 海色センサーによる海表面クロロフィル a（Chl.a）濃度．下段は対応する ROMS-NPZDモデル計算値．

ここに，𝑋!"#$% ：モデル出力， 𝑋"&' ：対応する観測デ333333
ータ，𝑋"&' ：観測データの空間平均値， N：比較を行うデータの総数であり，総和SはN個の要素に対して求める．RMSEは比較するデータと同じ次元を有し，値が小さいほどモデルの精度が良いことを表す．一方，skillは無次元量であり，1に近いほどモデルと観測値の一致度が高く，一般的には0.6程度以上で高精度とされる．
以下の解析では，観測データと計算結果の時刻を統一した（例えば2015年1年間平均値をそれぞれに対して求めた）うえで，観測データをROMS-NPZDの計算格子点に線形写像し，陸域を除外した全海洋計算格子点に対してRMSEとskillを計算し，評価に供した．

3. 計算結果と考察(1)モデルの再現性検証：現地観測データセット
ROMS-NPDZモデルの再現性評価を行うために，人工衛星データと，船舶観測による気象庁東経137度定線における水温・塩分データを用いる．衛星データは，気象庁によるMGDSST海表面水温データ（解像度1/4度），ECMWF CopernicusによるAVISO海面高度偏差データから理論的に計算される表層地衡流速分布（解像度1/3度），NASAのAqua衛星MODIS（中分解能撮像分光放射計）海色センサーから推定されるChl.a（約5 km間隔に格子化されたデータセット）を用いる．MGDSSTとAVISOはいずれも複数衛星の観測値に対して最適内挿法を適用して客観解析・格子化されたデータである．東経137度定線データは年数回の観測船によるCTDキャスト観測結果を季節平均気候値として独自に取りまとめた鉛直断面データである．
（2） モデルskill・RMSEによる再現性評価
ROMS-NPZDモデルの再現性を確認するために，2015年1年分の計算結果（年平均値）を対応する衛星観測データと比較した（図-2）．海表面水温（SST）と表面流速の絶対値に関しては黒潮流路やそれに付随する水温フロント構造などの大規模な空間構造をよく再現しており，全体的に低波数成分の一致度は極めて高い．一方，海面Chl.a濃度に関しては，太平洋沿岸，中国・ロシア沿岸，東北・北海道海域などでChl.aが増大するなどの全般的な空間パターンについては表現できているが，黒潮流路，黒潮続流域およびその北側などでChl.aを過大評価している．反対に日本海側では全体的に過小評価傾向にあり，北緯39〜40度のユーラシア大陸から本州にかけて発達する極前線周辺でのChl.aの再現性は必ずしも十分ではない．
次に，衛星観測データに対するROMS-NPZDモデルのRMSEと適合度指標skillを求めた（表-2）．海表面水温
SSTはRMSE：0.5982度，skill：0.9981，海表面流速は RMSE ： 0.1257 m/s ， skill ： 0.8933 ， Chl.a は
RMSE：0.8679 µg/L，skill：0.3228であり，特に物理場の再現性が非常に高いことが定量的に示されている．
一方で，図-2の空間分布から予見されるように，生物量であるChl.aの誤差は大きく，改善の余地がある．
鉛直構造の再現性確認のために，気象庁東経137度定線観測による水温・塩分データとの比較を行なった（図-3）．ここでは紙面の都合上，夏季（7−9月）に対する季節平均値の比較のみを示す．表層の季節

I_1011

表-2 衛星観測データに対する ROMS-NPZD モデルの再現性
指標値一覧．水温，流速絶対値，Chl.a 量は全て海表面での値である．

RMSE
skill

水温 SST

流速絶対値

0.5982 (oC)
0.99810.1257 (m/s)
0.8933(a) ROMS-NPZD

Chl.a
0.8679 (µg/L)
0.3228(b) JMA observation S [psu]T [oC]

表-3 気象庁東経 137度定線観測データに対する ROMS-NPZD
モデルの再現性指標値一覧（季節平均気候値）．
RMSE
skill

夏季水温0.4100 (oC)
0.9987

夏季塩分0.0350 (psu)
0.9843

冬季水温0.2862 (oC)
0.9989冬季塩分0.0331 (psu)0.9855成層の構造や岸方向への傾斜など，全体的な密度構造は両者ともに良好に一致していることが分かる．
図-3の結果に冬季の結果もあわせてRMSEとskillを求めた（表-3）．水温・塩分は夏季・冬季いずれにおいてもRMSEは非常に小さく，skillはほぼ1に近く，物理場に関しては水平分布だけではなく鉛直構造の再現性も極めて高いことが示されている．

図-3 東経 137度ラインにおける鉛直断面内の塩分（上段）と水温（下段）．（a）ROMS-NPZD モデル計算値，（b）
気象庁観測船による実測値，いずれも夏季（7−9 月）に対する季節平均気候値である．

ることを試みる．例えば，ブロック化された RMSEは以下のように定式化される．
/

$
,!01
!!*+,##*,!!*+,##*,'
𝑅𝑀𝑆𝐸%𝐱 !,# '＝( ∑+!01− 𝑋&-.
+2$ ∑,2$ %𝑋!&'()%（3） 生物量再現性向上への方策
生物量再現性の改善策としては，①生態系モデルの精緻化・複雑化，②境界条件の高精度化，③高解像度化などの方策が考えられる．①に関しては，本NPZDモデルの各種パラメータを調整したり，比較的低次の窒素収支モデルを超えた，より高度な物質循環モデルを用いること，②についてはWOAなどの気候値ではなく何らかのデータ同化再解析値を側方境界条件に用いることや，今回は考慮していない河川などからの陸域負荷を正確に導入すること，③については，例えば極前線構造の再現には水平解像度1 kmが必要と言われており3)，非地衡乱流の効果を導入すること，などの対策が精度向上つながると期待される．いずれも今後の重要な課題としたい．

4. 洋上風力発電発電施設適地選定に向けた評価(1) 新しいブロック評価指標の提案
3. で検討した観測データに対する RMSE や skill は，原義的には計算領域全体で空間平均された指標であり，計算領域全体でのモデルとデータの適合性を評価するものである．洋上風力発電施設の適地選定に応用するには，まずは再現性の空間分布を評価する方法を開発する必要がある．そこで本研究では，任意の水平座標を中心とした正方形ブロック状の分割領域に対して指標評価計算を行い，計算グリッド上に写像して空間的に表示・評価す(3)

ここで，𝐱!,* は ROMS-NPZD モデルの任意の水平格子点!,*!,*
座標， 𝑋!"#$%と𝑋"&'
は格子点番号(𝑚, 𝑛)における計算
値と観測値，𝑀 = 𝑖𝑚𝑎𝑥 ∙ 𝑗𝑚𝑎𝑥 , 𝑚𝑚 = 𝑚 − 𝑖𝑚𝑎𝑥/2,𝑛𝑛 = 𝑛 − 𝑗𝑚𝑎𝑥/2 である．skillについても同様な総和を取ることによりブロック化 skill を求めることができる．
この操作の意味するところは極めて単純であり，計算領域の全てではなく，任意の𝑖𝑚𝑎𝑥 × 𝑗𝑚𝑎𝑥のサイズにブロック化された分割領域を定義し，その中に含まれるデータだけで RMSE と skill を計算するということである．
以下の解析では，𝑖𝑚𝑎𝑥 = 100，𝑗𝑚𝑎𝑥 = 100とすることにより約300 km四方の分割領域を定義するとともに，格子点番号(𝑚, 𝑛)を 2 ステップずつ（6 km ずつ）移動させることにより，東西南北の全方向に対して 6 km おきにブロックをスライドさせながら各ブロックの周囲約300 km 四方についての RMSE と skill を求める．したがって，移動ステップや影響範囲を変化させることにより，任意の地点における任意の空間スケールでの「再現性」が評価可能となる．この考え方を RMSE や skill だけではなく適地性に関する様々な評価指標に対して適用することにより，モデルの結果を用いて洋上風力発電施設の適地性を定量的に空間評価することが可能になる．
(2) ブロック再現指標によるモデル再現性の空間評価
これまでモデルの RMSE や skill は表-2 のような評価しかできなかったが，ブロック指標の導入により再現性が空間評価できるようになる．まず，ブロック化された SST の skill 分布（図-4a）を見ると，全領域 skill

I_1012

(a) SST: skill(b) velocity |us|: skill(c) Chl.a: skill(d) SST: RMSE [oC](e) velocity |us|: RMSE [m/s](f) Chl.a: RMSE [µg/L]

図-4 図-2 の計算値と観測値から求められるブロック化 skill（上段 a-c）とブロック化 RMSE（下段 d-f）の空間分布. a, dは海表面水温 SST（℃），b, eは表面流速絶対値（m/s），c, fは海表面 Chl.a濃度（µg/L）
表-4 ブロック化されたモデル再現性指標 PI（PredictabilityIndex）の評価に用いた重みパラメータの一覧．

aibii = 1：SSTi = 2：流速絶対値

i = 3：Chl. a201020102515

Chl.a に対するブロック化 skill・RMSE 分布（図-4c, f）
から，黒潮流軸の南側領域で skill が低いことや，沿岸域で RMSE が大きく skill が小さいことが分かる．
Chl.a 分布計算値の比較（図-2c, f）より，計算値には約 0.2μg L-1 のバイアスがあるなど，再現性に課題が残されていることが改めて再確認される．
(3) ブロック結合指標の提案
OWF 設置海域の適性を ROMS-NPZD モデルの格子点ごとに数値化するための前準備として，モデルの再現性を総合的かつ空間的に評価する指数を提案する．ここでは将来実施する漁礁効果の計測を念頭に，水温・流速・Chl.a のブロック化 skill と，ブロック内平均値で無次元化した RMSE の重みつき線形和によって再現性指標（Predictability Index, PI）を定義する．

図-5 モデル再現性指標 PI（Predictability Index）の空間分布．
表-4 の重みを使用した．

値 0.9981 に対応して，大部分の海域で 0.9以上と評価されている．SST 分布の比較（図-4a, d）から，領域全体でskill が高い一方で，太平洋側の黒潮流路周辺で skill が低下するなど，ROMS-NPZD モデルの再現性の空間的な評価が可能となったことが確認される．
表層流速の比較（図-4b, e）を見ると，流速再現性は日本海側より太平洋側の方が良好であり，黒潮流路周辺で skill が最大値を取っている．しかしながらこの海域では RMSE が極めて大きくなっており，一致度（skill）は高いものの，流速値そのものが大きいため（図-2b, e）誤差の絶対値が大きくなることが示されている．また，

𝑃𝐼 = ( 𝛼+ 𝑆+ + ( 𝛽+ 𝐸+(4)ただし，式 (4) の下付き添字 i = 1, 2, 3 はデータの種類を表し，i=1：SST，i=2：表面流速，i=3：Chl.a である．また，Si は skill（0〜1）であるが，Ei は式 (5)のようにモデルと観測値の平均値で RMSE を標準化した無次元の指標で，平均値に対する RMSE の割合に関連し，この相対誤差が小さいほど 1 に近づく．
𝐸+ = 1 −𝑅𝑀𝑆𝐸
1∑(𝑋!"#$% + 𝑋"&' )2𝑁(5)ai, bi は Si と Ei に対する重みである．PI を 100 点満点で評価すること，物理量よりも生物量に，RMSE よりも skill に重みを置くことなどを念頭に，試行的に

I_1013

表-4 の組み合わせを用いる．これは単なる一例であり，当然ながらユーザーの PI に対する思想に応じて他の組み合わせを用いても全く構わない．
このような着想のもとに定義された PI の空間分布を求めたところ（図-5），太平洋側・日本海側を含む日本沿岸の大部分で PI は 60 点以上となり，関東沖で最大値を取る．物理-生物結合指標 PI は，広域モデルは良好な再現性を有しているものと評価される．

る海底勾配・流速振幅・沖合亜表層の高栄養塩層までの水平距離・鉛直混合強度などを合理的に組み合わせて指標化することなどが考えられる．このように，今後は本研究の成果を発展させ，どこの海域が OWF としての適地になりうるかを定量的に得点化して評価する枠組みに展開させていくことを考えている．
謝辞：本研究は新エネルギー・産業技術総合開発機構（NEDO）技術研究開発事業（課題番号：20356011），科学研究費（課題番号：18H03798）の援助を受けた．

5. おわりに
参考文献
本研究では，広域モデル・超高解像度モデル・モニタリングシステムを統合した洋上風力発電施設（OWF）の適地選定・評価システムのコア部分である広域モデルを開発し，現地観測データとの比較による再現性検証を通じて，海域ごとの空間評価方法を検討した．広域モデ
ルは JCOPE 2-ROMS ダウンスケーリング海洋循環モデルをベースに窒素循環モデルを統合した ROMS-NPZD モデルであり，日本沿岸の任意の海域での物理量および生物・化学量の定量化が可能となった．特に物理量の再現性は水平・鉛直方向ともに非常に高いものであった．
さらに，新たなブロック評価指標を提案し，計算モデルによる膨大なデータから様々なパラメータを抽出して格子点ごとに定量的に評価する簡便な枠組みを提示した．
提案されたブロック評価指標を応用することで，OWF構造物を敷設した際により高い漁礁効果が見込まれる海域を定量的に評価でき，さらにこれらを得点化することにより，適地を事前スクリーニングすることさえも可能になる．その際，例えば，施工や維持管理の容易さに関わる OWF に作用・関連する直接的な物理外力（水深・1)

2)
3)

4)

5)

6)

7)

風速・波高など）に加えて，生物・化学量分布に関連す中田英昭：洋上風力発電に伴う環境影響の調査と評価 ～欧州の事例をもとに今後を展望する～，Ocean Newsletter，
460号，海洋政策研究所，2019.
Uchiyama, Y., Suzue, Y. and Yamazaki, H.: Eddy-driven nutrienttransport and associated upper-ocean primary production along the Kuroshio, J. Geophys. Res. Oceans, Vol. 122, pp. 5,046–5,062, 2017.
Uchiyama, Y., Kanki, R., Takano, A., Yamazaki, H. and Miyazawa, Y.:
Mesoscale reproducibility in regional ocean modeling with a 3-D stratification estimate based on Aviso-Argo data, Atmosphere-Ocean, Vol. 56,
No. 4, pp. 212-229, 2018.
Fasham, M. J. R., H. W. Ducklow and S. M. McKelvie: A nitrogenbased model of plankton dynamics in the oceanic mixed layer, J. Mar.
Res., Vol. 48, pp. 591–639, 1990.
Gruber, N., H. Frenzel, S.C. Doney, P. Marchesiello, J.C. McWilliams,J.R. Moisan, J.J. Oram, G.-K. Plattner and K.D. Stolzenbach: Eddy-resolving simulation of plankton ecosystem dynamics in the California
Current System, Deep-Sea Res. I , Vol. 53, pp. 1483–1516, 2006.
Morel, A. and Berthon, J.F.: Surface pigments, algal biomass profiles, and potential production of the euphotic
layer: Relationships reinvestigated in view of remote-sensing applications, Limnol. Oceanogr., Vol. 34, pp. 15451562, 1989.
Willmott, C.J.: On the validation of models, Phys. Geogr.,Vol. 2, pp. 184-194, 1981.
(Received March 17, 2022)(Accepted July 21, 2022)

DEVELOPMENT OF AN EVALUATION SYSTEM FOR SUITABILITY OF
OFFSHORE WIND FARM SITES AROUND JAPAN BASED ON A REGIONALOCEAN CIRCULATION–BIOGEOCHEMICAL MODEL
Toshitake TAMAI, Yusuke UCHIYAMA, Xu ZHANG,
Eiji MASUNAGA, Kazuhiro TSUJI and Hidekatsu YAMAZAKI
Expansion of offshore wind power generation facilities (or offshore wind farm; OWF) has become a social
issue for achieving carbon neutrality. One of the problems unique to Japan is the necessity to obtain theconsent of fishing rights holders before constructing OWFs. The ultimate goal of the present study is todevelop a comprehensive consensus-building system by quantitatively evaluating the nutrient supply to theocean surface layer by enhanced vertical mixing associated with OWF facilities, which is often referred toas a "fishing bank effect" to provide gathering places for fishes. This manuscript particularly introduces acore component of the proposing system based on ROMS-NPZD regional circulation-biogeochemical
model that further drives an embedded high-resolution model around the OWFs. Some results on the spatial
evaluation method for each coastal and oceanic area to seek for suitable OWF sites are reported.

I_1014