学際大規模情報基盤共同利用・共同研究拠点 平成 23 年度共同研究 最終報告書 2012 年 5 月

11-NA15

環オホーツク圏を中心とした大気・海洋シミュレーション
中村知裕（北海道大学）
概要 環オホーツク圏では、オホーツク海起源の海洋熱塩循環・潮汐混合・栄養物質（特に鉄）循環に伴い、世界でも最大規模の基礎生産が生じ、豊富な水産資源や炭素循環に影響している。また、夏季の下層雲やオホーツク海高気圧の形成等を介して北日本の気候に影響を与えている。最近ではこれらに長期変動が生じていることも明らかになってきた。本研究では、こうした現象の理解と数値モデルによる再現性の向上に向け、環オ
ホーツク圏の大気・海洋シミュレーションと使用する数値モデルの高速化チューニングを行う。本年度は、先ず(I)計算機科学分野の課題として、大気モデルを水平２方向にMPI 並列処理を行い SR16000 での性能テストを行った。次に(II)環オホーツク圏大気・海洋シミュレーションとして、(a) 海洋内鉄循環シミュレーションと各種パラメタへの感度実験、(b) 潮汐による鉛直混合過程の数値実験、(c) オホーツク海高気圧の強い年と弱い年の夏季下層雲シミュレーションとその比較、を行った。

1. 研究の目的と意義での潮汐鉛直混合により、オホーツク海起源の熱研究の目的

塩循環（密度差に由来し、鉛直方向の輸送で特徴環オホーツク圏（オホーツク海とその周辺地域） 付けられる海洋循環）が駆動される [Nakamura etは、季節変動や気候変動に顕著な特徴を持ち、日al., 2006]。熱塩循環に伴い、大気に接していた本とりわけ北日本の気候に大きな影響を与えてい海水が海洋中層に潜り込むので、大気中の酸素、る。例えば夏季には、オホーツク海周辺では下層温室効果気体（二酸化炭素など）、フロンといった雲（高度の低い雲）や霧が頻繁に形成される。下様々な気体が北太平洋中層に取り込まれ、それら層雲は、日射を遮り赤外線を放射することで大気物質の循環に影響する [Uchimoto et al., 2011;－海洋間あるいは大気－陸面間の熱・水フラック

2009]。

スに大きく影響する [中村・古関, 2010]。その結

環オホーツク圏はまた、世界最大規模の基礎生

果、下層雲の下では気温が低下し、これが顕著に産（海洋植物プランクトン増殖）で知られている。

なると農作物の不作の原因となる。一方冬季には、 高い基礎生産は、食物連鎖を通し豊富な水産資源
シベリア東部は北半球の寒極（気温の最も低い地の基礎となるとともに、二酸化炭素の海洋内貯留域）となり、そこから吹き出す寒気によりオホーなど炭素循環にも重要な役割を果たす。この高いツク海では大規模な海氷生成が起こる。

基礎生産を支える上で、上述の熱塩循環によりアオホーツク海における海氷生成は北海道沿岸のムール川から運ばれて来た「鉄」が重要であるこ

気候や産業に影響するだけでなく、北太平洋ほぼとが最近の研究から分かってきた [Nishioka et全域の海洋中層循環に影響を与えている。海氷生

al., 2007]。すなわち、鉄は２価と３価のイオン成の際には、海水が結氷点まで冷却される上に、を持つことから光合成等における電子伝達に使わ海水が凍る際に不純物として塩が排出されるため、 れているが、海水には極めて溶けにくいため多く高密度の海水が生成される。生成された高密度水の海域において基礎生産を律速している。例外的はオホーツク海そして北太平洋の中層に広がる。に環オホーツク圏では、熱塩循環に伴う鉄供給のまた、千島列島域のように海底地形上を強い潮おかげで鉄律速がかかりづらく、このことが水産流が流れる海域では、内部重力波の生成と砕波に資源の豊さに繋がっている。

より激しい鉛直混合が生じる[Nakamura et al.,

環オホーツク圏におけるこれらの現象は、顕著

2000; 2010]。これらの高密度水生成と千島列島域な数年から数十年規模の変動を持つのに加えて、
1

学際大規模情報基盤共同利用・共同研究拠点 平成 23 年度共同研究 最終報告書 2012 年 5 月

シベリアを中心に顕著な温暖化が進行している。

にとって重要な地域である。しかしながら、大半温暖化の特に著しい地域では、平均して 10 年間にがロシア領ないしロシア経済水域に含まれている2℃のペースで冬季の気温が上昇している。温暖化ため、大気・海洋・陸面全てにおいて公開されてに伴うようにオホーツク海の海氷面積も減少傾向いる観測データが限られている。また、基礎生産にあることから、上述のオホーツク海起源の熱塩に重要な鉄は、海水中に比べ観測船内の方が遙か循環が弱まっていく可能性が危惧されているに量が多いため観測には特殊な技術を要するので[Nakanowatari et al., 2007; Matsuda et al.,

断片的知識しか得られていない。そのため環オホ

2009]。

ーツク圏の気候・環境研究には、数値シミュレー

このように、環オホーツク圏では科学的に興味ションとそれに基づくメカニズム解明、そして解深くかつ社会的に重要な気候・環境の形成および明されたメカニズムに基づく数値シミュレーショそれらの長期変動が生じており、これらのより良ンの改良が欠かせない。

い理解と数値シミュレーションによる再現・予測が求められている。そこで代表者らのグループで

2. 当拠点公募型共同研究として実施した意義

は、環オホーツク圏の気候および環境変動のより(1) 共同研究を実施した大学名と研究体制

良い理解のために、環オホーツク圏を対象とした

北海道大学

海洋および大気の数値シミュレーションを行って(2) 共同研究分野いる。本研究ではその一環として、環オホーツク

超大規数値計算系応用分野

圏の気候・環境に重要な３つの現象―(a) オホー(3) 当公募型共同研究ならではという事項など

ツク海を起源とする栄養物質循環と熱塩循環、(b)本共同研究により、計算科学の専門家と共同研

潮汐による鉛直混合過程、(c) オホーツク海周辺究を行えたことが本事項の第一である。このことの下層雲―について数値シミュレーション研究をは、数値モデルの高速化および、北海道大学情報行う。加えて、これらシミュレーションの効率化

基盤センターの大型計算機システム更新後の新シと大規模化に向けて、大型計算機更新後の新シスステムへの移行、新システムでの高速化チューニテムへの移植と新システムでの高速化チューニン

ングで特に顕著であった。

グを検討する。

また、北海道大学情報基盤センターの大型計算機システムを使用させて戴くことで、数値シミュ研究の意義レーションとその結果の解析を円滑に且つ高速に現在気候・現在環境のシミュレーションとそれ実施できるのも当共同研究の利点である。

らの形成・変動メカニズム解明は、地球温暖化に加えて、JHPCN シンポジウムにおいて他分野の伴う環オホーツク圏の変化の理解および信頼性の計算科学に関する専門家と交流し、その研究発表高い予測の必要条件である。中でも、環オホーツを拝聴できることは、当共同研究以外ではなかな

ク圏では近年表層の栄養物質そして基礎生産が減か得られない機会であった。

少傾向にあり、温暖化がその主な原因ではないかと示唆されている。もしこの示唆が正しければ、

3. 研究成果の詳細と当初計画の達成状況

温暖化が進むと基礎生産がさらに減少し、ひいて(1) 研究成果の詳細については水産資源・炭素循環も大きく影響を受けること(I)計算機科学分野での検討課題から、先行きが懸念されている。本研究の成果は、本課題では、昨年度に SMP 並列をチューニング

こうした変動の要因解明と予測の重要な基礎となした大気モデルの MPI 並列化、および北海道大学る。

情報基盤センター大型計算機システム更新後の新対象とする環オホーツク圏は上述のように日本システムでの高速化チューニングを行った。
2

学際大規模情報基盤共同利用・共同研究拠点 平成 23 年度共同研究 最終報告書 2012 年 5 月

本研究で使用している領域大気モデル(IPRC 領

本年度は、先ず、このモデルを図１b のように域気候モデル)は、図 1a に示すように南北方向（以東西方向（以下 j 方向）にも領域分割した。なお、

下 i 方向）に MPI 並列処理されている。SR11000

モデル・コードでは、j が最内ループ、次が i、最で計算する際は、分割した領域を SMP 並列とする外が k（鉛直方向）である。鉛直方向は反復計算

ハイブリッド並列を行っており、昨年度は、そのや陰解法に伴う通信が多いので、鉛直方向には領SMP 並列の高速化チューニングを行った。

域分割しない。

図１：(a)元々の領域分割。水平一方向で分割し、各領域内で SMP 並列する。(b)水平２方向に領域分割した。

図２：水平２方向に領域分割したことに伴う袖、及びそのデータ通信先。

領域分割に伴い、「袖」などと呼ばれるダミー・チューニングを行う。

グリッドを作成し、その通信を行う。今回の場合

性能プロファイルを調べた所、極端に性能が低の設定を図 2 に示す。モデル中の水平差分精度は

くチューニング幅の大きいループはないことが確最高で４次なので、袖は２グリッドですむ。袖の認された。そこで、測定条件をいろいろ変えて性作成に伴い do ループの範囲も変更した。

能を比較した。その結果を表 1 に示す。

次に、上記のように水平両方向に領域分割した

水平２方向で領域分割した効果を、先ず、MPIモデルを、北海道大学情報基盤センターの新しい並列と SMP 並列を併用した hybrid 並列処理で測定大型計算機システム(SR16000M1)に移植し、高速化した所、２方向で分割した場合の方が遅くなって3

学際大規模情報基盤共同利用・共同研究拠点 平成 23 年度共同研究 最終報告書 2012 年 5 月

しまった（表 1 の 1,2 段目を比較）。これは、SMP

また、このテストケースでは問題規模を小さくし並列高速化に有効であった「物理スキームのよりているが、分割数を 128（４倍）にしても約 3 倍上位ルーチンでの SMP 並列化」を妨げるルーチンの向上が見られた（表１の最下段とその上を比較）。

（積雲パラメタリゼーション内の通信）があるた

最後に、SR16000 の SMT 機能（一つの物理コアめだと判明した。

で複数のスレッド及びプロセスを同時に実行出来

そこで次に、MPI 並列のみを用いて比較した。

る機能）についてもテストを行った。本モデルでi,j 両方向に 8 分割すると、i 方向のみに 64 分割

は比較的大量のメモリを使用するため、残念ながした場合に比べ、計算時間が約 0.7 倍に減少した

ら、SMT 機能をオフにした方が性能が向上してい（表 1 の 4,5 段目を比較）。これは、hybrid の場る（表 1 の 1,3 段目、4,6 段目をそれぞれ比較）。

合よりも向上している（表 1 の 1,4 段目を比較）。

表１：水平両方向に領域分割した後の大気モデル性能測定結果。

以上のテストを元に、高分解シミュレーション(II) 環オホーツク圏大気海洋シミュレーションを開始した。図 3 は、実際に用いる高分解能モデ(a) オホーツク海を起源とする栄養物質循環と熱ルの設定での測定結果である。大気モデルを１日塩循環間積分して経過した実時間を測定し、プロセス数

本課題では、環オホーツク圏の栄養物質循環を（分割数）160 の場合を基準にした相対速度で示特徴付ける鉄に注目してシミュレーションを行う。

している。プロセス数 1280 までは向上が見られる

そのため、熱塩循環を概ね良好に再現している海が、並列化効率は徐々に落ち、1920 ではほとんど

洋大循環モデル[Uchimoto et al., 2011]に、鉄化向上しなかった。これは先に出て来た積雲パラメ学モデル[Parekh et al., 2005]を組み込む。本年

タリゼーション内の通信が原因であった。解決に度は、気候学的季節変動場について鉄化学モデル

は、パラメタリゼーションの計算方法またはパラのパラメタ・テスト等を行い、鉄・リン分布の再メタリゼーションを見直す必要があり、これは今現性向上を行った。

後の課題とする。

鉄化学モデルを組み込んだシミュレーションの結果得られた、溶存鉄濃度の中層における水平分布（図 4a）とオホーツク海西岸沿いの鉛直分布（図4b）、および海面でのリン酸(PO4)濃度（図 4c）を示す。これらは全球モデル用に調整されたパラメタをそのまま用いているが、定性的な特徴はかなり良好に再現されている。例えば図 4a に見られるように、鉄濃度が中層ではオホーツク海の方が太図 3：高分解能シミュレーションの設定を用いた測定結果。

4学際大規模情報基盤共同利用・共同研究拠点 平成 23 年度共同研究 最終報告書 2012 年 5 月

平洋より高く、オホーツク海の中でも陸棚高密度及び中層で高い。逆にリン酸はオホーツク海より水の経路であるオホーツク海北部の陸棚からオホ太平洋で高く、中でもアリューシャン列島域、亜ーツク海西岸にかけて高い。

寒帯域、および千島列島域で高い（図 4c）。これ
らの特徴は観測から指摘されている定性的特徴と(a) 溶存鉄濃度：中層一致している。加えて、オホーツク海内の鉄濃度の値も数少ない観測値から推測されている現実的な範囲に収まっていた。

(b) 溶存鉄濃度：

図 5：鉄化学モデルとその設定の模式図。

西岸沿いさらに再現性を向上するため、および鉄循環への理解を深めるため、各種パラメタに対する感度実験を行った。今回用いた鉄化学モデルの模式図(c)リン酸(PO4)濃度：海面（図 5）に示すように、重要な鉄化学過程として、海洋中の鉄のソースである(1)風送塵の沈着と溶解および(2)堆積物からの溶解または巻き上げ、シンクとなる(3)凝集・捕集とそれに関わる錯形成、そして内部の循環に関する(4)生物への取込および(5)その粒子としての沈降と粒子の再無機化がある。
これらの過程には、実は、現在でも良く分かっていない点あるいは専門家の間で意見の分かれる

図 4：鉄循環シミュレーション結果。気候学的季節変動場

点が幾つもある。中でも甚だしいのが、鉄のソーにおける夏季。(a)溶存鉄濃度(nM)の水平分布。中層の等スに関わる風送塵の溶解度と堆積物からの鉄供給ポテンシャル密度（26.8）面上のもの。(b) 溶存鉄濃度率である。溶解度は 0.1%～10%と研究者により大の鉛直分布。オホーツク海北西陸棚域（東経 140 度、北緯きく異なる値が提案されている。堆積物からの鉄57 度付近）からオホーツク海西岸流域を通り、千島列島中供給は本来、海底付近の流速や海底の堆積物の量央部に至る経路に沿ったもの。(c) 海面におけるリン酸濃と状態などに応じた海域依存性があるはずだが、度(M)。

観測がカリフォルニア沖で１度実施されたものし鉛直断面で見ると図 4b のように、鉄濃度は陸棚域かないため、そのときの観測値が全球で用いられ5

学際大規模情報基盤共同利用・共同研究拠点 平成 23 年度共同研究 最終報告書 2012 年 5 月ている。加えて、生物への取込も、本来は種によさいケースから順に左から並んでいる。一見してり取込率が異なるので種組成への依存性があり、分かるように、海面に比べると中層は溶解度へのしたがって海域によって異なる。そこで、これら依存性が非常に小さい。鉄の観測データは最適なのパラメタについてテストを行った。

溶解度を推定するには少なすぎるので、リン酸濃先ず、沈着した風送塵の溶解度について感度実度の海面分布から判断すると溶解度は 1%の時が最験を行った。その結果を図 6 に示す。上段は海面、も観測データに近かった（図は省略）。

下段は中層（26.8）での鉄濃度で、溶解度の小図 6：風送塵の溶解度に対する感度実験結果。海面（上段）と中層 26.8（下段）における鉄濃度。左から溶解度が 0,1,5,10%
の場合。

次に、堆積物からの鉄供給率に対する感度実験

さらに、どの海域の堆積物が重要かを調べるたを行った。堆積物からの供給をなくすと（図 7）、め、オホーツク海内の陸棚を図 8 の a-d のように中層の鉄濃度は劇的に減少し、「オホーツク海の方４つに分け、それぞれの海域のみで堆積物からのが太平洋より鉄濃度が高く特に中層で鉄濃度が高鉄供給を与えた。その結果が図 9 で、a-d はそれい」といった基本的特徴が消えた。すなわち、オぞれ図 8 の a-d に対応する。基本的にソース域とホーツク海周辺の中層鉄循環シミュレーションにそれより下流で鉄濃度が高いが、下流に行くにつ堆積物からの供給は欠かせない。

れて影響は弱くなっていく。
図 8 の各海域での堆積物からの供給率を変化させ、オホーツク海西部中層の鉄濃度に対する影響調べると、(1)供給域から、距離が伸びると濃度が減少していくこと、(2)どこが供給域でも供給率が0 から 1  mol m-2d-1 くらいまでは劇的に鉄濃度が高まることが分かった（図は省略）。これは、堆積物からの鉄供給が少ない場合は、リガンドが余っているので、溶存鉄が増加しやすいのに対し、鉄供給が増えて海中の鉄濃度が高くなるとリガン

図 7：堆積物からの鉄供給がない場合の、中層（26.8）

ドが少なくなり凝集・捕集されやすくなるためと鉄濃度。

考えられる。また、(3)海水中の鉄濃度が一定以上であれば、堆積物からの鉄供給率への応答は線形
6(a)

(b)

学際大規模情報基盤共同利用・共同研究拠点 平成 23 年度共同研究 最終報告書 2012 年 5 月に近くなった。このことを利用すると、各海域にしている潮汐による鉛直混合過程について数値的与えるべき供給率の粗い推定が得られる。に調べた。代表者らの海洋大循環モデルではその(a)(b)効果が簡単なパラメタ化により考慮されているが、超高分解能非静水圧海洋モデルによる鉛直混合過程シミュレーションを行い、パラメタ化の改良を図る。本年度は、昨年度に行った数値シミュレー
ションの詳細な解析、および理論化を目指しより(c)(d)理想化された条件下での数値実験を行った。現在、査読付国際誌に投稿準備中であり、詳細の記述は投稿後または受理後に行いたい。
(c) オホーツク海周辺の下層雲
昨年度に高速化チューニングした領域大気気候

図 8：堆積物からの鉄供給を与える海域。

(a)

モデルを用いて、夏季の下層雲のシミュレーショ

(b)ンを行い、下層雲の形成・維持機構を調べた。昨年度は、オホーツク海高気圧の発達した 2003 年 7月に注目した。本年度は、オホーツク海高気圧の発達しなかった年についてもシミュレーションを(c)行い、比較することで下層雲形成・維持のメカニ(d)ズムと再現における鍵を調べた。
例として、オホーツク海高気圧が発達した別の年 1998 年とオホーツク海高気圧が発達せず、海面気圧が平年値より有意に低かった 2004 年について 7 月月平均の下層雲分布を図 10 に示す。

図 9：堆積物からの鉄供給を図 8 中に赤で示した海域のみ

下層雲はオホーツク海高気圧が顕著だったで与えた場合の、中層（26.8）鉄濃度。

2003 年に最も多いが、2004 年にも形成されてい最後に、生物への取込過程に関するパラメタのる。1998 年はむしろ比較的少ない。分布も年によテストを行った。取込の量はリン酸濃度(PO4)・鉄って大きく異なる。それでも詳細に見ると、どの濃度(Fe)・光量(I)を用いて、次式でパラメタ化さ年も定性的には、陸域や太平洋などオホーツク海れている。この式中のαと KFe についてテストした。 の外から風の吹き込む海域で下層雲が少なく、気(a)(b)PO 4FeIPO 4  K PO4 Fe  K Fe I  K I

塊がオホーツク海上に長く止まる海域で下層雲が基本的に、海面のリン酸濃度はαが大きいと減燥した気塊が、オホーツク海上で低い海面水温に少し、KFe が大きいと増加する（図は省略）。そこ伴う顕熱フラックスで冷却されて、相対湿度が高で、海面リン酸濃度が気候学的分布に近くなるよまり雲ができて気温が下降し、放射冷却と海面熱うにこれらパラメタを調整した。

フラックスがバランスして下層雲が形成された状 多い。つまり、オホーツク海外の比較的高温で乾態で安定に近づくという過程は多くの年で共通し(b) 潮汐による鉛直混合過程ている。より定量的な解析は今後の課題である。

本課題では、熱・塩・栄養物質の鉛直輸送を介して熱塩循環および物質循環に重要な役割を果た
7

学際大規模情報基盤共同利用・共同研究拠点 平成 23 年度共同研究 最終報告書 2012 年 5 月

件下での数値実験を終了した。成果は国際誌に投稿準備中である。
(c)オホーツク海下層雲：計画通りオホーツク海高気圧の発達した年とそうでない年について夏季下層雲シミュレーションを行った。当初計画に加えて、(I)で MPI 並列化した大気モデルを用いた高分解能シミュレーションを開始した。
図 10：下層雲の4. 今後の展望

分 布 。 (a) 1998

(I)計算機科学分野については、今年度の成果で年、(b) 2004 年、大気モデルの高速化が一段落したので、次は海洋(c) 2003 年の７モデルの高速化チューニングについて取りかかり月月平均。(a)(b)たい。の矢印は下層風。

(II) 環オホーツク圏大気海洋シミュレーショ
ンについては、今年度の成果を踏まえて次は、(a)(2) 当初計画の達成状況について鉄循環の経年変動シミュレーション、(b)鉛直混合

課題申請時の計画は全て達成された。それに加に寄与する３次元的過程の数値実験、(c)下層雲のえて来年度に向けた準備（大気シミュレーション）

高解像シミュレーションを実施したいと考えていを開始している。以下、具体的項目毎に達成状況る。

を記す（詳細は 3 節参照）。

5. 研究成果リスト(I)計算機科学分野の課題：
本年度に予定していた水平双方向領域分割によ(1) 学術論文（投稿中のものは「投稿中」と明記）

る MPI 並列処理及び新スパコンへの移植と高速化

Koseki, S., T. Nakamura, H. Mitsudera, Y. Wang

を全て終了した。また、当初計画に加えて、新ス(2012). Modeling low-level clouds over the Okhotsk

パコンで数 10～100 ノードを用いる大規模計算のSea in summer: Cloud formation and its effects onテストを行った。なお、今回判明した SMP および

the Okhotsk high. Journal of Geophyscal Research,

さらなる分割化の妨げになっているルーチンは、

117, D05208, doi: 10.1029/2011JD016462.

コードの変更では解決出来ない、物理スキームにUchimoto, K., T. Nakamura, J. Nishioka, H. Mitsudera,関係するものなので今後の課題とする。それ以外

M. Yamamoto-Kawai, K. Misumi, D. Tsumuneの点については大気モデルの高速化チューニング

(2011). Simulations of chlorofluorocarbons in andはほぼ煮詰まったので、来年度以降は鉛直混合過

around the Sea of Okhotsk: Effects of tidal mixing

程の数値実験に用いているモデルなど、別のモデand brine rejection on the ventilation. J. Geophys.

ルの高速化も検討していきたい。

Res., 116, C02034, doi:10.1029/2010JC006487.

(II)環オホーツク圏大気海洋シミュレーション：

中村知裕, 古関俊也, 三寺史夫: 2012. オホーツク(a)鉄循環シミュレーション：計画していた気候

海における大気海洋相互作用：夏季の下層雲－値場の再現を実施した。またパラメタ・テストに海面水温フィードバック. 沿岸海洋研究（受理）基づく最適パラメタの推定を行い、再現性を向上

Abe, S. and T. Nakamura. Processes of Breaking ofした。得られた成果は、2012 年 3 月の国内学会でLarge-Amplitude Unsteady Lee Waves Leading to

発表済みである。

Turbulence. J. Geophys. Res.（投稿中）(b)潮汐による鉛直混合過程：計画通り理想的条

Nakamura, T., J. P. Matthews, T. Awaji, and H.
8

学際大規模情報基盤共同利用・共同研究拠点 平成 23 年度共同研究 最終報告書 2012 年 5 月

Mitsudera. Submeso-scale eddies near the Kuril(4) 国内会議発表

Straits: Asymmetric generation of clockwise and内本圭亮, 中村知裕, 西岡純, 三寺史夫, 三角和

counterclockwise eddies by barotropic tidal flow. J.

弘, 津旨大輔. オホーツク海のどこの海底に鉄Geophys. Res.（投稿中）

ソースがあれば高密度陸棚水の高濃度の鉄が実

中村知裕, 磯田豊, 三寺史夫, 高木省吾, 長澤真現され得るか? 2012 年度 日本海洋学会 春季大樹：2011. アムチトカ海峡で観測された大振幅会, 2012 年 3 月 26-30 日, 筑波大学, 茨城県つ内部波の砕波とグローバル評価. 月刊海洋, 通くば市.

巻 494 号, Vol.43, No.12, 699-704.

松田淳二, 笹島雄一郎, 三寺史夫, 中村知裕, 羽内本圭亮, 中村知裕, 西岡純, 三寺史夫, 川合美

角博康. 高解像度太平洋・オホーツク海モデル

千代, 三角和宏, 津旨大輔 (2011): オホーツの開発. 2012 年度 日本海洋学会 春季大会,ク海物質循環モデリング. 月刊海洋, 通巻 493

2012 年 3 月 26-30 日, 筑波大学, 茨城県つくば号, Vol.43, No.11, 682-687.

市.

(2) 国際会議プロシーディングス

内本圭亮, 中村知裕, 西岡純, 三寺史夫, 三角和なし

宏, 津旨大輔. オホーツク海における鉄移送モ(3) 国際会議発表

デリングの試み. 2011 年度 日本海洋学会 秋季

Nakamura, T, J. Nishioka, T. Ono, H. Mitsudera.

大会, 2011 年 9 月 26-30 日, 九州大学 筑紫キャンパス, 福岡県春日市.

Winter mixed layer development by subtropicalwater intrusion over subarctic water and associated松田淳二, 笹島雄一郎, 三寺史夫, 中村知裕, 羽iron supply. 2012 Ocean Sciences Meeting, Feb 20 -角博康. 高解像度太平洋オホーツク海海洋-海24, 2012, Salt Lake City, USA.

氷モデルの構築. 2011 年度 日本海洋学会 秋季大会, 2011 年 9 月 26-30 日, 九州大学 筑紫キ

Abe, S. and T. Nakamura. Transition process frombreaking

large-amplitudeinternalwaves

ャンパス, 福岡県春日市.

to三寺史夫,中村知裕. オホーツク海における大気

turbulence. 2012 Ocean Sciences Meeting, Feb 20 -海洋相互作用. 2011 年度 日本海洋学会 秋季大24, 2012, Salt Lake City, USA.
Abe, S. and T. Nakamura. Diapycnal mixing
generatedbybreakingof会, シンポジウム B, 東アジア縁辺海における

tideinduced

大気海洋相互作用と海洋生態系への影響－沿岸

large-amplitude internal waves. IUGG2011,

海洋シンポジウム－ 2011 年 9 月 26-30 日, 九28 Jun-7 Jul, 2011, P03_28PP028, Melbourne,

州大学 筑紫キャンパス, 福岡県春日市.

Australia.
Uchimoto, K., Nakamura, T., Nishioka, J.,(5) その他（特許，プレス発表，著書等）

Mitsudera, H., Misumi, K. and Tsumune, D.

三寺史夫, 中村知裕: 数値モデルを用いた環オホ

Toward a simulation of iron circulation fromーツク地域の環境研究―将来予測へ向けて. 環

the Okhotsk Sea to the Pacific. ESSAS 2011オホーツク海地域の環境と経済」スラブ・ユー

OSM, May 22-26, 2011, Seattle, WA, USA.

ラシア叢書 11, 61-88，田畑伸一郎、江淵直人

松田淳二, 笹島雄一郎, 三寺史夫, 中村知裕, 羽編、北海道大学出版会.

角博. Modeling of thermo-haline circulation of theSea of Okhotsk and North Pacific. 第 27 回 北方圏国際シンポジウム「オホーツク海と流氷」
2012 年 2 月 20-24 日, 紋別市文化会館, 紋別市.
9