史飛飛，肖建設(shè)，李林，曹曉云，李曉東，李紅梅

1.青海省氣象科學(xué)研究所，西寧 810008

2.青海省防災(zāi)減災(zāi)重點實驗室，西寧 810008

3.青海省氣象局，西寧 810008

4.青海師范大學(xué)地理科學(xué)學(xué)院，西寧 810008

5.青海省氣候中心，西寧 810008

引 言

青藏高原氣候惡劣、災(zāi)害頻發(fā)，是氣候變化的敏感區(qū)和生態(tài)脆弱帶，青藏高原牧區(qū)作為我國重要的特色畜牧業(yè)生產(chǎn)地之一，覆蓋了中國四大牧區(qū)中的青海和西藏高寒牧區(qū)，平均海拔在4000米以上，共包括108個牧區(qū)半牧區(qū)縣[1-3]。牧區(qū)雪災(zāi)是由于積雪過厚，且長時間掩埋草場致使牲畜無法正常采食，導(dǎo)致大量牲畜掉膘和死亡的自然災(zāi)害[4]。長久以來青藏高原大部分牧區(qū)內(nèi)牛羊仍以草原粗獷放牧形式為主，一般當(dāng)雪層深度較厚，積雪持續(xù)較長時會對家畜采食產(chǎn)生影響，并在無飼草儲備或飼草儲備不足的牧區(qū)極有可能出現(xiàn)雪災(zāi)[5-6]。當(dāng)前在全球氣候變化背景下雪災(zāi)發(fā)生頻率和程度均呈增加趨勢，并且隨著牧區(qū)畜牧規(guī)模的擴大，其成災(zāi)損失和社會影響也明顯加大，已嚴(yán)重制約了青藏高原畜牧業(yè)的可持續(xù)發(fā)展[7-8]。例如在 2019年三江源地區(qū)發(fā)生雪災(zāi)期間，受其影響的家畜數(shù)量達到了231.20萬羊單位，占到全區(qū)受災(zāi)家畜標(biāo)準(zhǔn)羊單位總量的78.30%[9]。

目前，在牧區(qū)雪災(zāi)監(jiān)測與預(yù)警方面，張學(xué)通等[10]和王瑋等[11]先后在青藏高原地區(qū)系統(tǒng)開展了牧區(qū)雪災(zāi)的監(jiān)測與預(yù)警研究工作，所構(gòu)建的雪災(zāi)風(fēng)險預(yù)警模型精度能達到85%以上；周秉榮等[12-13]研制了青海高原牧區(qū)雪災(zāi)綜合預(yù)警評估模型并在雪災(zāi)服務(wù)中進行了應(yīng)用。而在雪災(zāi)風(fēng)險評估和區(qū)劃方面，李紅梅等[14]和韓炳宏等[15]進行了青海高原牧區(qū)雪災(zāi)的風(fēng)險區(qū)劃，其結(jié)果與實際雪災(zāi)發(fā)生狀況較為相符。上述研究中的雪災(zāi)案例均主要通過查閱《中國氣象災(zāi)害年鑒》[16]和《中國氣象災(zāi)害大典》[17]等資料進行收集，因歷史災(zāi)情數(shù)據(jù)的發(fā)布往往具有滯后性，同時不同資料間對災(zāi)情信息統(tǒng)計的詳實程度也存在差異，將會對雪災(zāi)風(fēng)險評估理論模型構(gòu)建和雪災(zāi)風(fēng)險區(qū)劃等研究工作產(chǎn)生不利影響[18]。針對該問題，氣象部門發(fā)展出一種基于地面氣象站點資料進行青藏高原牧區(qū)雪災(zāi)實時監(jiān)測的方法，并已在牧區(qū)雪災(zāi)等級監(jiān)測和預(yù)警等方面得到較好應(yīng)用，同時其監(jiān)測結(jié)果也可作為雪災(zāi)數(shù)據(jù)的有益補充。

數(shù)據(jù)集構(gòu)建過程中綜合考慮了雪災(zāi)致災(zāi)因子危險性和雪災(zāi)承災(zāi)體易損性，利用了青藏高原 109個地面氣象臺站逐日積雪觀測資料和地面災(zāi)情資料，參照《氣象災(zāi)害分級指標(biāo)》[19]（DB63/T372-2018）中的雪災(zāi)監(jiān)測標(biāo)準(zhǔn)以及牲畜死亡狀況制定了雪災(zāi)等級分級指標(biāo)，并據(jù)此統(tǒng)計了近20年來青藏高原牧區(qū)發(fā)生輕度、中度、重度和特重度等級雪災(zāi)的日數(shù)，形成的雪災(zāi)數(shù)據(jù)集可為青藏高原雪災(zāi)風(fēng)險評估、雪災(zāi)風(fēng)險預(yù)警和風(fēng)險區(qū)劃提供數(shù)據(jù)支撐。

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

1.1 數(shù)據(jù)生產(chǎn)流程

從中國氣象局綜合氣象信息共享平臺（China Integrated Meteorological Information Sharing System,CIMISS）獲取青藏高原109個氣象臺站2000-2020年逐日雪深觀測資料，并已在CIMISS平臺中按照氣象行業(yè)規(guī)范進行了質(zhì)量控制[20]，地面氣象站點分布狀況見圖1。地面災(zāi)情資料主要包括：牲畜存欄數(shù)、牲畜死亡率、受災(zāi)人口數(shù)、財產(chǎn)損失等，上述資料來自于《中國氣象災(zāi)害大典》[17]、《青海統(tǒng)計年鑒》[21]、《西藏統(tǒng)計年鑒》[22]以及各地氣象部門收集和調(diào)查資料等。根據(jù)文中構(gòu)建的雪災(zāi)等級劃分標(biāo)準(zhǔn)，在Matlab2018b軟件中通過編寫程序，逐年統(tǒng)計了近20年來青藏高原各地發(fā)生輕度、中度、重度和特重度等級雪災(zāi)的日數(shù)，最后依據(jù)政府部門發(fā)布的災(zāi)情數(shù)據(jù)和現(xiàn)有文獻資料對數(shù)據(jù)集質(zhì)量進行評估，具體流程如圖2所示。

圖1 青藏高原氣象站點分布Figure 1 Distribution of meteorological stations on the Tibetan Plateau

圖2 數(shù)據(jù)處理流程Figure 2 Data processing flow

1.2 雪災(zāi)等級分級標(biāo)準(zhǔn)

在《氣象災(zāi)害分級指標(biāo)》[19]（DB63/T372-2018）中主要采用了積雪深度和積雪維持時間兩項指標(biāo)用于評估雪災(zāi)的異常程度，而李紅梅等[14]研究指出高原內(nèi)牧區(qū)雪災(zāi)對牲畜影響程度最大，因而可使用牲畜死亡率來表征雪災(zāi)的損失程度，并且在積雪深度、積雪維持時間和牲畜死亡率間具有較好相關(guān)性，綜合上述指標(biāo)能更好地開展青藏高原牧區(qū)雪災(zāi)評估。本文在上述研究成果的基礎(chǔ)上制定了青藏高原牧區(qū)雪災(zāi)等級分級指標(biāo)（具體見表1），其中冬季指當(dāng)年10月至次年2月，而春季指次年3月至5月，積雪深度是指從積雪表面到地面的垂直深度，積雪維持時間是指地面積雪穩(wěn)定維持的連續(xù)日數(shù)，牲畜死亡率是指受雪災(zāi)造成的大牲畜（主要包括牛、羊和馬）死亡數(shù)占年初的大牲畜存欄數(shù)的比例，上述指標(biāo)均可通過各地氣象部門的觀測資料和氣象災(zāi)害數(shù)據(jù)庫統(tǒng)計獲得。

表1 雪災(zāi)等級指標(biāo)Table 1 Indices of snow disaster levels

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

2000-2020年青藏高原牧區(qū)雪災(zāi)等級數(shù)據(jù)集包含4個數(shù)據(jù)文件，分別為輕度雪災(zāi)發(fā)生日數(shù)、中度雪災(zāi)發(fā)生日數(shù)、重度雪災(zāi)發(fā)生日數(shù)和特重度雪災(zāi)發(fā)生日數(shù)，格式為shp，數(shù)據(jù)量42 KB。文件屬性表中“sta”為站點編號，“總和”為2000-2020年對應(yīng)級別雪災(zāi)發(fā)生的總?cè)諗?shù)，而以年份（例如2000）命名的屬性列中記錄了該年內(nèi)對應(yīng)級別雪災(zāi)發(fā)生的日數(shù)，圖3為數(shù)據(jù)集的制圖示例。

圖3 2000-2020年青藏高原牧區(qū)雪災(zāi)發(fā)生日數(shù)Figure 3 Days of snow disasters in pastoral areas of the Tibetan Plateau during 2000-2020

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

青藏高原內(nèi)有109個有人值守的氣象臺站，均分布在有從事畜牧業(yè)生產(chǎn)活動的牧業(yè)和半牧業(yè)市縣內(nèi)。在雪災(zāi)發(fā)生后地方政府會及時通報災(zāi)情狀況，而當(dāng)?shù)貧庀蟛块T則會通過氣象災(zāi)害管理系統(tǒng)(http://10.1.64.146/disaster/)對政府發(fā)布的災(zāi)情信息以及實地調(diào)查的災(zāi)損情況進行上報，主要包括災(zāi)害類別、災(zāi)害時間、經(jīng)濟損失和人員傷亡等情況。為驗證本數(shù)據(jù)集的精度和可靠性，篩選出117個有詳細地面災(zāi)情調(diào)查資料的雪災(zāi)案例，逐個對比案例中雪災(zāi)發(fā)生地區(qū)內(nèi)政府部門發(fā)布的雪災(zāi)發(fā)生日數(shù)和本數(shù)據(jù)集中統(tǒng)計的雪災(zāi)發(fā)生日數(shù)，總體上兩者相關(guān)程度較高，R2達到了0.84，平均絕對誤差保持在5 d。同時從2000-2020年青藏高原不同等級雪災(zāi)發(fā)生日數(shù)空間分布圖來看（圖3），輕度雪災(zāi)發(fā)生日數(shù)較多的是清水河、甘德、嘉黎、達日、石渠、錯那和德令哈等地區(qū)，累計輕度雪災(zāi)日數(shù)超過100 d；中度雪災(zāi)發(fā)生日數(shù)較多的是瑪多、甘德、清水河、聶拉木、普蘭和錯那等地區(qū)，累計中度雪災(zāi)日數(shù)超過75 d；重度雪災(zāi)發(fā)生日數(shù)較多的是清水河、聶拉木、普蘭、瑪多、雜多和德令哈地區(qū)，累計重度雪災(zāi)日數(shù)超過75 d；特重度雪災(zāi)發(fā)生日數(shù)最多的是聶拉木、普蘭、錯那、清水河和瑪多地區(qū)，累計特重度雪災(zāi)日數(shù)超過75 d。總體上三江源地區(qū)易發(fā)生中度及以上牧區(qū)雪災(zāi)，青藏高原的喀喇昆侖山脈、昆侖山脈、喜馬拉雅山脈和阿尼瑪卿山等地易發(fā)生特重度雪災(zāi)，而柴達木盆地、青海湖流域、青海省東部農(nóng)業(yè)區(qū)和藏南谷地雪災(zāi)風(fēng)險較低，本數(shù)據(jù)集展現(xiàn)的雪災(zāi)分布特征與李紅梅等[14]和韓炳宏等[15]的研究結(jié)果一致。

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

本數(shù)據(jù)以2000-2020年青藏高原地面氣象臺站積雪觀測資料和災(zāi)情資料為基礎(chǔ)，以雪災(zāi)等級分級指標(biāo)作為判別標(biāo)準(zhǔn)統(tǒng)計了青藏高原牧區(qū)發(fā)生輕度、中度、重度和特重度等級雪災(zāi)的日數(shù)，經(jīng)驗證數(shù)據(jù)精度較高，并能較好反映青藏高原不同等級雪災(zāi)的分布空間特征，將對農(nóng)牧業(yè)生產(chǎn)和牧區(qū)雪災(zāi)防治等方面工作發(fā)揮重要作用。