彭九慧， 王 多， 童 俊， 楊慶紅

(1.承德市氣象局，河北 承德 067000； 2.河北石油職業(yè)技術(shù)大學(xué)，河北 承德 067000)

引 言

霜凍是我國發(fā)生范圍較廣的一種農(nóng)業(yè)氣象災(zāi)害,其發(fā)生頻次多、影響范圍大，對農(nóng)業(yè)、林果業(yè)危害嚴(yán)重，常造成嚴(yán)重的經(jīng)濟損失。因此，進行霜凍災(zāi)害風(fēng)險區(qū)劃，制定綜合管理對策，對霜凍災(zāi)害高風(fēng)險地區(qū)進行風(fēng)險管理和制定相應(yīng)的應(yīng)對措施，是防災(zāi)減災(zāi)的有效途徑[1-3]。近年來，各地的科技工作者，針對霜凍災(zāi)害的風(fēng)險評估與風(fēng)險區(qū)劃作了大量工作，積累了寶貴的經(jīng)驗。如李軍玲等[4]選取冬小麥晚霜凍有害積寒和頻度、冬小麥種植面積、灌耕比和單位提灌面積機井?dāng)?shù)等風(fēng)險評價指標(biāo),建立了冬小麥晚霜凍風(fēng)險評價模型。陳凱奇等[5]研究了遼寧省初霜凍日期的變化特征及其對玉米種植的影響，并對霜凍災(zāi)害風(fēng)險進行了評估。 嚴(yán)應(yīng)存等[6]采用圖表法和指數(shù)乘積法，研究了青海省霜凍災(zāi)害變化特征及風(fēng)險區(qū)劃。還有很多學(xué)者[7-10]研究了霜凍出現(xiàn)時間、霜凍強度與作物受災(zāi)程度的關(guān)系。氣象行業(yè)標(biāo)準(zhǔn)(作物霜凍害等級，2008)規(guī)定, 霜凍害是由于日最低氣溫下降，致使植株莖、葉溫度下降到0 ℃或0 ℃以下，使正在生長發(fā)育的作物受到凍傷，從而導(dǎo)致減產(chǎn)、品質(zhì)下降或絕收的現(xiàn)象[11]。以往的研究多以地面0 cm的最低溫度或日最低氣溫作為霜凍災(zāi)害等級的評價標(biāo)準(zhǔn), 把霜凍特別是終霜凍出現(xiàn)的時間作為霜凍災(zāi)害風(fēng)險評價因子的較少。通過分析承德市1965-2018年54年的終霜凍災(zāi)情發(fā)現(xiàn),即使最低氣溫或地面0 cm的溫度較高(在0 ℃左右),如果終霜出現(xiàn)的時間夠晚, 也會造成嚴(yán)重的經(jīng)濟損失，有時出現(xiàn)時間晚的終霜甚至比溫度低造成的損失更大。

如2005年5月6-8日,平泉市遭受嚴(yán)重的終霜凍災(zāi)害, 受災(zāi)人口9.5萬人，作物凍死、凍傷面積6500公頃，直接經(jīng)濟損失2000多萬元，當(dāng)日平泉市氣溫0 ℃以上, 地面0 cm最低氣溫在-1 ℃以上，終霜日期比多年平均日期晚10天。

承德市地處山區(qū)，幾乎每年春天都有因終霜凍造成的經(jīng)濟損失。精細(xì)化的終霜災(zāi)害風(fēng)險區(qū)劃將是調(diào)整農(nóng)業(yè)種植計劃和品種選擇的重要依據(jù)。本研究通過構(gòu)建以時間系數(shù)和強度系數(shù)加權(quán)組合的綜合霜凍指數(shù)評價終霜凍災(zāi)害風(fēng)險等級的方法，增強霜凍災(zāi)害等級劃分的科學(xué)性。在此基礎(chǔ)上，以致災(zāi)因子的危險性、孕災(zāi)環(huán)境敏感性、承載體的脆弱性為評價因子，基于GIS空間分析技術(shù)，制作精細(xì)化的終霜凍災(zāi)害的風(fēng)險區(qū)劃圖[12],以期為農(nóng)業(yè)防災(zāi)減災(zāi)提供決策支撐。

1 資料來源

選取1965-2018年承德市近54年9個縣區(qū)春季逐日地面0 cm最低溫度觀測資料，承德市民政局提供的205個鄉(xiāng)鎮(zhèn)的54年各區(qū)縣89個終霜凍災(zāi)情損失資料，基準(zhǔn)年為2014年精細(xì)到鄉(xiāng)鎮(zhèn)的全市玉米播種面積和蘋果、梨種植面積資料，1∶50000比例尺的承德市數(shù)字高程資料。在進行ArcGIS空間分析時，為了方便對每個圖層賦予不同權(quán)重進行計算疊加，將每個數(shù)據(jù)序列按公式(1)進行歸一化處理，所有的數(shù)據(jù)值均在0～1，柵格單元大小均設(shè)置為1 km×1 km。

2 資料處理方法

氣象行業(yè)標(biāo)準(zhǔn)是利用氣溫的高低對作物霜凍的強度進行分級。分析承德市多年歷史災(zāi)情數(shù)據(jù)發(fā)現(xiàn)，終霜出現(xiàn)的時間晚比氣溫強度對作物造成的損失更大。因此在對終霜凍災(zāi)害等級劃分時，增加了時間判別系數(shù)，由強度、時間兩項指標(biāo)根據(jù)災(zāi)情賦予不同的權(quán)重，形成綜合霜凍指數(shù)來對終霜凍災(zāi)害危險等級進行劃分。

2.1 終霜凍災(zāi)害的時間系數(shù)處理

以每年春季地面0 cm最低溫度最后一次達(dá)到0 ℃或以下出現(xiàn)的日期作為終霜日進行分析。將近54年逐年的終霜出現(xiàn)日期進行排序，以終霜出現(xiàn)的最早日期為1，從小到大排序，根據(jù)當(dāng)年出現(xiàn)的日期對應(yīng)的序號位置，得到當(dāng)年終霜災(zāi)害對應(yīng)的時間系數(shù)值。對逐年的時間系數(shù)值按公式(1)進行歸一化處理，得到為0～1的每年的終霜凍災(zāi)害時間系數(shù)(出現(xiàn)時間越晚，數(shù)值越大)：

(1)

其中，xi′與xi分別表示數(shù)據(jù)的量化值和原始值，xmax表示指標(biāo)i在所有數(shù)據(jù)中的最大值。

2.2 終霜凍災(zāi)害的強度系數(shù)處理

將54年終霜日地面0 cm最低溫度進行排序，因為大多數(shù)數(shù)值為負(fù)值，所有的數(shù)值均加上最小值的絕對值，將所有數(shù)據(jù)變?yōu)檎担賹λ械臄?shù)據(jù)取倒數(shù)后歸一化處理，從而得到介于0～1的每年的終霜凍災(zāi)害強度系數(shù)。地面0 cm最低溫度越低，數(shù)值越大。為了防止負(fù)值轉(zhuǎn)正值計算中出現(xiàn)數(shù)字0，最小值直接賦值0.01。

2.3 綜合霜凍指數(shù)的制定與分級

加權(quán)評價法是一種多因子評價法，基本原理是根據(jù)各因子相對于評價對象的重要程度確定其權(quán)重，加權(quán)綜合后用一個總的量化指數(shù)來表達(dá)評價對象優(yōu)劣的評價方法[13-14]。在本研究中，終霜凍災(zāi)害危險等級采用終霜凍時間系數(shù)和強度系數(shù)加權(quán)組合的綜合霜凍指數(shù)來表征。

分析54年災(zāi)情數(shù)據(jù)發(fā)現(xiàn)，承德市終霜凍出現(xiàn)的時間對霜災(zāi)損失影響較大，終霜日出現(xiàn)得越晚，造成的損失越大。因此參照災(zāi)情的實際情況，采用專家打分法,將時間系數(shù)的權(quán)重設(shè)定為0.7，強度系數(shù)的權(quán)重設(shè)定為0.3。

I=0.7X1+0.3X2

(2)

I為綜合霜凍指數(shù)，X1為時間系數(shù)，X2為強度系數(shù)。

計算近 54年各縣區(qū)綜合霜凍指數(shù)的期望值和標(biāo)準(zhǔn)差，采用期望值加標(biāo)準(zhǔn)差的正負(fù)0.5倍將霜凍災(zāi)害分為3級，利用標(biāo)準(zhǔn)差分級法將霜凍危險級別分為重霜凍、中霜凍和輕霜凍(表1)。

表1 承德市終霜凍災(zāi)害分級標(biāo)準(zhǔn)(I為綜合霜凍指數(shù))

計算綜合霜凍指數(shù)與89個災(zāi)情個例農(nóng)作物受災(zāi)面積的簡單相關(guān)系數(shù)為0.67, 通過了ɑ=0.01 的顯著性水平檢驗(圖1)。因此，采用綜合霜凍指數(shù)評價法，能根據(jù)災(zāi)情科學(xué)地劃分霜凍災(zāi)害危險等級。

圖1 綜合霜凍指數(shù)與農(nóng)作物受災(zāi)面積相關(guān)分析圖

3 終霜凍災(zāi)害風(fēng)險區(qū)劃

3.1 危險性評價

致災(zāi)因子的危險性往往是由致災(zāi)因子的強度和頻次兩方面的因素決定的[15-16]，而綜合霜凍指數(shù)反映了終霜出現(xiàn)的時間、強度的綜合特征，某地的綜合霜凍指數(shù)大，說明霜凍災(zāi)害致災(zāi)因子較活躍，災(zāi)害發(fā)生的概率大。根據(jù)綜合霜凍指數(shù)，分別計算各縣區(qū)的終霜重霜凍、中霜凍、輕霜凍的出現(xiàn)頻次，根據(jù)災(zāi)損率確定輕、中、重3級霜凍害的危險性權(quán)重比為0.3∶0.7∶1.0[17],歸一化處理后, 確定輕、中、重3級霜凍害的危險性權(quán)重分別為0.15、0.35和0.50。根據(jù)公式(3)，計算各地的終霜災(zāi)害危險性指數(shù)H，得到承德市終霜災(zāi)害的危險性分布(圖2)。

(3)

式中,H為終霜災(zāi)害危險性指數(shù)，Pi為終霜重霜凍、中霜凍、輕霜凍的出現(xiàn)頻次，Ii為輕、中、重3級霜凍害的危險性權(quán)重。

采用自然斷點分級(圖2)后可以看出，承德市終霜災(zāi)害的高危險區(qū)與次髙危險區(qū)主要分布在圍場縣、灤平縣 、隆化縣和豐寧縣中南部, 終霜凍的中危險區(qū)主要分布在豐寧縣北部、市區(qū)、平泉市、承德縣和興隆縣北部，終霜凍的低危險區(qū)分布在承德市南部的寬城縣和興隆縣大部。

圖2 承德市終霜災(zāi)害危險性指數(shù)分布圖

3.2 孕災(zāi)環(huán)境評價

孕災(zāi)環(huán)境主要指致災(zāi)因子、承災(zāi)體所處的自然環(huán)境和社會經(jīng)濟環(huán)境。分析承德市霜災(zāi)發(fā)生頻次與海拔高度關(guān)系發(fā)現(xiàn)，隨著海拔高度的升高，霜災(zāi)發(fā)生的次數(shù)增多。因此，確定孕災(zāi)環(huán)境敏感性主要考慮海拔高度與霜災(zāi)頻次2個評價因子(權(quán)重見表2)。將數(shù)字高程數(shù)據(jù)投影后柵格化、歸一化處理后，與歸一化的霜災(zāi)頻次加權(quán)求和，計算孕災(zāi)環(huán)境敏感性指數(shù)。對敏感性指數(shù)分級，得到終霜凍孕災(zāi)環(huán)境敏感性分布(圖3)。

圖3 承德市終霜災(zāi)害孕災(zāi)環(huán)境敏感性分布圖

從圖3可以發(fā)現(xiàn)，承德市終霜凍災(zāi)害高敏感區(qū)主要分布在海拔高度比較高的圍場縣和受災(zāi)頻次較高的灤平縣及隆化縣的東北部地區(qū)。

3.3 承災(zāi)體的易損性

承災(zāi)體的物理暴露性(亦稱為易損性)是指暴露在致險要素影響范圍內(nèi)的承險體數(shù)目或價值量。承德市54年的終霜凍災(zāi)情資料分析顯示，承德市霜凍受災(zāi)的對象主要是農(nóng)作物與經(jīng)濟林果，因此確定承災(zāi)體易損性主要考慮耕地面積、經(jīng)濟林面積。全市耕地種植面積排在第一位的是玉米，占總種植面積的70%左右。承德市大部區(qū)縣玉米在4月中下旬播種，9月中下旬收獲，終霜災(zāi)害有2/3出現(xiàn)在5月上旬，正值玉米的苗期，因此受終霜害影響最大的也是玉米。經(jīng)濟林中正值花期的蘋果和梨受終霜凍影響最大[18],因此選擇玉米種植面積、蘋果種植面積、梨種植面積作為終霜災(zāi)害易損性評價指標(biāo)，利用層次分析法[19-22]，得到各指標(biāo)的權(quán)重(表2)。將數(shù)據(jù)歸一化處理后,在ArcGIS中進行空間插值，通過柵格計算，加權(quán)求和[23]計算，得到承德市終霜災(zāi)害承災(zāi)體綜合易損性分布圖(圖4)(分級標(biāo)準(zhǔn)見表3)。

圖4 承德市終霜災(zāi)害承災(zāi)體易損性分布圖

從圖4中可以看出,承德市終霜凍災(zāi)害高易損區(qū)與次高易損區(qū)主要分布在玉米種植面積比較大的平泉市、豐寧縣中南部、灤平縣北部與蘋果、梨種植面積較大的承德縣、寬城縣北部。

3.4 制作承德市終霜災(zāi)害風(fēng)險區(qū)劃圖

氣象災(zāi)害風(fēng)險區(qū)劃由致災(zāi)因子危險性、孕災(zāi)環(huán)境敏感性、承災(zāi)體易損性、防災(zāi)減災(zāi)能力4個方面決定。霜凍災(zāi)害主要對大田和露天作物造成影響，受經(jīng)濟、技術(shù)、成本諸多因素影響，即使預(yù)測到霜凍災(zāi)害，也無法采取經(jīng)濟有效的措施防御[24]，因此在作終霜災(zāi)害風(fēng)險區(qū)劃時，選擇致災(zāi)因子危險性、孕災(zāi)環(huán)境敏感性、承災(zāi)體易損性3個指標(biāo)作為評價因子[25-26]。終霜災(zāi)害風(fēng)險評價函數(shù)表示為：終霜災(zāi)害風(fēng)險指數(shù)(R)=f(危險性、易損性、敏感性)。本文采用層次分析法通過構(gòu)造3階判斷矩陣,確定終霜災(zāi)害風(fēng)險3個評價指標(biāo)的權(quán)重并通過一致性檢驗(表2)。

表2 承德市終霜災(zāi)害風(fēng)險指數(shù)各評價因子與評價指標(biāo)權(quán)重

承德市每個格點的終霜災(zāi)害風(fēng)險指數(shù)=0.6522×歸一化危險性指數(shù)+0.107×歸一化易損性指數(shù)+0.265×歸一化敏感性指數(shù)。柵格計算后，利用自然斷點方法進行風(fēng)險等級劃分(分級標(biāo)準(zhǔn)見表3)，得到承德市終霜凍災(zāi)害風(fēng)險區(qū)劃(圖5)。

表3 終霜災(zāi)害風(fēng)險評價因子歸一化指數(shù)分級標(biāo)準(zhǔn)

從圖5中可以看出,承德市終霜災(zāi)害的高風(fēng)險區(qū)與次高風(fēng)險區(qū)主要位于平泉市、圍場縣、豐寧縣大部、灤平縣大部和隆化縣北部；中風(fēng)險區(qū)位于隆化縣大部、承德縣東部、興隆縣北部、灤平縣中部和豐寧縣部分地區(qū)；低風(fēng)險區(qū)位于市區(qū)、營子區(qū)、寬城縣、承德縣中西部和興隆縣中南部地區(qū)。

圖5 承德市終霜災(zāi)害風(fēng)險區(qū)劃圖

4 結(jié)論與討論

4.1 結(jié) 論

(1)本文在傳統(tǒng)劃分霜凍等級方法的基礎(chǔ)上，提出了根據(jù)歷史災(zāi)情、基于終霜出現(xiàn)時間和地面0 cm最低溫度，確定終霜出現(xiàn)時間系數(shù)與強度系數(shù)，構(gòu)建綜合霜凍指數(shù)，劃分山區(qū)終霜災(zāi)害危險等級的評價標(biāo)準(zhǔn)，綜合霜凍指數(shù)與農(nóng)作物受災(zāi)面積的相關(guān)系數(shù)為0.67。

(2)從致災(zāi)因子的危險性、承災(zāi)體的脆弱性、孕災(zāi)環(huán)境的敏感性出發(fā)，建立了承德市終霜凍災(zāi)害的風(fēng)險評價模型。分析終霜災(zāi)害風(fēng)險區(qū)劃圖發(fā)現(xiàn)，終霜災(zāi)害高風(fēng)險區(qū)主要位于海拔較高或受災(zāi)頻次較大的地區(qū)，承德市中北部區(qū)縣的終霜凍災(zāi)害風(fēng)險大于承德市南部地區(qū)的災(zāi)害風(fēng)險。

4.2 討 論

(1)在實際工作中,可以根據(jù)每年終霜出現(xiàn)的時間相對于最早日期的天數(shù)和預(yù)報(或?qū)崨r)的最低地溫值，計算出綜合霜凍指數(shù)，確定終霜災(zāi)害等級，再根據(jù)終霜災(zāi)害風(fēng)險區(qū)劃圖，確定受災(zāi)重點區(qū)域，為霜災(zāi)評估和政府防災(zāi)決策提供依據(jù)。

(2)由于霜凍觀測的空間分辨率較低，對于地形復(fù)雜、地域?qū)拸V的山區(qū) ，隨著氣象監(jiān)測網(wǎng)的逐步建立及監(jiān)測站點增多，可以將模型進一步細(xì)化。另外，受終霜災(zāi)情資料的限制及每年種植結(jié)構(gòu)變化的影響，區(qū)劃模型考慮的因子、風(fēng)險指標(biāo)的選擇，還需要進一步檢驗并動態(tài)更新。

(3)對于終霜災(zāi)害高風(fēng)險區(qū)或終霜頻次較高的地區(qū)，在生長期不變的情況下，可以選擇比較耐寒的玉米品種，以減少農(nóng)業(yè)損失。