郝曉華，趙琴，紀(jì)文政，王建，李弘毅

1.中國(guó)科學(xué)院西北生態(tài)環(huán)境資源研究院，蘭州 730000

2.中國(guó)科學(xué)院大學(xué)，北京 100049

引 言

積雪是冰凍圈的重要組成部分之一，也是全球氣候系統(tǒng)的重要組成部分[1-2]。全球約有98%的積雪分布于北半球，每年冬季，北半球陸地最大積雪范圍約為47×106km2，占北半球陸地面積的近50%[3-5]。積雪這種每年增加和消融的周期性變化規(guī)律以及年際間的變化通過(guò)表面能量收支平衡對(duì)氣候系統(tǒng)、碳匯水平、生態(tài)環(huán)境、能量交換、人文經(jīng)濟(jì)等產(chǎn)生直接或間接的影響[6-8]。因此準(zhǔn)確識(shí)別地表積雪信息、探究積雪物候變化對(duì)全球生態(tài)系統(tǒng)、氣候變化具有重大意義[9-11]。積雪物候（SCP, Snow Cover Phenology）是表征季節(jié)性積雪隨季節(jié)周期變化趨勢(shì)和變化規(guī)律的信息，包括積雪日數(shù)（Snow cover days, SCD）、積雪初日（Start of snow cover days, SCS）、積雪終日（Melt of snow cover days，SCM）等指標(biāo)[12]。我國(guó)包含三大季節(jié)性積雪區(qū)：東北內(nèi)蒙古積雪區(qū)，青藏高原積雪區(qū)和新疆北部積雪區(qū)，是中低緯度地區(qū)積雪分布最為廣泛的國(guó)家，因而獲取中國(guó)積雪物候數(shù)據(jù)集，進(jìn)而探究中國(guó)積雪物候變化，具有重要的氣候，水文、生態(tài)和人文意義[13-14]。

積雪物候研究主要是需要日降雪的資料，通常通過(guò)地面臺(tái)站觀測(cè)和衛(wèi)星遙感兩種技術(shù)手段獲取。地面臺(tái)站是早年獲取氣象數(shù)據(jù)、積雪空間分布的主要途徑，很多學(xué)者使用地面站點(diǎn)數(shù)據(jù)進(jìn)行了積雪物候相關(guān)研究[15-18]。但是氣象臺(tái)站多位于城市，代表性差，一些高海拔地區(qū)積雪分布廣泛但觀測(cè)點(diǎn)稀少，實(shí)測(cè)數(shù)據(jù)難以代表整個(gè)區(qū)域的積雪空間分布特征。遙感被廣泛應(yīng)用在積雪各項(xiàng)參數(shù)的反演研究中，用于獲取積雪物候的主要遙感產(chǎn)品可以分為光學(xué)遙感產(chǎn)品和微波遙感產(chǎn)品。光學(xué)遙感產(chǎn)品時(shí)空分辨率較高，但受云的影響較大，微波雪水當(dāng)量產(chǎn)品可以排除云的影響，但是空間分辨率較低，難以滿足積雪監(jiān)測(cè)[2]。目前更多的積雪產(chǎn)品主要是針對(duì)MODIS，雖然其空間分辨率較高，但時(shí)間序列較短，從2000年以后才可獲取數(shù)據(jù)。HAO等[19]基于NOAA的AVHRR CDR反射率數(shù)據(jù)，發(fā)展了一套中國(guó)1980-2020年5 km逐日無(wú)云積雪面積數(shù)據(jù)集，驗(yàn)證總體精度達(dá)到85%以上，與知名的JASMES AVHRR產(chǎn)品相比，整體精度提高約15%，漏分誤差從 60.8%下降到19.7%，錯(cuò)分誤差從31.9%下降到21.3%，能夠滿足積雪物候制備的需求。

因此，本文基于HAO等人的中國(guó)AVHRR逐日無(wú)云積雪面積產(chǎn)品[19]，獲取了中國(guó)水文年主要的積雪物候參數(shù)，積雪日數(shù)、積雪初日和積雪終日三個(gè)數(shù)據(jù)集，并利用地面臺(tái)站實(shí)測(cè)雪深資料對(duì)結(jié)果進(jìn)行精度評(píng)估。本數(shù)據(jù)集可服務(wù)于中國(guó)積雪時(shí)空變化分析，為氣候研究，水文管理，生態(tài)環(huán)境，人文經(jīng)濟(jì)等科學(xué)研究、工程建設(shè)以及社會(huì)服務(wù)提供基礎(chǔ)數(shù)據(jù)資料。

1 數(shù)據(jù)采集和處理方法

1.1 數(shù)據(jù)來(lái)源

中國(guó)1980-2020年5 km AVHRR逐日無(wú)云積雪面積產(chǎn)品[20]來(lái)自國(guó)家冰川凍土沙漠科學(xué)數(shù)據(jù)中心(www.ncdc.ac.cn)。該產(chǎn)品利用NOAA CDR AVHRR第四版本反射率數(shù)據(jù)作為輸入，通過(guò)云判別算法、雪判別算法和空缺值插補(bǔ)算法，最終獲取積雪面積產(chǎn)品（NIEER-CGF-AVHRR-SCE）。該數(shù)據(jù)集詳細(xì)信息如表1所示，其中1、2、3表示積雪，0表示非雪，4表示水體，255表示填充值。

表1 AVHRR逐日無(wú)云積雪面積產(chǎn)品詳細(xì)信息Table 1 NIEER-CGF-AVHRR-SCE dataset details

驗(yàn)證數(shù)據(jù)來(lái)自中國(guó)氣象局（http://data.cma.cn）地面臺(tái)站觀測(cè)的1980-2020年每日地面氣候積雪資料數(shù)據(jù)集。具體觀測(cè)方法為當(dāng)視野面積超過(guò)50%，利用米尺在北京時(shí)間8點(diǎn)人工測(cè)量氣象站雪深。觀測(cè)精確到1 cm，雪深小于1 cm的記錄為無(wú)雪，同時(shí)也記錄了缺測(cè)信息。最終，有積雪記錄的362個(gè)氣象站被用來(lái)對(duì)積雪物候數(shù)據(jù)集進(jìn)行精度評(píng)估。

1.2 處理方法

1.2.1 積雪物候參數(shù)定義

水文年定義為每年的9月1日到次年的8月31日，將原數(shù)據(jù)重分類成有雪無(wú)雪二值產(chǎn)品。

積雪日數(shù)定義為一個(gè)水文年上觀測(cè)到積雪的天數(shù)之和，公式如式（1）：

其中，i表示水文年中從9月1日起的第i天，n為一個(gè)水文年的總天數(shù)，范圍為0~365或366，t表示柵格像元值，t = 0表示無(wú)雪，t = 1表示有雪；

積雪初日定義為一個(gè)水文年中第一次出現(xiàn)連續(xù)5天是雪的首日對(duì)應(yīng)日期，即像元如果滿足公式（2）

則可得到

即該像元的積雪初日為從水文年9月1日起的第s天，其中，s表示水文年中從9月1日起的第s天，若存在不滿足公式（2）的像元，則該像元為0，因此SCS范圍為0~365/366。

積雪終日定義為一個(gè)水文年中最后出現(xiàn)連續(xù)5天是雪的終日對(duì)應(yīng)日期[18]，即像元如果滿足公式（4）

則可得到

即該像元的積雪終日為從水文年9月1日起的第n-m天，其中m為一個(gè)水文年中倒著數(shù)的第m天，n為全年天數(shù)，為365/366；若存在不滿足公式（2）的像元，則該像元為0，因此SCM范圍為0~365/366。

1.2.2 計(jì)算流程

積雪物候參數(shù)計(jì)算流程如下圖1所示，利用逐日無(wú)云AVHRR積雪面積產(chǎn)品作為輸入值，第一步對(duì)產(chǎn)品進(jìn)行預(yù)處理，將產(chǎn)品的有雪柵格值（t）賦為1，無(wú)雪柵格值賦為0，獲取二值積雪產(chǎn)品（有雪/非積雪）；第二步根據(jù) 1.2.1中的公式，計(jì)算每一個(gè)水文年的中國(guó)積雪日數(shù)、積雪初日和積雪終日。

圖1 積雪物候數(shù)據(jù)集制備流程示意圖Figure 1 Schematic of producing the dataset of snow phenology

2 數(shù)據(jù)樣本描述

本數(shù)據(jù)集包含積雪日數(shù)、積雪初日、積雪終日3個(gè)目錄，每個(gè)目錄下包含40個(gè)逐水文年積雪物候參數(shù)文件，數(shù)據(jù)量為47.1 MB。文件命名規(guī)則示例：NIEER_AVHRR_TTT_5000m_YYYY-YYYY.tif，其中NIEER表示中國(guó)科學(xué)院西北生態(tài)環(huán)境資源研究院，AVHRR表示先進(jìn)的超高分辨率輻射計(jì)，TTT表示不同的積雪物候參數(shù)（SCD為積雪日數(shù)，SCS為積雪初日，SCM為積雪終日），5000m表示空間分辨率，YYYY-YYYY表示水文年。積雪物候產(chǎn)品值代表的像元意義如表2所示，積雪物候參數(shù)取值范圍均為0~365或366，無(wú)效值為-1，水體為-255。其中積雪初日、終日0值為不討論區(qū)，其數(shù)值代表從每年9月1日開(kāi)始第n天（nth）代表的日期。圖2-4是以1980-1981水文年為例的產(chǎn)品示意圖，分別表示積雪日數(shù)、積雪初日和積雪終日。

表2 中國(guó)積雪物候產(chǎn)品詳細(xì)信息Table 2 Details of Chinese snow phenology dataset

圖2 中國(guó)積雪日數(shù)產(chǎn)品示意圖（審圖號(hào)：GS(2022)747號(hào)）Figure 2 Schematic of Chinese SCD dataset

圖3 中國(guó)積雪初日產(chǎn)品示意圖（審圖號(hào)：GS(2022)747號(hào)）Figure 3 Schematic of Chinese SCS dataset

圖4 中國(guó)積雪終日產(chǎn)品示意圖（審圖號(hào)：GS(2022)747號(hào)）Figure 4 Schematic of Chinese SCM dataset

3 數(shù)據(jù)質(zhì)量控制和評(píng)估

本數(shù)據(jù)集是基于積雪物候定義計(jì)算的遙感產(chǎn)品，為了保證結(jié)果的可靠性，使用站點(diǎn)實(shí)測(cè)雪深數(shù)據(jù)對(duì)積雪物候產(chǎn)品進(jìn)行精度驗(yàn)證，數(shù)據(jù)評(píng)估指標(biāo)包括相關(guān)系數(shù)（R2）、均方根誤差（RMSE）和平均絕對(duì)誤差（MAE）均方根誤差。站點(diǎn)篩選原則為：不考慮整個(gè)雪季期間有雪天數(shù) < 20天的站點(diǎn)。共篩選出代表性氣象站點(diǎn)100個(gè)，主要分布在三大積雪區(qū)，本研究雪深> 2 cm則認(rèn)為有雪，去除站點(diǎn)水文年存在nodata值的數(shù)據(jù)，站點(diǎn)雪深物候參數(shù)與本遙感數(shù)據(jù)集參數(shù)分布如圖5所示，圖5（a）所示，積雪日數(shù)的R2為0.86、RMSE為22.78天，MAE為13.26天；積雪初日和積雪終日由于定義的限制，剔除異常值，圖5（b）可以看到積雪初日的R2為0.76、RMSE為17.87天，MAE為7.51天；圖5（c）積雪終日的R2為0.94、RMSE為16.39天，MAE為7.76天，精度可靠。

圖5 積雪物候站點(diǎn)驗(yàn)證圖Figure 5 Site verification diagram of snow phenology

4 數(shù)據(jù)價(jià)值

數(shù)據(jù)集以中國(guó)1980-2019年積雪面積5公里逐日無(wú)云產(chǎn)品為基礎(chǔ)資料，計(jì)算了逐年的積雪物候參數(shù)，地面站點(diǎn)驗(yàn)證表明數(shù)據(jù)精度相對(duì)較高，可用于表征中國(guó)積雪時(shí)空分布及變化規(guī)律，可為中國(guó)境內(nèi)各區(qū)域氣候系統(tǒng)、碳匯水平、生態(tài)環(huán)境、能量交換、人文經(jīng)濟(jì)等提供基礎(chǔ)數(shù)據(jù)。

5 數(shù)據(jù)使用方法和建議

本數(shù)據(jù)集可從國(guó)家冰川凍土沙漠科學(xué)數(shù)據(jù)中心（http://www.ncdc.ac.cn）下載獲取。數(shù)據(jù)文件均為GeoTIFF格式，可以通過(guò)GIS與遙感相關(guān)的軟件如ENVI、GRASS、ArcGIS等直接進(jìn)行查看與應(yīng)用，或者使用編程語(yǔ)言等相應(yīng)的軟件進(jìn)行編譯讀取、計(jì)算分析等。對(duì)多年數(shù)據(jù)進(jìn)行空間疊加分析，可以得到積雪物候區(qū)域時(shí)空分布及年際變化趨勢(shì)，可結(jié)合氣象因素、人類活動(dòng)等進(jìn)行區(qū)域積雪變化的驅(qū)動(dòng)力分析，以期為生產(chǎn)及災(zāi)害預(yù)警等提供信息服務(wù)。

致 謝

感謝國(guó)家冰川凍土沙漠科學(xué)數(shù)據(jù)中心提供的1980-2020年積雪面積5公里逐日無(wú)云產(chǎn)品，感謝中國(guó)氣象局提供的地面氣候積雪資料日值數(shù)據(jù)集。