北海道の雪氷 No.39（2020）

Annual Report on Snow and Ice Studies in Hokkaido

AMSR2 データを用いた北極海一年氷の海氷厚推定アルゴリズムの改良
Improving the algorithm to estimate sea ice thickness in Arctic by AMSR2
渡辺 由梨加 1，舘山 一孝 2, 泉山 耕 3, 佐藤 和敏 2
Yurika Watanabe1, Kazutaka Tateyama 2, Koh Izumiyama3, Kazutoshi Sato2
Corresponding author: m2052200226@std.kitami-it.ac.jp (Y.Watanabe)
Satellite data shows a decrease trend in sea ice extent and the increase in fraction of first-year ice in the Arctic Ocean.
Therefore, sea ice observation by satellite remote sensing is becoming more important for monitoring the global climatechange and using the Arctic Ocean routes. This study focuses on thickness of the first-year ice, and improve the accuracyof the algorithm to estimate sea ice thickness in Arctic by AMSR2.
１．はじめに北極海の海氷は，地球の平均気温が上昇していることに対し減少傾向にあり，特に海氷が融解する夏季の減少傾向は著しい．気象庁 1)によると海氷域面積の年最小値は，1979 年から 2019 年の間に 1 年当たり 8.9 万 km2 の速度で減少しており，これは北海道の面積（8.3 万 km2）に匹敵する．
このような北極海の変化により，北極海には以下に記述する 2 つのことが期待されるようになった．
第一に北極海航路の実利用である．海氷の少なくなる夏季には，北極海の海氷が融解した海域を船が航行できる．北川ら 2)によると，南周り航路に対し北極海航路は 60%程度の航程であるため，経済効果が大きいとされている．
第二に北極海で石油事業の発達が見込まれることである．アメリカ地質調査所（USGS：United StatesGeological Survey）3)によると，北極には約 900 億バレルの石油，約 50 兆 m3 の天然ガス，また約 440 億バレルの液体天然ガスが未発掘の状態で存在すると推定されている．
北極海航行にせよ石油採掘にせよ，これらの事業を安全に執り行うには海氷を避けることが必須である．そのためには海氷の場所，種類，形状を常に観測する必要がある．また，海氷は気候の影響を敏感に受けるため，それらの海氷に関する情報は，気候を評価する際のパラメータとして使用することができる．北極海は極めて広大であり，環境も過酷であるため，常に現場で観測を行うことは困難である．そのため，人工衛星に搭載したマイクロ波リモートセンシングによる観測が最適であると考える．
Krishfiled ら 4)は宇宙航空研究開発機構（JAXA：
北見工業大学大学院 社会環境工学専攻
北見工業大学 地球環境工学科
3
北海道大学 北極域研究センター
12Japan Aerospace Exploration Agency）の水循環変動観
測 衛 星 「 し ず く 」（ GCOM-W1 ： Global Change
Observation Mission 1st-Water "Shizuku"）搭載の高性能マイクロ波放射計（AMSR-E：Advanced Microwave
Scanning Radiometer for EOS）で観測された 6GHz から 89GHz までの 7 つの周波数の輝度温度データを使用して海氷厚推定式を作成した．以下のような 2つの輝度温度比（PR：Polarization Ratio）（GR：GradientRatio）を用いて海氷厚を算出する．1 つ目は 36GHzの垂直偏波𝑇𝐵 ，36GHz 水平偏波𝑇𝐵 による比𝑃𝑅 である．
𝑃𝑅𝑇𝐵

𝑇𝐵

𝑇𝐵/ 𝑇𝐵(1)

2 つ目は 36GHz の垂直偏波𝑇𝐵 と 18GHz の垂直偏波𝑇𝐵 による比𝐺𝑅である．
𝐺𝑅𝑇𝐵

𝑇𝐵/ 𝑇𝐵𝑇𝐵(2)

0.025を閾値
これらの比を用いて，𝐺𝑅とし，それぞれ以下の式で推定海氷ドラフト𝐷を算出する．一年氷（𝐺𝑅0.025）のとき
𝐷m

2.34 ∙ exp

.
.

0.085

(3)

多年氷（𝐺𝑅0.025）のとき
𝐷m
0.162 0.244 ∙ exp 20.79 ∙ 𝐺𝑅(4)5)
舘山ら はこれらのアルゴリズムの精度を高めるため，AMSR-E の後継機である AMSR-2（Advanced
Microwave Scanning Radiometer 2）を使用し，𝐷Graduate School of Civil Engineering, Kitami Institute of Technology
School of Earth, Energy and Environmental Engineering, Kitami Institute of Technology
Arctic Research Center, Hokkaido University－ 23 －Copyright©2020 公益社団法人日本雪氷学会
The Japanese Society of Snow and Ice

北海道の雪氷

No.39（2020）

Annual Report on Snow and Ice Studies in Hokkaido

を全氷厚に換算する式及び海氷厚を表面温度によって補正する式を開発した．
泉山 6)は以上のアルゴリズムの精度を検証するため，Beaufort Gyre Exploration Project（BGEP）による，ボーフォート海の 4 地点において Upward LookingSonar（ULS）により計測された海氷ドラフトの日平均データと，先行研究で完成したアルゴリズムのドラフトデータにおけるデータ数のヒストグラムを作成した（図 1）．横軸が海氷のドラフトの大きさ，縦軸がそのドラフトに対応するデータの数である．
ヒストグラムでは，約 60cm 付近の推定ドラフトのデータ数が ULS による計測ドラフトに対して大幅に欠けており，反対に約 70～90cm の部分にデータが集中している．このことから約 1m 以下の海氷，すなわち 1 年氷の推定精度にばらつきがあると言える．

図 1 BGEP 日平均ドラフト（上）と先行研究による推定ドラフト（下）の分布．泉山 7)により作成された.
本研究の目的は，推定精度が低く，気候変動の影響により割合が増加し重要性の増した一年氷に着目し，推定精度を向上させるための新たな推定式を提案することにある．
２．一年氷の海氷厚推定アルゴリズムの改良
２．１ 使用したデータ
輝度温度のデータは JAXA の AMSR2 により計測されたものを使用した．
AMSR2 の輝度温度データと比較し推定式を作成するため，定着氷上でドリル掘削により測定された氷厚データ（以下，「実測氷厚」という）を使用した．

これらの測定データには，山崎哲秀氏が 2012 年から 2019 年までの期間にグリーンランドで測定したデータと，カナダの Ice Thickness Program により2012 年から 2016 年の期間に測定されたデータを使用した．8 地域，10 地点での測定データが含まれている．
２．２ 海氷厚推定式の導出方法
AMSE2 には，フットプリントという周波数毎のビーム幅に応じた観測エリアがある．このフットプリントの中に陸地が入ると海氷厚の推定に影響するため，，AMSR2 の輝度温度データ取得場所を，少なくとも陸地から 30km は離して設定する必要がある．
山崎氏のデータやカナダの Ice Thickness Programの実測氷厚データは，いずれも陸地付近で海氷厚を計測している.そのため，可能な限り現地観測地点から離れないよう衛星データ取得箇所を陸地から30km 離した．このとき，陸地付近の氷厚と衛星データ取得地点の氷厚は同じであると仮定している．
また陸地に囲まれており 30km 離すことのできない地点のデータは，本研究では除外した（図 2）．

図 2 AMSR2 データ（緑）と実測氷厚データ（赤）を取得した地点．数字は各地点間の距離．数字の記載がないものは除外した地点．
はじめに，海氷厚推定アルゴリズムが海氷の密接度や融解度，種類による影響を受けないよう，データの選別を行う．海氷密接度は AMSR2 が提供している海氷密接度を使用し，100%未満のデータを排除した．融解度は田中ら 7)によるメルトポンド割合推定式を使用し，30%以上は融解とみなして排除した．また，海氷の種類は式 (2) を使用し，一年氷と多年氷に区別した後，多年氷のデータを排除した．
上記の手法によって選別された各周波数の輝度
温度データを実測氷厚と比較し，各周波数の輝度温－ 24 －Copyright©2020 公益社団法人日本雪氷学会
The Japanese Society of Snow and Ice

北海道の雪氷 No.39（2020）

Annual Report on Snow and Ice Studies in Hokkaido

度データを図 3 に示した．横軸が実測氷厚であり，縦軸が各周波数の輝度温度である．輝度温度は，赤色が水平偏波によるもの，青色が垂直偏波によるものである．10GHz 水平偏波の輝度温度が実測氷厚に対して最も相関が高く，ヒストグラムのばらつきも比較的少ないため，一年氷の海氷厚の変動に敏感であると言える．

𝑋𝑃𝑅𝑇𝐵

𝑇𝐵/ 𝑇𝐵

𝑇𝐵

輝度温度比𝑋𝑃𝑅と実測氷厚との関係を図 4に示した．横軸が実測氷厚，縦軸が度温度比
𝑋𝑃𝑅の値，データの色の違いは氷厚の計測地点の違いである．赤色のケンブリッジ・ベイのデータが散乱していることが分かる．ケンブリッジ・ベ
イの AMSR2 データ取得地点は氷厚計測地点から約92.81km と著しく距離が離れている．そのため，ケンブリッジ・ベイのデータは両者が離れたことによる海氷厚の変化が誤差を生んだと考えられるため，次の解析では除去した．その他の地点は約 36.38～62.50km 離れているが，すべて許容範囲内であるものとした．

図 4 実測氷厚と輝度温度比𝑋𝑃𝑅

𝑋𝑃𝑅

0.0136∙ exp

図 3 AMSR2 による各周波数の輝度温度と実測氷厚の関係．赤：水平偏波，青：垂直偏波
図中の数字は相関係数を表す．
この結果をもとに，先行研究で使用された輝度温度比を用いて海氷厚推定式を作成する．10GHz の水平偏波の輝度温度𝑇𝐵 と各周波数各偏波の輝度温度を組み合わせて輝度温度比を作成し，それぞれ実測氷厚と比較して決定係数 R2 を算出した．最も決定係数が高い輝度温度比は 6GHz 垂直偏波の輝度温度𝑇𝐵 との比 XPR（Cross Polarization Ratio）であるため，この輝度温度比𝑋𝑃𝑅をもとに推定式を作成する．

(5)の関係図

0.04860.00381 ∙ 𝑥

図 5 実測氷厚と輝度温度比𝑋𝑃𝑅
回帰曲線の関係図とケンブリッジ・ベイのデータを除去し，データに回－ 25 －Copyright©2020 公益社団法人日本雪氷学会
The Japanese Society of Snow and Ice

北海道の雪氷

No.39（2020）

Annual Report on Snow and Ice Studies in Hokkaido

帰直線を描いたグラフを図 5 に示した．横軸が実測氷厚，縦軸が輝度温度比𝑋𝑃𝑅の値である．この回帰曲線から海氷厚推定式を算出する．
３．海氷厚推定結果
以上の手法をもとに，実測氷厚と輝度温度比𝑋𝑃𝑅の回帰曲線を使用して作成した一年氷の推定式が以下である．の海氷厚𝐻

𝐻m.

∙ 𝑙𝑛.
.

(6)

𝐻の推定精度を検証するため，図 6 を作成した．横軸が実測氷厚，縦軸が𝐻である．データの色の違いは，計測地点の違いである．平均二乗誤差は 1.77 m である．

つきが大きいため，距離が離れて海氷厚が変化していることがばらつきの主な原因であると考えられる．しかし，両者の距離がさほど遠くないホール・ビーチのデータが著しくばらついていることから，距離だけでなく地点ごとの特性も影響している可能性が考えられる．
【謝辞】
本研究は山崎哲秀氏によって取得された現地観測データと JAXA から提供された AMSR2 データを使用した．
【参考文献】
1) 気象庁，海氷域面積の長期変化傾向（全球）URL:http://www.data.jma.go.jp/gmd/kaiyou/shinda

2)

3)

4)
図 6 実測氷厚と推定氷厚𝐻の関係図

４．考察
本研究で作成した推定式は，図 6 によると海氷厚が 1m 以下である一年氷については概ね測定値と推定値が同じである．しかし，特定の地点（ホール・ビーチ）で著しく過大評価している．このような過大評価する原因として，本研究で採用した回帰曲線の線形が原因の１つと考えられる．海氷厚の値が小さい範囲の回帰曲線は傾きが大きく，回帰曲線の変化に応じて適切な推定氷厚を示すことができる．しかし，海氷厚の値が大きくなるにつれ，回帰曲線は傾きが小さくなり，わずか輝度温度の変化で急激に推定氷厚が増加してしまう．
また，図６から地点ごとにばらつき度合いが異なっていることが示された．実測地点と AMSR2 デー
タ取得地点が比較的近いアラート 1，2 ほどばらつ
きが少なく，両者の距離が遠いその他の地点はばら5)

6)

7)－ 26 －

n/a_1/series_global/series_global.html（2020 年 6
月 18 日閲覧）北川弘光，小野延雄，山口一，泉山耕，亀崎一彦，2000：「北極海航路 ―東アジアとヨーロッパを結ぶ最短の海の道―」，財団法人シップ・アンド・オーシャン財団，東京，扇興社，1-3.
United States Geological Survey，2008：CircumArctic Resource Appraisal ， Estimates ofUndiscovered Oil and Gas North of the Arctic Circle
URL:https://pubs.usgs.gov/fs/2008/3049/（2020 年
3 月 23 日閲覧）
R. A. Krishfield，A. Proshutinsky，K. Tateyama,，
W. J. Williams，E. C. Carmack，F. A. McLaughlin，
M.-L. Timmermans ， 2014 ： Deterioration ofperennial sea ice in the Beaufort Gyre from 2003 to2012 and its impact on the oceanic freshwater cycle，Journal of Geophysical Research，101，1297-1303
Kazutaka TATEYAMA ， Jun INOUE ， Seita
HOSHINO，Shota SASAKI，Yasuhiro TANAKA，
2018：Development of a new algorithm to estimate
Arctic sea-ice thickness based on Advanced
Microwave Scanning Radiometer 2 data，Okhotsk
Sea and Polar Oceans Research，2，16-17
泉山耕，2019：KIT アルゴリズムによる海氷厚
推定 ―現地観測結果との比較と高度化への考察―，私信，46-51
Tanaka, Y., K. Tateyama, T. Kameda, and J. K.
Hutchings，2016：Estimation of melt pond fraction
over high-concentration Arctic sea ice using
AMSR-E passive microwave data.，
Oceans, 121（9），7065Copyright©2020 公益社団法人日本雪氷学会
The Japanese Society of Snow and Ice