冰川冻土, 2020, 42(4): 1135-1144 doi: 10.7522/j.issn.1000-0240.2018.1146

冰冻圈与全球变化

基于FY-3 MWRI数据的北极海冰密集度反演研究

吴展开,1, 王星东,1,2, 王峰1

1.河南工业大学 信息科学与工程学院,河南 郑州 450001

2.中国科学院 遥感与数字地球研究所,北京 100094

Inversion of Arctic sea ice concentration based on FY-3 MWRI data

WU Zhankai,1, WANG Xingdong,1,2, WANG Feng1

1.College of Information Science and Engineering,Henan University of Technology,Zhengzhou 450001,China

2.Institute of Remote Sensing and Digital Earth,Chinese Academy of Sciences,Beijing 100094,China

通讯作者: 王星东, 副教授, 从事基于主被动微波的极地冰盖/海冰研究. E-mail: zkywxd@163.com.

编委: 武俊杰

收稿日期: 2018-09-13   修回日期: 2019-03-04  

基金资助: 冰冻圈科学国家重点实验室开放基金项目.  SKLCS-OP-2020-6
河南省教育厅自然科学研究项目.  19B420001
河南省重点研发与推广专项.  202102310334
中国科学院数字地球重点实验室开放基金项目.  2018LDE005

Received: 2018-09-13   Revised: 2019-03-04  

作者简介 About authors

吴展开(1993-),男,河南兰考人,2017年在河南工业大学获学士学位,现为河南工业大学在读硕士研究生,从事极地海冰研究.E-mail:spreadwu@126.com , E-mail:spreadwu@126.com

摘要

以ASI算法(ARTIST sea ice algorithm)为基础, 得到基于风云3C气象卫星(FY-3C)微波辐射计(MWRI)数据的纯水与纯冰系点值, 利用插值方法确定基于FY-3 MWRI数据的ASI海冰密集度计算公式, 采用大津法(Otsu算法)得到基于MWRI数据的天气滤波器阈值。以2016年1月数据为例, 对北极海冰密集度进行反演, 并与美国国家冰雪数据中心(NSIDC)以及德国不莱梅大学提供的海冰密集度产品进行对比验证。结果表明: 基于MWRI数据得到的1月平均海冰面积以及平均密集度均介于二者之间, 其中平均密集度与不莱梅产品更接近, 仅相差1.310%。与风云卫星空间分辨率为250 m的中分辨率光谱成像仪(MERSI)数据得到的结果进行对比, 发现二者的海冰外缘线基本一致, MERSI数据得到的海冰密集度以及海冰面积比MWRI数据得到的结果分别高出5.029%、 9.318%。因此, 应用该方法可有效推进MWRI数据反演北极海冰密集度, 进而监测北极海冰分布和变化。

关键词: FY-3 ; MWRI ; 海冰密集度 ; ASI算法 ; 北极

Abstract

Based on the ASI (ARTIST sea ice) algorithm, the pure water and pure ice point values of the FengYun-3C meteorological satellite (FY-3C) microwave radiometer (MWRI) data are obtained. The ASI sea ice concentration calculation formula based on FY-3 MWRI data is determined by interpolation method, and the weather filter threshold is obtained by Otsu method. Taking the January, 2016 data as an example, the inversion results of the Arctic sea ice concentration were obtained. It was compared with sea ice concentration products provided by National Snow and Ice Data Center (NSIDC), USA, and University Bremen, Germany. The results show that the average sea ice area and average concentration in January based on MWRI data are between the two products, with a discrepancy of average concentration of only 1.310% different with Bremen products. The results are compared with those obtained from MERSI (medium resolution spectral imaging) data with 250 m spatial resolution of FengYun satellite. The results show that the outer line of sea ice is basically the same. The sea ice concentration and area obtained from optical data are 5.029% and 9.318% higher than those from MWRI data, respectively. Therefore, this method can effectively promote the inversion of Arctic sea ice concentration by MWRI data, and can be effectively used to monitor Arctic sea ice distribution and changes.

Keywords: FY-3 ; MWRI ; sea ice concentration ; ASI algorithm ; Arctic

PDF (1913KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

吴展开, 王星东, 王峰. 基于FY-3 MWRI数据的北极海冰密集度反演研究[J]. 冰川冻土, 2020, 42(4): 1135-1144 doi:10.7522/j.issn.1000-0240.2018.1146

WU Zhankai, WANG Xingdong, WANG Feng. Inversion of Arctic sea ice concentration based on FY-3 MWRI data[J]. Journal of Glaciology and Geocryology, 2020, 42(4): 1135-1144 doi:10.7522/j.issn.1000-0240.2018.1146

0 引言

近年来北极海冰面积和体积均急剧减少1-2, 北极海冰的消融已成为地球变暖的重要证据之一, 是全球生态系统中重要的指示和衡量标志3。通过大量的反馈机制与气候系统的相互作用, 北极海冰对北极地区以及全球气候都有显著影响4-6。此外, 准确了解北极海冰的状态对天气预报以及北极相关作业(钻探和导航)也至关重要7。其中, 海冰密集度是研究海冰变化的重要参量, 它是指单位区域内海冰覆盖率的百分比8。探测北冰洋海冰密集度及分布状况, 可为开拓北冰洋航路和全球气候的研究提供可靠的基础数据和科学依据。

目前比较成熟的利用星载微波辐射计亮温数据来反演海冰密集度的方法主要有: Cavalieri等9利用微波辐射计(Special Sensor Microwave/Image, SSM/I)19.4 GHz垂直和水平方向极化, 以及37 GHz垂直方向极化的亮温数据, 提出NASA Team算法, 该算法可用于一年冰以及多年冰密集度的反演; Comiso10-11利用基本能量运输方程和海冰发射率在不同频段的不同特征, 提出Bootstrap算法来计算海冰密集度, 该算法主要利用19(18.7)GHz和37(36.5)GHz两个频段的亮温数据来反演, 最终得到空间分辨率为25 km的海冰密集度产品; Markus等12在NASA Team的基础上, 加上89 GHz垂直和水平极化亮温数据, 提出NASA Team 2算法; Hao13通过引入AMSR-E 6.9 GHz数据, 对已有的NASA Team算法进行改进, 提高了多年冰结果的计算精度; Lomax等14利用SSM/I 85.5 GHz数据, 提出Lomax算法来计算海冰密集度, Lomax算法是基于Spencer等15在1989年提出的极化纠正亮温PCT (polarization corrected temperature)概念的基础上得到的。而近年来又有学者基于其他不同层面或其他数据对海冰密集度算法做出新的研究, 如Kern等16通过对熔池的研究进而分析其对夏季微波亮温数据以及海冰密集度结果的影响, 并改善海冰密集度的反演精度; Korosov等17提出一种基于欧拉平流方案对不同冰龄的海冰分布进行估算的算法; Gabarro等18提出一种用最大似然估计法结合亮度温差的SMOS估算北极海冰密集度的新方法; Ye等19通过引入再分析空气温度数据提高海冰密集度反演结果精度; 王欢欢20根据89 GHz频段亮温数据针对一年冰、 多年冰以及海水的不同特性, 提出了一种反演多年密集度的方法; Karvonen21提出一种使用神经网络框架内的双极化SAR数据来估算海冰密集度; Kasapoğlu22的一项研究中, 针对每个入射角范围选择SAR图像特征, 并使用线性混合模型通过被动微波成像来反演海冰密集度; ASI算法(ARTIST sea ice algorithm)产生于1998年“北极辐射和湍流交换的研究”项目中, 该算法是Kaleschke等23在Svendsen等24提出的接近90 GHz频段海冰反演算法的基础上改进而来的, 最初用于对北极海冰边缘的大气边界层做中尺度数值模拟, 它的一个优势在于相比于其他使用85 GHz频段数据的算法, 并不需要另外的数据输入, 可以直接根据测得的亮温数据进行结果反演25, 且与利用其他通道的海冰密集度算法有相似的结果26。ASI算法最初应用于SSM/I数据, 随着AMSR-E (Advanced Microwave Scanning Radiometer for the Earth Observing System)数据的产生, Spreen等27将ASI算法应用到AMSR-E数据, 并得到相应的海冰密集度计算公式; 苏洁等28对基于AMSR-E 89 GHz频段微波数据的ASI算法中北极纯冰和纯水代表区域的P1P0年平均值进行插值推算, 并对海冰密集度反演公式进行了修正。Lu等29通过辐射传输模型, 评估了风、 水蒸气以及液态水对地表亮温的影响, 进而对原有的ASI算法进行改进, 改进后的算法大大提高了冰水混合区域的海冰密集度的精度。张翔等30提出一种利用HY-2卫星扫描微波辐射计数据反演极地海冰密集度的算法; 石立坚等31利用HY-2卫星微波辐射计数据, 基于NASA Team算法反演了北极区域海冰密集度。罗丽程等32利用高分辨率SAR数据, 并采用k均值聚类算法提取海冰信息, 其与传统的海冰密集度产品结合, 可构造信息更丰富的多尺度冰情图。闻斌等33基于SMAP卫星雷达数据, 分析L波段雷达后向散射系数、 极化比和归一化极化差与海冰密集度之间的相关性, 并建立人工神经网络算法对海冰密集度进行反演研究, 结果表明SMAP海冰密集度产品与现有业务化海冰密集度产品具有很好的一致性。宋翔宇等34将亮温极化梯度率和光谱梯度率引入基于全约束最小二乘法(fully constrained least squares, FCLS)的海冰密集度估算方法, 结果表明改进的FCLS算法具有较高的精确度, 且对多年冰的估算优于一年冰。周颖等35利用FY-3 MERSI数据提出一种提取北极海冰分布信息的算法。

但目前应用FY-3 MWRI数据反演得到的海冰密集度产品较少, 更没有将ASI海冰密集度算法应用到MWRI数据上的研究。与SSM/I数据相比, MWRI数据拥有更宽的频率域; 与AMSR-E数据相比, MWRI在频率分布与空间分辨率上与前者相差不大。随着我国极地战略的推进以及冰上丝绸之路的提出, 将目前比较成熟的海冰密集度算法应用到国产卫星数据中具有极其重要的意义。本研究为利用国产FY-3 MWRI数据反演北极海冰密集度做出部分工作, 可为我国极地战略的推进以及北极航道的研究提供技术保障和数据支持。

1 数据来源

本文使用2016年1月FY-3卫星MWRI数据对北极地区的海冰密集度进行分析研究。FY-3卫星能够获取全球范围内全天时间段多光谱通道的大气、 地表和海表特性参数。MWRI是装载于FY-3卫星上的传感器, 该仪器为全功率双极化微波成像仪器36。MWRI主要探测地表的热辐射, 对地表参数比较敏感37, 在多个特定的频率上接收来自地球表面与大气的水平和垂直两个极化的电磁辐射, 提供全球性信息资料, 主要用于反演降雨、 云和大气中水汽含量、 地表植被、 海冰、 雪覆盖及海面油污等分布情况, 为天气预报和环境监测等提供资料38。FY-3 MWRI数据相对于SSM/I数据拥有更多的通道数, 同样满足ASI算法的应用需求, 用于极地海冰密集度反演(http://satellite.nsmc.org.cn/PortalSite/Data/Satellite.aspx#)。表1比较了MWRI、 SSM/I (搭载卫星Defense Meteorological Satellite Program, DMSP)以及AMSR-E (搭载卫星Aqua)三种传感器参数。

表1   微波传感器MWRI、 SSM/I及AMSR-E部分参数比较

Table 1  Some parameters of microwave sensor MWRI, SSM/I and AMSR-E

传感器搭载卫星频率/GHz地面分辨率/(km×km)极化方式轨道高度/km幅宽/km
MWRIFY-3C10.651×85V/H8361 400
MWRIFY-3C18.730×50V/H8361 400
MWRIFY-3C23.827×45V/H8361 400
MWRIFY-3C36.518×30V/H8361 400
MWRIFY-3C899×15V/H8361 400
SSM/IDMSP19.369×43V/H8331 394
SSM/IDMSP22.360×40V8331 394
SSM/IDMSP3737×29V/H8331 394
SSM/IDMSP85.515×13V/H8331 394
AMSR-EAqua6.943×75V/H7051 445
AMSR-EAqua10.6529×51V/H7051 445
AMSR-EAqua18.716×27V/H7051 445
AMSR-EAqua23.818×32V/H7051 445
AMSR-EAqua36.58×14V/H7051 445
AMSR-EAqua899×10V/H7051 445

新窗口打开| 下载CSV


2 算法

2.1 ASI算法

ASI算法是在1998年“北极辐射和湍流交换的研究”项目中产生的, 最开始提出这个算法的目的是为了更好地利用更高分辨率的SSM/I 85 GHz数据, 对北极海冰边缘的大气边界层做中尺度数值模拟24。ASI算法根据极化差异(polarization difference)来计算海冰密集度。

P=Tbv-Tbh

式中: Tbv为垂向极化的亮温; Tbh为水平极化的亮温。考虑到大气对极化差异的影响acP可以写为

P=Pse- τ(1.1e- τ-0.11)=Psac

式中: τ为大气的不透明度; Ps为被观测物体表面的极化差异。

Ps=CPs,i+ (1-CPs,w

式中: Ps,iPs,w分别为海冰表面和海水面的极化差异; C为海冰密集度。对于开阔水面, C=0, 极化差异P0=a0Ps,w; 对于100%海冰覆盖区域, C=1, 极化差异P1=a1Ps,i。对式P=[CPs,i+(1-CPs,wacC=0和C=1处做泰勒展开, 并忽略高阶项可得

P=a0CPs,i-Ps,w)+P0C→0)
P=a1C-1)(Ps,i-Ps,w)+P1C→1)

于是得到

C=(PP0-1)(Ps,wPs,i-Ps,w) C→0)
C=PP1+(PP1-1)(Ps,wPs,i-Ps,w)C→1)

对于北极地区而言, Ps,wPs,i-Ps,w的典型取值为-1.1423

为了更详细地反演从0到100%所有的海冰密集度, 需要对式(6)~(7)的解进行合理插值。因此, 选择一个三阶的多项式来拟合从0到100%的海冰密集度。

C=d3p3+d2p2+d1p+d0

利用式(6)~(7)以及两式对P的一阶导数, 计算得到d0d1d2d3

P03P02P01P13P12P113P02P1103P122P110d3d2d1d0=01- 1.14- 0.14

所以, 如果已知无冰海面和100%海冰覆盖区域的特征点的P0P1, 就可以得到d0d1d2d3, 最终得到海冰密集度C

2.2 基于FY-3 MWRI数据的ASI算法

ASI算法针对不同的数据核心参数会有所变化, 将ASI算法应用到MWRI数据中需要重新确定海冰密集度计算公式中的相关参数。通过选取典型海域的样本点, 获取纯水与纯冰的系点值, 进而确定海冰密集度的计算公式, 其中确定纯水系点值的典型区域须满足该区域全年均为开阔水域, 确定纯冰系点值的典型区域须满足该区域全年均为海冰。因此选取图1中两个矩形区域, 其中海冰区域位于加拿大北极群岛以北的区域, 海水区域位于格陵兰海冰外缘线以南区域。

图1

图1   样本点区域位置选取(来源: https://seaiceuni-bremen.de/start/)

Fig.1   Location of the sample points (source: https://seaice.uni-bremen.de/start/)


所选样本点会受到天气、 数据本身以及选取区域的影响, 因此在整个选取样本点的过程中, 对2016年全年的MWRI日数据89 GHz频段极化差进行样本点统计。通过分别计算两个区域内所有样本点的极化差异、 进行概率分布统计, 选取逐日最大概率发生的值作为当天的极化差异值, 再对全年的极化差异值进行求均值, 得到如图2~3所示的统计图。最终确定基于FY-3 MWRI数据89 GHz的纯水(P0)系点值平均结果为47.6 K, 纯冰(P1)系点值平均结果为10.8 K。对式(9)进行插值计算, 得到基于FY-3 MWRI数据的ASI海冰密集度C的表达式为

C=1.29×10- 5p3- 1.28×10- 3p2+1.01×10- 2p+1.02

图2

图2   海水样本点统计结果

Fig.2   The statistical results of seawater sample points


图3

图3   海冰样本点统计结果

Fig.3   The statistical results of sea ice sample points


使用高频数据时, 大气中云层液态水和海面水蒸气对亮温的影响较为为严重。特别是无冰海面上的气旋会减弱海水表面的极化差异, 使这部分海水像元接近海冰的极化差异, 可能导致该部分海水被误认为海冰。因此, 使用天气滤波器来去除在水面因误差计算出的海冰是非常有必要的。截至目前, 所有的天气滤波器基本上都是利用了低分辨率的低频段数据。Comiso10在1986年提出通过使用36.5 GHz和18.7 GHz两个频段的梯度率GR(gradient ratio)来减弱天气的影响, 因为在水面上GR(37/19)大于0, 而在冰面上GR(37/19)接近0或者小于0。1995年, Comiso11对这个天气滤波器进行了改进, 除了应用原来的GR(37/19)外, 又增加了GR(23/19), 因为23与19频段的梯度率对大气中的水蒸气更为敏感, 而GR(37/19)可以弱化海冰的边缘由于表面温度的变化而造成的影响。步骤如下:

①利用36.5 GHz以及18.7 GHz的梯度率GR过滤掉云中液态水的影响。

GR(37/19)=[Tb (37V)-Tb (19V)]/[Tb (37V)+Tb (19V)]
GR(37/19)≥0.045 (C→0)

②使用23.8 GHz以及18.7 GHz频段的梯度率GR去除海水表面大量水蒸气的影响。

GR(23/19)=[Tb (23V)-Tb (19V)]/[Tb (23V)+Tb (19V)]
GR(23/19)≥0.04 (C→0)

经过以上两次过滤, 天气对海冰造成误判的影响大部分已经消除掉, 此时只有某些特别极端的天气因素才会造成小部分影响。

天气滤波器的阈值也会随数据的变化而改变, 使用Otsu算法来确定FY-3 MWRI数据的天气滤波器阈值。Otsu算法, 又被称为大津法, 该方法是根据图像数据值之间的方差选取阈值。其原理如下:

首先假设阈值具有L级数据的图像划分为两类, 即C00,TC1T+1,  L-1。图像各数据级对应的概率为

pi=niN

式中: N为原图的总像素个数; ni 为数据级为i的像素个数; 且有pi ≥0, i=0L-1pi=1C0C1类的概率分别为

w0=pr(C0)=i=0Tw(T)
w1=pr(C1)=i=T+1L-1pi=1-w(T)

C0C1类的均值分别为

u0=i=0Tipiw0=u(T)w(T)
u1=i=TL-1ipiw1=u¯-u(T)1-w(T)

其中, u¯=i=0L-1ipi, 为图像的均值。C0C1的类间方差为

σg2=w0w1u1-u02

海冰和海水的最佳分类阈值T*应使类间方差最大, 即

T*=argmax0TL-1σg2

通过运算, 得出GR(37/19)天气滤波器的阈值为0.05, 即当GR(37/19)≥0.05时, 令其海冰密集度为0; GR(23/19)天气滤波器的阈值为0.045, 即当GR(23/19)≥0.045时, 令其海冰密集度为0。通过天气滤波器处理, 绝大部分受天气影响而导致计算错误的海冰都被过滤掉了, 被纠正为海水。

3 结果与验证

以2016年1月FY-3 MWRI数据为例对北极海冰密集度进行结果反演, 并经过天气滤波器的处理, 去除掉受天气影响的无效海冰像元, 得到每日平均海冰密集度结果。通过对1月31天每日海冰密集度结果进行累加并求取平均结果, 最终得到如图4(c)所示的月平均海冰密集度结果。为了充分验证本文结果的精度以及基于FY-3 MWRI数据提取的海冰与海水系点值的准确性, 分别从美国国家冰雪数据中心(National Snow and Ice Data Center, NSIDC,https://nsidc.org/data)以及德国不莱梅大学(University Bremen, Germany,https://seaice.uni-bremen.de/start/)获取同期的海冰密集度产品。其中, 图4(b)为NSIDC提供的由SSM/I数据基于NASA TEAM算法反演出的结果; 图4(a)为不莱梅提供的由AMSR-E数据基于ASI算法反演出的结果。

图4

图4   三种算法下2016年1月北极平均海冰密集度

Fig. 4   Averaged monthly sea ice concentration in Arctic in January, 2016: ASI algorithm based on University Bremen (a), NASA Team algorithm based on NSIDC (b), ASI algorithm based on FY-3 MWRI (c)


MWRI数据得到的结果与不莱梅提供的结果均使用了ASI算法, 二者的海冰密集度分布状况基本一致, 北冰洋中心区域均为高密集度海冰; NSIDC提供的结果基于NASA Team算法使用SSM/I数据反演得到, 而NASA Team算法使用的是低频数据, 明显低于ASI算法, 海冰边缘区域的低密集度海冰分布较多。为了更加明确地对比分析以上结果, 以及更为直观地对比三种海冰面积结果在1月的变化趋势以及差异, 制作了海冰面积变化趋势图(图5)和统计表(表2)。

图5

图5   三种算法下2016年1月北极海冰面积日变化

Fig. 5   Daily variations of sea ice area in January, 2016, quoted from FY-3 MWRI, University Bremen, and NSIDC


表2   三种算法下日平均海冰面积和海冰密集度

Table 2  Daily sea ice area and sea ice concentration resulted from three algorithms

日期

(年-月-日)

海冰面积/106 km2海冰密集度FY-3与不莱梅海冰 面积同比差值/%FY-3与NSIDC海冰 面积同比差值/%
FY-3不莱梅NSIDCFY-3不莱梅NSIDC
2016-01-0111.21511.49311.1910.9260.9270.855-2.4800.210
2016-01-0211.28511.47011.2430.9230.9270.851-1.6410.369
2016-01-0311.31811.45411.2050.9220.9270.853-1.1951.000
2016-01-0411.45411.48711.2430.9160.9290.855-0.2811.843
2016-01-0511.49911.39411.2620.9120.9240.8520.9122.060
2016-01-0611.47211.45511.3420.9140.9210.8470.1481.130
2016-01-0711.59211.59811.4210.9190.9230.848-0.0541.470
2016-01-0811.71211.70411.5580.9250.9270.8560.0681.313
2016-01-0911.72311.86711.7290.9120.9310.862-1.226-0.057
2016-01-1011.81012.04911.8680.9130.9300.860-2.025-0.495
2016-01-1111.90512.05111.8900.9120.9290.856-1.2260.128
2016-01-1211.91612.10811.8190.9150.9310.857-1.6080.818
2016-01-1312.00412.12811.8780.9240.9280.850-1.0331.055
2016-01-1412.04812.20611.8940.9130.9340.857-1.3101.280
2016-01-1512.05412.17911.9290.9190.9320.858-1.0331.034
2016-01-1612.01112.21912.0160.9230.9320.859-1.729-0.045
2016-01-1712.01212.31112.0650.9190.9340.863-2.488-0.443
2016-01-1812.02012.21912.1000.9140.9310.861-1.652-0.666
2016-01-1911.92112.15712.0890.9140.9270.862-1.987-1.410
2016-01-2012.01012.20312.1210.9100.9250.862-1.606-0.928
2016-01-2112.05412.23312.1770.9080.9270.862-1.484-1.021
2016-01-2212.06112.24412.1970.9130.9260.864-1.520-1.131
2016-01-2312.08412.27912.2000.9200.9290.866-1.617-0.960
2016-01-2412.09812.37612.1690.9180.9250.859-2.297-0.581
2016-01-2512.09112.46612.1480.9180.9280.857-3.100-0.470
2016-01-2612.12512.50112.1780.9130.9270.850-3.097-0.434
2016-01-2712.14612.51512.2100.9150.9320.854-3.036-0.526
2016-01-2812.18112.46712.2270.9100.9270.854-2.349-0.383
2016-01-2912.23512.45112.2640.9170.9240.850-1.768-0.239
2016-01-3012.33112.52512.3090.9150.9260.854-1.5750.177
2016-01-3112.30312.50712.2980.9110.9260.855-1.6590.039
平均值11.89412.07511.7490.9160.9280.856-1.5211.219
极大值12.33112.52512.3090.9260.9340.866-1.5730.178

新窗口打开| 下载CSV


表2可知:

(1) 平均海冰面积

MWRI数据得到的1月平均海冰面积为11.894×106 km2, 与NSIDC的11.749×106 km2同比差值为1.219%, 与不莱梅的12.508×106 km2同比差值为-1.521%。

(2) 单日最大海冰面积

MWRI数据得到的1月单日海冰面积最大为12.331×106 km2, 与NSIDC的12.309×106 km2同比差值为0.178%, 与不莱梅的12.525×106 km2同比差值为-1.573%。海冰极大值日均为1月30日。

(3) 单日最小海冰面积

MWRI数据得到的1月单日海冰面积最小为11.215×106 km2, 与NSIDC的11.191×106 km2同比差值为0.213%, 与不莱梅的11.394×106 km2同比差值为-1.596%。MWRI结果与NSIDC结果海冰极小值日均为1月1日, 不莱梅结果极小值日为1月5日。

(4) 平均海冰密集度

MWRI数据得到的1月平均海冰密集度为0.916, 与NSIDC的0.856同比差值为6.553%, 与不莱梅的0.928同比差值为-1.310%。MWRI结果更接近不莱梅的结果, 平均海冰面积三者基本一致。

图5看出, 本文基于FY-3 MWRI数据并使用ASI算法反演得到的结果, 与NSIDC以及不莱梅提供的结果都比较接近。从1月31天的海冰面积变化趋势来看, 三种结果均呈现出平稳增加的趋势。MWRI结果显示, 海冰面积平均以0.0313×106 km2·d-1的速度增加, NSIDC与不莱梅的结果分别以0.0389×106 km2·d-1、 0.0388×106 km2·d-1的速度增加。

为进一步验证本文结果反演精度, 引用空间分辨率为250 m的风云卫星中分辨率光谱成像仪(MERSI)数据进行结果对比。基于本文方法得到的2018年6月10日北极地区的海冰密集度结果如图6所示, 选取大利亚霍夫岛附近海域(图中矩形区域)作为研究样本, 用MERSI数据进行对比验证。

图6

图6   基于FY-3 MWRI数据得到的2018年6月10日北极海冰密集度

Fig.6   Sea ice concentration in Arctic on June 10, 2018, based on FY-3 MWRI data


从风云气象卫星中心获取2018年6月10日MERSI数据, 经预处理后并进行研究区域裁剪, 得到如图7(a)所示的可见光原始图像。本文采用目前较为普遍使用的单通道系点法反演海冰密集度信息, 算法公式39

C=0AAw(A-Aw)/(Ai-Aw)Aw<A<Ai100%AAi

式中: C为海冰密集度; A为像元的实际反照率; AiAw分别为纯冰与纯水的反照率。

图7

图7   北极大利亚霍夫岛附近海域光学验证结果

Fig.7   Optical verification results of sea area near Daliahov Island, Arctic: original image of MERSI data band 1 (a), sea ice concentration based on MERSI data (b), and sea ice concentration based on MWRI data (c)


通过对MERSI波段1的数据进行处理, 得到如图7(b)所示的海冰密集度结果。并对基于MWRI数据的2018年6月10日的海冰密集度结果进行裁剪, 得到如图7(c)所示的研究区域的海冰密集度结果。

对比图7(b)与图7(c)可知, 两者的海冰外缘线基本一致, 其中由光学数据得到的结果, 在陆地以及海岛边缘区域反演得更为精确, 而MWRI数据得到的结果对于浮冰、 流冰的识别较弱。经统计, 基于MWRI得到的平均海冰密集度为0.678, 海冰面积为160 929 km2; 基于MERSI得到的平均海冰密集度为0.714, 海冰面积为177 465 km2。二者海冰密集度相差5.029%, 海冰面积相差9.318%。综上可知, 本文基于MWRI数据运用ASI算法得到的海冰密集度结果, 其海冰边缘线与光学数据得到的结果基本一致, 海冰密集度以及海冰面积略低于光学数据得到的结果。

4 结论

以2016年1月共31天的FY-3 MWRI数据为例, 通过选取典型区域的样本点进行筛选统计得到基于MWRI数据的海冰与海水的系点值, 进而确定基于ASI算法的海冰密集度计算公式。再引入Otsu算法确定基于MWRI数据的天气滤波器阈值, 最终得到2016年1月北极海冰密集度以及海冰面积结果。利用美国国家冰雪数据中心(NSIDC)与德国不莱梅大学业务化海冰密集度产品进行对比验证, 发现不莱梅月平均海冰密集度以及海冰面积均最高, NSIDC结果均最低, MWRI结果介于两者之间。但三种结果的月平均海冰面积差别均不大, 而MWRI结果与不莱梅较为接近。此外, 以2018年6月10日数据为例, 截取大利亚霍夫岛附近海域为研究区域, 引入MERSI光学数据反演了海冰密集度。结果表明, 两者海冰外缘线基本一致, MWRI数据得到的海冰密集度以及海冰面积略低于光学数据结果。因此, 基于FY-3 MWRI数据所得到的ASI海冰密集度算法计算公式以及结果是可信的, 将国产FY-3 MWRI数据应用于ASI算法进行海冰密集度反演是可行的。

参考文献

Comiso J C.

Large decadal decline of the Arctic multiyear ice cover

[J]. Journal of Climate, 2011254): 1176 - 1193.

[本文引用: 1]

Stroeve J CSerreze M CHolland M Met al.

The Arctic’s rapidly shrinking sea ice cover: a research synthesis

[J]. Climatic Change, 20121103/4): 1005 - 1027.

[本文引用: 1]

Budikovan D.

Role of Arctic sea ice in global atmospheric circulation: a review

[J]. Global and Planetary Change, 2009683): 149 - 163.

[本文引用: 1]

Holland M MBitz C M.

Polar amplification of climate change in coupled models

[J]. Climate Dynamics, 2003213/4): 221 - 232.

[本文引用: 1]

Cohen JScreen J AFurtado J Cet al.

Recent Arctic amplification and extreme mid-latitude weather

[J]. Nature Geoscience, 201779): 627 - 637.

Vihma T.

Effects of Arctic sea ice decline on weather and climate: a review

[J]. Surveys in Geophysics, 2014355): 1175 - 1214.

[本文引用: 1]

Thomas D NDieckmann G S. Sea ice[M]. 2nd ed. Hoboken, New Jersey, USAWiley-Blackwell Publishing2010.

[本文引用: 1]

Ulaby F.

Microwave remote sensing: active and passive: microwave remote sensing fundamentals and radiometry

[J]. Remote Sensing A, 198125): 1223 - 1227.

[本文引用: 1]

Cavalieri D JGloersen PCampbell W J.

Determination of sea ice parameters with the NIMBUS 7 SMMR

[J]. Journal of Geophysical Research: Atmospheres, 198489D4): 5355 - 5369.

[本文引用: 1]

Comiso J C.

Characteristics of arctic winter sea ice from satellite multispectral microwave observations

[J]. Journal of Geophysical Research: Oceans, 198691C1): 975 - 994.

[本文引用: 2]

Comiso J C.

SSM/I sea ice concentrations using the Bootstrap algorithm: NASA RP-1380

[R]. Washington, D.C.National Aeronautics and Space Administration1995.

[本文引用: 2]

Markus TCavalieri D.

An enhancement of the NASA Team sea ice algorithm

[J]. IEEE Transactions on Geoscience and Remote Sensing, 2000381387 - 1398.

[本文引用: 1]

Hao Guanghua.

A study of multiyear ice concentration retrieval algorithms using AMSR-E data

[J]. Acta Oceanologica Sinica, 2015349): 102 - 109.

[本文引用: 1]

Lomax A SLubin DWhritner R H.

The potential for interpreting total and multiyear ice concentrations in SSM/I 85.5 GHz imagery

[J]. Remote Sensing of Environment, 19955413 - 26.

[本文引用: 1]

Spencer R WGoodman H MHood R E.

Precipitation retrieval over land and ocean with SSM/I: identification and characteristics of the scattering signal

[J]. Journal of Atmospheric and Oceanic Technology, 19896254 - 273.

[本文引用: 1]

Kern SRösel APedersen L Tet al.

The impact of melt ponds on summertime microwave brightness temperatures and sea-ice concentrations

[J]. Cryosphere Discussions, 2016105): 2217 - 2239.

[本文引用: 1]

Korosov A ARampal PPedersen L Tet al.

A new tracking algorithm for sea ice age distribution estimation

[J]. Cryosphere Discussions, 2018126): 2073 - 2085.

[本文引用: 1]

Gabarro CTuriel AElosegui Pet al.

New methodology to estimate Arctic sea ice concentration from SMOS combining brightness temperature differences in a maximum-likelihood estimator

[J]. Cryosphere Discussions, 2017114): 1987 - 2002.

[本文引用: 1]

Ye YHeygster GShokr M.

Improving multiyear ice concentration estimates with reanalysis air temperatures

[J]. IEEE Transactions on Geoscience and Remote Sensing, 2016545): 2602 - 2614.

[本文引用: 1]

Wang Huanhuan.

Multiyear ice retrieval using passive microwave remote sensing radiometer AMSR-E 89 GHz data

[D]. Qingdao, ShandongOcean University of China2009.

[本文引用: 1]

王欢欢.

利用遥感卫星AMSR-E 89 GHz频段的数据反演北极多年冰密集度

[D]. 山东青岛中国海洋大学2009.

[本文引用: 1]

Karvonen J.

Baltic sea ice concentration estimation based on C-band HH-polarized SAR data

[J]. IEEE Transactions on Geoscience and Remote Sensing, 201256): 1874 - 1884.

[本文引用: 1]

Kasapoğlu N G.

Sea ice concentration retrieval using composite scan SAR features in a SAR data assimilation process

[J]. IEEE Geoscience and Remote Sensing Letters, 20141112): 2085 - 2089.

[本文引用: 1]

Kaleschke LLiipkes CVihma Tet al.

SSM/I sea ice remote sensing for mesoscale ocean-atmosphere interaction analysis

[J]. Canadian Journal of Remote Sensing, 2001275): 526 - 537.

[本文引用: 2]

Svendsen EMauler CGrenfell T C.

A model for retrieving total sea ice concentration from a spaceborne dual-polarization passive microwave instrument operating near 90 GHz

[J]. International Journal on Remote Sensing, 1987810): 1479 - 1487.

[本文引用: 2]

Kern S.

A new method for medium-resolution sea ice analysis using weather-influence corrected special sensor microwave/imager 85 GHz data

[J]. International Journal of Remote Sensing, 20042521): 4555 - 4582.

[本文引用: 1]

Kem SKaleschke LClausi D A.

A comparison of two 85 GHz SSM/I ice concentration algorithm with AVHRR and ERS-2SAR imagery

[J]. IEEE Transactions on Geoscience and Remote Sensing, 20034110): 2294 - 2306.

[本文引用: 1]

Spreen GKaleschke LHeygster G.

Sea ice remote sensing using AMSR-E 89-GHz channels

[J/OL]. Journal of Geophysical Research: Oceans, 2008113(C2) [2019-03-04]. .

[本文引用: 1]

Su JieHao GuanghuaYe Xinxinet al.

The experiment and validation of sea ice concentration AMSR-E retrieval algorithm in polar region

[J]. Journal of Remote Sensing, 2013173): 495 - 513.

[本文引用: 1]

苏洁郝光华叶鑫欣.

极区海冰密集度AMSR-E数据反演算法的试验与验证

[J]. 遥感学报, 2013173): 495 - 513.

[本文引用: 1]

Lu JHeygster GSpreen G.

Atmospheric correction of sea ice concentration retrieval for 89 GHz AMSR-E observations

[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2018115): 1442 - 1457.

[本文引用: 1]

Zhang XiangWang ZhenzhanChen Hua.

A sea ice concentration algorithm based on HY-2 scanning radiometer data

[J]. Remote Sensing Technology and Application, 2012276): 912 - 918.

[本文引用: 1]

张翔王振占谌华.

一种利用HY-2卫星扫描微波辐射计数据反演极地海冰密集度的算法

[J]. 遥感技术与应用, 2012276): 912 - 918.

[本文引用: 1]

Shi LijianWang QimaoZou Binet al.

Arctic sea ice concentration retrieval using HY-2 radiometer data

[J]. Chinese Journal of Polar Research, 20144): 410 - 417.

[本文引用: 1]

石立坚王其茂邹斌.

利用海洋(HY-2)卫星微波辐射计数据反演北极区域海冰密集度

[J]. 极地研究, 20144): 410 - 417.

[本文引用: 1]

Luo LichengZhang WenqiHu Yonget al.

Information extraction of arctic sea ice based on Sentinel-1 satellite

[J]. Aerospace Shanghai, 20183): 16 - 22.

[本文引用: 1]

罗丽程张文奇胡勇.

基于Sentinel-1卫星的北冰洋海冰信息提取

[J]. 上海航天, 20183): 16 - 22.

[本文引用: 1]

Wen BinZhou XuanZhong Jinsonget al.

Sea ice concentration retrieval based on SMAP radar data

[J]. Acta Oceanologica Sinica, 2018406): 29 - 39.

[本文引用: 1]

闻斌周旋种劲松.

基于SMAP卫星雷达资料的海冰密集度反演技术研究

[J]. 海洋学报, 2018406): 29 - 39.

[本文引用: 1]

Song XiangyuLiu TingtingWang Zeminet al.

Estimation of Antarctic sea ice concentration based on an improved FCLS algorithm

[J]. Chinese Journal of Polar Research, 20181): 67 - 76.

[本文引用: 1]

宋翔宇刘婷婷王泽民.

基于改进FCLS算法的南极海冰密集度估算及算法比较

[J]. 极地研究, 20181): 67 - 76.

[本文引用: 1]

Zhou YingGong CailanHu Yonget al.

Extraction of sea ice information in the Arctic Ocean by using FY-3/MERSI data

[J]. Journal of Atmospheric and Environmental Optics, 201381): 54 - 59.

[本文引用: 1]

周颖巩彩兰胡勇.

风云三号MERSI数据提取北冰洋海冰信息方法研究

[J]. 大气与环境光学学报, 201381): 54 - 59.

[本文引用: 1]

Yang HuLi XiaoqingYou Ranet al.

Environmental data records from FengYun-3B microwave radiation imager

[J]. Advances in Meteorological Science and Technology, 201334): 136 - 143.

[本文引用: 1]

杨虎李小青游然.

风云三号微波成像仪定标精度评价及业务产品介绍

[J]. 气象科技进展, 201334): 136 - 143.

[本文引用: 1]

Zou Xiaolei.

Introduction to microwave imager radiance observations from polar-orbiting meteorological satellites

[J]. Advances in Meteorological Science and Technology, 201223): 45 - 50.

[本文引用: 1]

邹晓蕾.

极轨气象卫星微波成像仪资料

[J]. 气象科技进展, 201223): 45 - 50.

[本文引用: 1]

Guan MinYang Zhongdong.

Research on remote sensing image location method based on FY-3 microwave imager

[J]. Journal of Remote Sensing, 2009133): 469 - 474.

[本文引用: 1]

关敏杨忠东.

FY-3微波成像仪遥感图像地理定位方法研究

[J]. 遥感学报, 2009133): 469 - 474.

[本文引用: 1]

Steffen KMaslanik J A.

Comparison of Nimbus 7 scanning multichannel microwave radiometer radiance and derived sea ice concentrations with Landsat imagery for the north water area of Baffin Bay

[J]. Journal of Geophysical Research: Oceans, 198893C9): 10769 - 10781.

[本文引用: 1]

/