青藏高原风云3C卫星积雪产品评估

青藏高原FY-3C卫星积雪产品评估

闵文彬,1, 彭骏1, 李施颖2

1.中国气象局成都高原气象研究所/高原与盆地暴雨旱涝灾害四川重点实验室,成都 610071

2.四川省气象探测数据中心,成都 610071

选取2018年10月1日—2019年4月30日作为积雪研究期,利用青藏高原地区118个气象台站的积雪观测数据,对风云三号C星(FY-3C)多仪器融合的雪盖(multi-sensor synergy snow cover, MULSS_SNC)和微波成像仪的雪水当量(microwave radiation imager snow water equivalent, MWRIX_SWE)产品进行评估,以了解产品的区域可靠性。结果表明: MULSS_SNC和MWRIX_SWE雪盖判识准确率分别为87.18%和72.32%,召回率分别为66.67%和49.63%,误判率分别为12.81%和27.68%,漏判率分别为33.33%和50.37%; 对于积雪混合像元或积雪深度不足0.5 cm的卫星像元,MULSS_SNC和MWRIX_SWE都倾向于无雪判识,1 cm以下雪盖漏判率高达60%以上; MULSS_SNC在积雪深度达2 cm以上时,召回率可达到89.09%,而MWRIX_SWE产品在积雪深度达5 cm以上时,积雪判识召回率才达到63.37%。MWRIX_SWE反演的青藏高原地区雪深与地面观测值误差很大,不存在线性正相关关系,不建议直接使用。

0 引言

积雪作为冰冻圈的重要组成之一,在气候变化和水资源中扮演着重要角色。积雪的累积和消融影响水资源循环,积雪的高反射率、低导热率特性影响地表能量平衡和气候系统变化。因此,准确、及时掌握积雪的时空变化对于气候变化预测、水资源管理等至关重要。目前,积雪观测数据主要来自传统的地面站点观测和卫星遥感观测。在我国青藏高原地区,积雪观测站点少,空间分布不均匀,观测时间不连续,不足以提供区域积雪分布和变化特征信息。卫星遥感成为该地区积雪监测的理想手段,能够实现对积雪宏观、快速和周期性的监测。人们可利用光学遥感数据提取雪盖信息[1,2,3,4,5,6],利用微波遥感数据反演雪深和雪水当量[7,8,9,10],及时掌握研究区的积雪动态,认识研究区的积雪和雪水当量的时空分布特征[11,12,13]。然而,可见光遥感依赖于太阳光,受云层影响严重,在天空有云状况下无法工作; 被动微波遥感虽具有穿透性及全天候工作特点,不受天气状况的影响,但空间分辨率低,主要用于全球雪深、积雪覆盖范围和雪水当量的研究,在区域性的积雪动态监测和积雪研究中,卫星遥感积雪产品精度还有许多需要改进的地方,不同传感器、不同反演算法之间存在较大偏差[14,15,16,17,18,19,20]。所以,目前遥感积雪产品在区域性应用之前,必须经过真实性检验和精度评估,它是评价遥感产品质量、可靠性和适用性的唯一手段,是提高遥感产品精度、改善遥感产品质量的主要依据,更是推动遥感产品应用范围和应用水平的重要保障。

我国风云三号系列气象卫星不仅提供遥感观测资料,还提供大量的业务化产品,包括积雪监测产品,并实现了广泛的数据共享。本研究以地面积雪观测数据为参考,开展FY-3C多仪器融合数据的积雪覆盖产品(multi-sensor synergy snow cover, MULSS_SNC)和微波成像仪的雪深雪水当量产品(microwave radiation imager snow water equivalent, MWRIX_SWE)在青藏高原地区的精度检验和评估分析,这对于风云三号气象卫星积雪产品在研究区的应用是必不可少的,具有非常重要的意义。

1 数据与评估方法

1.1 数据来源

本文使用的卫星反演积雪产品取自国家气象科学数据共享服务平台(http: //data.cma.cn/site/index.html),要检验的FY-3C多仪器融合的积雪覆盖产品(MULSS _SNC)空间分辨率为1 km,FY-3C微波成像仪的雪深和雪水当量产品(MWRIX_SWE)空间分辨率为25 km。作为检验参考值的地面积雪观测数据取自全国综合气象信息共享平台(China integrated meteorological information sharing system, CIMISS),包括西藏、青海和四川境内的118个气象台站,台站分布见图1,海拔跨度较大,从1 370 m到4 800 m。

图1   青藏高原积雪观测站点分布示意图

1.2 评估方法

1.2.1 地面观测参考值

根据地面雪深观测规范,当气象站四周视野地面被雪覆盖超过一半时观测雪深,当雪深达到或超过5 cm时观测雪压。台站雪深观测记录为99990表示有少量积雪,地面积雪深度观测记录为0表示气象站四周视野地面被雪覆盖超过一半,但雪深不足0.5 cm。本文把台站雪深观测记录为99990和0,即有少量积雪或雪深不足0.5 cm的气象站对应的卫星像元看作部分积雪像元,把雪深观测记录大于或等于0.5 cm的气象站对应的卫星像元看作积雪像元,把无积雪气象站对应的卫星像元看作陆地像元。

1.2.2 卫星积雪覆盖产品评估

1.2.3 卫星反演雪深和雪水当量产品评估

2 FY-3C积雪产品评估分析

2.1 MULSS_SNC雪盖产品

表1为MULSS_SNC判别结果,2018年10月—2019年4月118个台站共有25 016个地面观测样本。

表1   MULSS_SNC判识结果

假设MULSS_SNC雪盖产品的云像元判识准确,不考虑卫星误将积雪判识为云的情况,也不评估部分积雪的混合像元情况,只对2 270个积雪和陆地卫星像元进行评估。根据公式(1)-(5)可计算出,MULSS_SNC雪盖产品总体精度为75.37%、雪盖判识准确率为87.18%、召回率为66.67%、雪盖误判率为12.81%、漏判率为33.33%。部分积雪像元只有2.63%被判识为积雪像元,97.37%被判识为无雪的陆地像元。

2.2 MWRIX_SWE积雪判识

将MWRIX_SWE产品的反演值大于0的卫星像元作为积雪像元看待,利用地面积雪观测数据对其有/无积雪判识准确率进行分析。从表2可见,有9 121个卫星像元缺数据,参与评估的只有15 895个像元。MWRIX_SWE产品判识有4 085个积雪像元,其中998个积雪像元与地面观测一致,误判382个,判识的2 705个积雪像元对应地面为部分积雪; MWRIX_SWE产品判识的陆地像元有11 810个,其中漏判了1 013个积雪像元,正确判识陆地像元613个,将10 184个部分积雪像元判识为了陆地。去掉缺数据的卫星像元,只统计参与评估的数据,根据公式(1)-(5)计算出,MWRIX_SWE产品总体精度为53.59%; 积雪像元判识准确率为72.32%、召回率为49.63%、误判率为27.68%、漏判率为 50.37%。79.01%的部分积雪像元被判识为陆地。

表2   MWRIX_SWE判识结果

2.3 不同深度的积雪判识

表3给出了FY3C积雪产品不同积雪深度的检验结果,表中不同深度积雪判识的召回率、漏判率和漏判百分比都是针对地面观测样本而言的,漏判百分比为某一深度的卫星漏判的积雪像元数占卫星漏判的积雪像元总数的比例。从分析结果可见,积雪深度在[0.5,1]cm的有460个像元,MULSS_SNC产品正确判识出168个,漏掉292个,积雪判识召回率只有36.52%,漏判率高达63.48%,漏判数占总漏判像元数的67.59%; 积雪深度在(1,2] cm之间的185个像元中,MULSS_SNC产品正确判识出116个,漏掉69个,判识召回率为62.70%,漏判率为37.30%,漏判百分比为15.97%; 当积雪深度达到2 cm以上时,MULSS_SNC产品对不同深度的雪盖判识召回率都达到80%以上,漏判率小于20%。积雪深度2 cm以上的651个像元中,MULSS产品正确判识出580个,漏掉71个,判识召回率达到89.09%,漏判率为10.91%,漏判百分比为16.44%。

表3   不同积雪深度的FY-3C积雪产品检测结果

由表3对参与评估的2 011个MWRIX_SWE像元的分析结果可见,积雪深度在[0.5,1] cm的有723个像元,MWRIX_SWE产品正确判识出238个,漏掉485个,判识召回率只有32.92%,漏判率高达67.08%,漏判数占总漏判像元数的47.88%; 积雪深度在(1,2] cm之间的286个像元中,MWRIX_SWE产品正确判识出125个,漏掉161个,判识召回率为43.71%,漏判率为56.29%,漏判百分比为15.89%; 积雪深度在(2,3] cm之间的214个像元中,MWRIX_SWE产品正确判识出101个,漏掉113个,判识召回率为47.20%,漏判率为52.80%,漏判百分比为11.15%。MWRIX_SWE对积雪深度在3 cm以下的积雪像元漏报率很高,超过了50%。积雪深度5 cm以上1 002个像元中,MWRIX_SWE产品正确判识出635个,漏掉367个,判识召回率达到63.37%,漏判率36.63%,漏判百分比36.23%。

通过表4和图2的MULSS_SNC和MWRIX_SWE产品检验结果对比可见,对于相同深度的积雪像元,MULSS_SNC产品雪盖判识召回率都比MWRIX_SWE产品的高,漏报率低; 当积雪深度达到2 cm以上时,MULSS_SNC产品雪盖判识召回率接近90%,比MWRIX_SWE高出25%,漏报百分比低50%。

表4   MULSS_SNC和MWRIX_SW合成的雪盖产品的判识结果

2.4 MULSS_SNC与MWRIX_SWE积雪合成产品雪盖判识准确率

假如MWRIX_SWE卫星像元地表单一,将其25 km像元反演值均匀插值到1 km,与MULSS_SNC产品空间匹配,利用MWRIX_SWE监测的积雪像元替代MULSS_SNC的云像元,减小MULSS_SNC产品的云影响,得到MULSS_SNC和MWRIX_SWE合成产品。从表4可见,合成积雪产品总体精度为64.85%、雪盖判识准确率为83.08%、召回率为61.54%、误判率为16.92%、漏判率为38.46%,将89.77%的部分积雪像元判为陆地像元。虽然雪盖合成产品的精度不比MULSS_SNC高,但能判识出的积雪像元明显多于单一用MULSS_SNC和MWRIX_SWE的判识结果,可在一定程度上弥补MULSS_SNC的云和MWRIX_SWE无数据的影响,获得更多的积雪信息。当然,MULSS_SNC和MWRIX_SWE产品不同分辨率的空间匹配势必造成判识结果差异,这里将微波与光学产品简单结合应用,只是下一步工作的思路测验而已。

2.5 MWRIX_SWE的雪深和雪水当量定量分析

由图3的地面观测与卫星反演雪深散点图可知,当地面雪深达到60 cm以上时,MWRIX_SWE反演雪深值却不足10 cm,严重偏小; 反之,MWRIX_SWE反演雪深值达到60 cm以上时,地面观测值也不足10 cm。即,不论是地面观测还是卫星反演雪深达到60 cm以上,二者偏差都达到50 cm以上,偏差很大。而且地面观测与卫星反演雪深的判定系数R2为0.022 8,说明二者没有正相关关系。将地面雪深数据作为参考值,误差分析得知,MWRIX_SWE反演雪深的平均偏差为6.05 cm,偏差值范围-151 cm~72.4 cm; 绝对误差为11.76 cm; 均方根误差18.41 cm; 平均相对误差高达486.16%。

接下来,在998组MWRIX_SWE积雪对比数据中,剔除出卫星反演或地面观测积雪深度大于60 cm的25组数据,对余下973组数据进行统计分析,偏差范围在-37~27.4 cm,平均偏差为7.54 cm; 绝对误差为9.78 cm; 相关系数为-0.04。图4为60 cm以下地面观测与卫星反演雪深散点图,从中可以看出,地面观测与卫星反演雪深值之间相关性仍然不好。误差和相关分析结果都表明MWRIX_SWE的雪深产品在青藏高原地区不可直接使用。分析结果表明,被动微波积雪产品不受天气状况的影响,但空间分辨率低,主要用于全球雪深、积雪覆盖范围和雪水当量的研究,在区域性的积雪动态监测中还存在较大偏差[12]。

鉴于以上分析结果,目前FY3C卫星反演雪深产品可靠较差,而且地面雪压观测值与MWRIX_SWE的雪水当量值有效配对样本数较少,只有22个,本文就不再定量评估雪水当量产品。

3 结论

利用青藏高原118个气象站积雪观测数据,对FY-3C的MULSS_SNC雪盖产品和MWRIX_SWE雪深与雪水当量产品进行评估分析,得出如下结论:

1)MULSS_SNC产品的总体精度为75.37%,积雪判识准确率为87.18%,召回率为66.67%,雪盖误判率为12.81%,雪盖漏判率为33.33%,97.37%的部分积雪混合像元被判识为无雪。

2)MWRIX_SWE产品对像元有无积雪判识的总体精度为53.59%,积雪像元判识准确率为72.32%,召回率为49.63%,雪盖误判率为27.68%,雪盖漏判率为50.37%,79.01%部分积雪混合像元被判识为无雪。

3)利用MWRIX_SWE产品对MULSS_SNC产品进行消云处理后的雪盖合成产品总体精度为64.85%、积雪准确率为83.08%、雪盖召回率为61.54%、误判率为16.92%、漏判率为38.46%,89.77%的部分积雪混合像元判识为无雪。

4)积雪深度在1 cm以下时,MULSS_SNC产品判识召回率只有36.52%,漏判率高达63.48%,漏判数占总漏判像元数的67.59%; 积雪深度达2 cm以上时,MULSS_SNC产品判识召回率达到89.09%,漏判率为10.91%,漏判百分比为16.44%。

5)积雪深度在1 cm以下时,MWRIX_SWE产品判识召回率只有32.92%,漏判率高达67.08%,漏判百分比为47.88%; 积雪深度在3 cm以下时,漏报率都很高,超过了50%。积雪深度达到5 cm以上时,判识召回率才达到63.37%,漏判率为36.63%,漏判百分比为36.23%。

6)MWRIX_SWE卫星反演雪深的平均偏差为6.05 cm,偏差值范围-151~72.4 cm; 绝对误差为11.76 cm; 均方根误差18.41 cm; 平均相对误差高达486.16%; 卫星反演雪深值与地面观测值之间相关系数为-0.15; 不论是地面观测值,还是MWRIX_SWE卫星反演雪深值达到60 cm以上,二者偏差都达到50 cm以上。