基于脆弱性曲線的河北省蘋果低溫凍害風(fēng)險評價

李婷,孫玉龍,陳笑娟,2*,徐藝芙

基于脆弱性曲線的河北省蘋果低溫凍害風(fēng)險評價

李婷1,孫玉龍1,陳笑娟1,2*,徐藝芙1

1. 河北省氣象災(zāi)害防御和環(huán)境氣象中心, 河北 石家莊 050021 2. 河北省氣象與生態(tài)環(huán)境重點實驗室, 河北 石家莊 050021

基于1981-2020年河北省142個國家站氣溫數(shù)據(jù)，利用優(yōu)勢分析法和信息熵賦權(quán)法分別構(gòu)建了蘋果開花—幼果期的低溫凍害致災(zāi)強度指數(shù)和致災(zāi)危險性指數(shù)，分析不同概率下的致災(zāi)強度?；陲L(fēng)險系統(tǒng)理論，通過“致災(zāi)強度—減產(chǎn)率”脆弱性曲線對蘋果低溫凍害受影響風(fēng)險進行評價。-2.0 ℃～-1.1 ℃之間的氣溫截點總體準確度在80%以上，確定致災(zāi)閾值為-1.4 ℃。致災(zāi)強度和致災(zāi)因子危險性自西北向東南逐漸降低，致災(zāi)因子危險性高等級分布在張家口西北部，其他大部分地區(qū)處于低和較低等級。蘋果低溫凍害高風(fēng)險區(qū)主要分布在單產(chǎn)水平較高的縣區(qū)。10年一遇致災(zāi)強度水平下，衡水阜城風(fēng)險最大，受影響產(chǎn)量為16.1 t/hm2。30年一遇時大部分種植區(qū)受影響產(chǎn)量增至2.1～4.0 t/hm2。50年一遇時，承德圍場風(fēng)險最高且呈自西北向東南呈逐級遞減，受影響產(chǎn)量最高為19.9 t/hm2。受致災(zāi)因子危險性、承災(zāi)體暴露度和脆弱性共同作用，承德圍場風(fēng)險較高，張家口東部和承德中西部地區(qū)致災(zāi)強度較高，但其暴露水平較低，風(fēng)險也較低。對于承德和張家口蘋果種植較多的縣區(qū)，應(yīng)考慮栽種抗逆性強的品種。

蘋果; 低溫凍害; 風(fēng)險評價

河北省是我國蘋果生產(chǎn)大省之一，蘋果作為河北省具有較強競爭力的優(yōu)勢果品之一，種植面積和產(chǎn)量均位于國內(nèi)前列。春季低溫凍害是影響產(chǎn)量、品質(zhì)和商品率的主要氣象災(zāi)害之一[1,2]。蘋果在不同生育期的耐寒力由弱到強依次為：幼果期<開花期<現(xiàn)蕾期[3]。春季冷空氣活動較為頻繁，加之海拔、環(huán)境等因素，可能有不同程度的凍害發(fā)生，尤其在高海拔種植區(qū)，低溫凍害對開花授粉或幼果生長有較大影響，會導(dǎo)致座果率及產(chǎn)量的下降[4]。在全球氣候變暖背景下，蘋果花期普遍提前，增加了遭遇低溫凍害的幾率和強度，導(dǎo)致蘋果生產(chǎn)風(fēng)險進一步加大[5-7]。

根據(jù)災(zāi)害風(fēng)險系統(tǒng)理論[8,9]與IPCC AR6第二工作組報告[10]，風(fēng)險的大小不僅取決于致災(zāi)因子發(fā)生的可能性及嚴重程度，也取決于近期的溫升水平、承災(zāi)體暴露度、脆弱性和社會經(jīng)濟發(fā)展水平及適應(yīng)措施等。目前的凍害研究主要側(cè)重于根據(jù)氣象指標劃分凍害等級，并且大多在陜西和山東開展，而結(jié)合承災(zāi)體脆弱性的河北蘋果低溫凍害風(fēng)險研究較少。在風(fēng)險研究方法上，近年來國內(nèi)不少成果采用基于指標體系的半定量方法，但是該方法的合理性存在一定爭議[11,12]。目前，探索定量化評價方法成為研究熱點，脆弱性曲線是通過建立不同致災(zāi)強度與承災(zāi)體損失的定量關(guān)系，反演可能的損失[13]。高歌等[14]以2016年9月漳州市山洪溝暴雨洪澇災(zāi)害為例，定量刻畫了淹沒水深與室內(nèi)財產(chǎn)損失的關(guān)系，驗證了該方法的可靠性。王丹丹等[15]構(gòu)建了不同過程累積降雪量與受災(zāi)人口、農(nóng)作物受災(zāi)面積之間的定量函數(shù)關(guān)系。

蘋果受凍既與低溫強度有關(guān)，也與持續(xù)時間密切相關(guān)，一般地，溫度越低，持續(xù)時間越長，受凍率越高，凍害越嚴重[16]。前人進行了氣候影響、凍害指標和風(fēng)險分析等相關(guān)研究。張山清等[17]分析了影響新疆蘋果生產(chǎn)的氣象要素變化特征，研究了氣候變化對蘋果種植氣候適宜性及分布的影響；李健等[18]以2 ℃、0 ℃、-2 ℃為不同凍害等級臨界氣溫，采用統(tǒng)計結(jié)合試驗等方法確定陜西蘋果花期凍害主要危害地區(qū)和時段；亦有國家氣象行業(yè)標準[19]通過日最低氣溫和低溫持續(xù)時間對富士系蘋果花期凍害等級進行劃分；屈振江等[20]以2013年4月我國一次典型凍害過程為基礎(chǔ)，通過最大熵模型建立了花期凍害與氣候因子的關(guān)系，提取主要致災(zāi)因子及閾值范圍，對蘋果主產(chǎn)區(qū)花期凍害風(fēng)險等級進行評估。以往研究中結(jié)合凍害指標的定量風(fēng)險研究報道較少。

本文基于災(zāi)情數(shù)據(jù)判斷了河北省蘋果開花—幼果期低溫凍害致災(zāi)閾值，并利用優(yōu)勢分析法和信息熵賦權(quán)法分別構(gòu)建了致災(zāi)強度指數(shù)和致災(zāi)危險性指數(shù)，分析了不同概率下的致災(zāi)強度與致災(zāi)因子危險性?；跒?zāi)害風(fēng)險系統(tǒng)理論，通過建立開花—幼果期低溫凍害致災(zāi)強度與蘋果受損的脆弱性曲線，開展蘋果低溫凍害風(fēng)險研究，以期為河北省蘋果低溫凍害風(fēng)險預(yù)警服務(wù)與評估業(yè)務(wù)提供技術(shù)支撐，對穩(wěn)定生產(chǎn)具有現(xiàn)實意義。

1 資料來源與處理

致災(zāi)因子數(shù)據(jù)為河北省氣象局1981—2020年142個國家氣象臺站氣溫資料，包括日最低氣溫和逐小時氣溫。根據(jù)河北省蘋果氣象中心規(guī)定，將河北省劃分為四個區(qū)域，分別為南部、中部、東部和北部（表1）。蘋果生育期數(shù)據(jù)來源于不同地區(qū)調(diào)查和觀測數(shù)據(jù)，確定花期和幼果期的平均物候期，并按四個區(qū)的相應(yīng)時段統(tǒng)計氣溫數(shù)據(jù)。

表1 河北省蘋果種植分區(qū)及花期、幼果期平均物候期

災(zāi)情數(shù)據(jù)為河北省各縣（市、區(qū)）低溫凍害災(zāi)情信息，包括開始時間、結(jié)束時間、發(fā)生地點、天氣過程描述、災(zāi)情描述等，來自文獻史料、民政部門、互聯(lián)網(wǎng)災(zāi)情和實地調(diào)查。對多源災(zāi)情進行匯總和除重，提取受凍率和減產(chǎn)率等損失信息，匹配當日當?shù)靥鞖庑畔?，用于承?zāi)體脆弱性曲線建立。

2 研究方法

2.1 致災(zāi)閾值判別

根據(jù)歷史區(qū)縣的災(zāi)情信息統(tǒng)計對應(yīng)區(qū)縣氣象臺站在凍害發(fā)生前及過程中的日最低氣溫構(gòu)建受災(zāi)、無災(zāi)樣本序列。如：災(zāi)情描述“2013年4月19-21日，低溫致使井陘縣6個鄉(xiāng)鎮(zhèn)的蘋果、梨及正值花期的其它林果和蔬菜遭受較大損失?！?，根據(jù)對應(yīng)區(qū)縣的氣象臺站查詢得到4月19-21日井陘縣日最低氣溫依次為0.3 ℃、0 ℃和3.6 ℃，確定受災(zāi)樣本為災(zāi)害發(fā)生時段的首日最低氣溫，為0.3 ℃，無災(zāi)樣本為前5日（4月14-18日）最低氣溫的最小值6.6 ℃，共反演得到樣本160對。每隔0.1 ℃進行判斷，計算總體準確度（Over Accuracy，OA）[21]，OA越大，對災(zāi)害的描述性越強。

式中：為總體準確度；為受災(zāi)和無災(zāi)樣本總數(shù)；為判斷正確的受災(zāi)樣本數(shù)，判識有災(zāi)實際有災(zāi)；為判斷錯誤的受災(zāi)樣本數(shù)，判識無災(zāi)而實際有災(zāi)；為判斷正確的無災(zāi)樣本數(shù)，判識無災(zāi)實際無災(zāi)；為判斷錯誤的無災(zāi)樣本數(shù)，判識有災(zāi)實際無災(zāi)。

采用感受性曲線（Receiver Operating Characteristic curve，ROC）[22]對蘋果開花—幼果期低溫致災(zāi)閾值進行評價。以假陽性率（）為橫坐標，真陽性率（）為縱坐標，ROC曲線越靠近左上角，準確性越高。以曲線覆蓋面積（Area Under Curve，AUC）作為評價統(tǒng)計量，范圍在0.5～1.0之間視為該判別方法有效，AUC越大，分類效果越好。和計算方法如下：

2.2 致災(zāi)強度指數(shù)

利用統(tǒng)計軟件SPSS，對單因子（最低氣溫）和兩因子（最低氣溫、持續(xù)時間）進行層級回歸，并分析因子與受凍率的相關(guān)性。受凍率為低溫凍害造成的蘋果受凍花朵占總花朵數(shù)的百分率[19]，通過災(zāi)情信息提取得到。基于日最低氣溫和低于閾值的持續(xù)小時數(shù)，構(gòu)建低溫凍害致災(zāi)強度指數(shù)（式4），通過優(yōu)勢分析法[23]判定兩因子的權(quán)重。優(yōu)勢分析方法是對全部致災(zāi)因子進行對比，分析相對重要程度，以致災(zāi)因子對致災(zāi)強度的平均貢獻之和占模型方差的百分比作為判別標準，不會夸大或忽視任一因子的重要性[24]。

=11+22(4)

式中：為致災(zāi)強度指數(shù)；1為歸一化的日最低氣溫；2為歸一化的持續(xù)小時數(shù)；1、2分別為最低氣溫與持續(xù)小時的權(quán)重。

提取1981-2020年四個種植分區(qū)相應(yīng)時段內(nèi)歷次低溫過程的日最低氣溫和低于致災(zāi)閾值的持續(xù)小時數(shù)，歸一化[25]后計算致災(zāi)強度指數(shù)，構(gòu)建年極值序列?；贕umbel、Weibull、Wakeby等多種分布函數(shù)，通過最大似然法[26]與KS檢驗[27]確定最優(yōu)分布函數(shù)，計算不同年遇水平的低溫凍害致災(zāi)強度指數(shù)。通過Kriging插值得到致災(zāi)強度與危險性分布，并利用自然斷點法[25]劃分等級。

網(wǎng)格模型如圖2(a)和圖2(b)所示，外界壓力為一個標準的大氣壓，圖示的氣壓入口可在Fluent中輸入相應(yīng)的氣壓參數(shù)，球狀可以以壁面(wall)處理。

2.2 致災(zāi)危險性

致災(zāi)強度指數(shù)可表征低溫凍害強度，當?shù)蜏仡l次越多時，發(fā)生凍害的可能性也越大，因此致災(zāi)危險性需要考慮強度和頻次的綜合作用。將頻次樣本累計值與致災(zāi)強度指數(shù)樣本平均值歸一化處理，采用信息熵賦權(quán)法[28,29]分別賦予權(quán)重，得到致災(zāi)危險性指數(shù)，公式如下：

=WX+WX(5)

式中：為致災(zāi)危險性指數(shù)；X為歸一化的致災(zāi)強度；X為歸一化的低溫頻次；W、W分別為致災(zāi)強度與頻次的權(quán)重。

2.3 承災(zāi)體脆弱性曲線

以脆弱性曲線來描述承災(zāi)體遭受災(zāi)害損失程度的變化規(guī)律，利用致災(zāi)強度指數(shù)與蘋果損失數(shù)據(jù)，分別使用線性、多項式、冪、指數(shù)和Logistic函數(shù)進行擬合，選擇最優(yōu)擬合函數(shù)建立致災(zāi)因子與災(zāi)害損失的對應(yīng)關(guān)系，繪制蘋果低溫凍害脆弱性曲線。

2.4 低溫凍害風(fēng)險

根據(jù)自然災(zāi)害風(fēng)險系統(tǒng)理論，風(fēng)險是由于致災(zāi)因子導(dǎo)致的承災(zāi)體可能的期望損失[9]，計算公式：

()=(())×(6)

式中：()為年一遇的損失風(fēng)險；()為年一遇的致災(zāi)強度；為承災(zāi)體的脆弱性函數(shù)；為承災(zāi)體暴露度。

3 結(jié)果與分析

3.1 致災(zāi)閾值判別

對構(gòu)建的160對受災(zāi)、無災(zāi)樣本序列進行統(tǒng)計分析，以0.1 ℃為間距，計算得到不同基準截點的總體準確度（表2）?？梢钥闯觯?1.0 ℃判斷正確的受災(zāi)樣本最多，為148個，但其判識正確的無災(zāi)樣本數(shù)量最少，其總體準確度僅有79.76%。-2.0 ℃～-1.1 ℃之間的總體準確度在80%以上，其中，-1.4 ℃的總體準確度最高，為83.95%，其判斷正確的受災(zāi)樣本數(shù)（TP）為144個，判斷正確的無災(zāi)樣本數(shù)（TN）為119個。其次為-1.6 ℃，總體準確度為83.40%。

表2 正誤樣本數(shù)量及總體準確度

建立ROC感受性曲線，對基于日最低氣溫的蘋果開花—幼果期低溫凍害閾值進行評價。由圖1可見，建立的曲線上的各點均處于(0,0)和(1,1)連線的上方，經(jīng)計算，該曲線的評價統(tǒng)計量AUC值為0.87，表明本方法整體性能較好，若妥善設(shè)定閾值，具有判別價值。-1.4 ℃為最佳的臨界閾值，其對應(yīng)的真陽性率（）為0.918，假陽性率（）為0.240。

圖1 ROC曲線

3.2 致災(zāi)危險性

基于SPSS進行層級回歸分析，僅考慮最低氣溫一個因子時（表3），最低氣溫與受凍率的相關(guān)系數(shù)為0.614（<0.01）。將持續(xù)小時因子納入，2改變值為0.138，方差分析顯著，具有統(tǒng)計學(xué)意義，說明排除最低氣溫的影響，持續(xù)小時也是一個重要影響因素。

表3 層級回歸分析結(jié)果

注：模型1因子為最低氣溫；模型2因子為最低氣溫和持續(xù)小時。

Note: Factor of model 1 is the minimum temperature; The factors of model 2 are minimum temperature and duration hours.

通過優(yōu)勢分析法構(gòu)建低溫凍害致災(zāi)強度指數(shù)（表4）。當各因子被加入到模型后，帶來的增值貢獻為Δ2。以受凍率為因變量，當模型中只有最低氣溫因子時，判定系數(shù)2為0.377；當加入持續(xù)小時因子時，2增加0.138，為0.515。以總平均貢獻表示每個因子的重要程度，可以看出，最低氣溫和持續(xù)小時在對受凍率的影響權(quán)重分別為76.9%和23.1%。

表4 各因子的增值貢獻及總平均貢獻

圖2為不同年遇水平下的蘋果開花—幼果期低溫凍害致災(zāi)強度分布。致災(zāi)強度自西北向東南逐漸降低，10年一遇下，較低強度分布在張家口中部和西南部、承德中部和保定西北部，中等強度及以上分布在張家口北部和承德西北部，其中，張家口北部處于較高強度，無高強度。30年一遇，較高強度的范圍有所擴大，分布在張家口中北部和承德西北部，50年一遇時，張家口北部康保縣處于高強度等級，保定西北部達到中等強度等級。

圖2 不同年遇水平的低溫凍害致災(zāi)強度

信息熵賦權(quán)法確定致災(zāi)強度和頻次的權(quán)重分別為0.624和0.376，計算致災(zāi)危險性指數(shù)，得到蘋果開花—幼果期低溫凍害致災(zāi)危險性分布（圖3）。西北部的致災(zāi)危險性高于東南部，張家口中北部和承德西北部處于中等危險性等級，張家口西北部為高危險性，全省大部分地區(qū)處于低危險性和較低危險性等級。

3.3 承災(zāi)體暴露度和脆弱性

2017年河北省蘋果種植面積為24.53萬hm2，結(jié)果率為77.3%，品種以紅富士和國光為主，紅富士占比52.6%，國光占比13.2%。承德市的蘋果種植面積最大，為5.3728萬hm2。東北部產(chǎn)量較高，其次是南部，承德圍場產(chǎn)量最高，達30.30萬t，產(chǎn)量在2.0萬噸以上的縣區(qū)集中在河北省中南部。圖4為河北省蘋果單產(chǎn)分布。全省平均單產(chǎn)為17.78 t/hm2。南部和東部的單產(chǎn)水平較高。單產(chǎn)在15～25 t/hm2的有48個縣區(qū)，25～40 t/hm2的有35個縣區(qū)，40～80 t/hm2的有5個縣區(qū)，分別是任縣、定州、隆堯、樂亭、運河區(qū)。衡水阜城縣的單產(chǎn)水平最高，為124.10 t/hm2。

圖3 蘋果開花—幼果期低溫凍害致災(zāi)危險性分布

圖4 2017年河北省各縣蘋果單產(chǎn)分布

以致災(zāi)強度為軸，以減產(chǎn)率為軸，采用線性、指數(shù)、多項式、冪和Logistic函數(shù)進行擬合，采用最優(yōu)模型構(gòu)建蘋果低溫凍害脆弱性曲線，見圖5。通過擬合優(yōu)度檢驗（表5），指數(shù)函數(shù)的確定系數(shù)2最接近1，其剩余平方和（RSS）和均方根誤差（RMSE）均為最小，擬合度最佳，因此選取指數(shù)函數(shù)來表達由低溫導(dǎo)致的蘋果減產(chǎn)率脆弱性。

圖5 河北省蘋果低溫凍害脆弱性曲線

表5 曲線方程及擬合優(yōu)度

3.4 低溫凍害風(fēng)險

圖6為不同致災(zāi)強度水平下河北省蘋果受低溫凍害的受影響產(chǎn)量分布，其數(shù)值越高，風(fēng)險越大。各個低溫凍害致災(zāi)強度下蘋果受影響產(chǎn)量的空間分布趨勢大體一致，與暴露水平密切相關(guān)。10年一遇低溫凍害致災(zāi)強度水平下，南部蘋果受影響產(chǎn)量主要在2.1～4.0 t/hm2，承德最北部圍場縣的東南地區(qū)、唐山東部遷安和灤縣以及南部樂亭和曹妃甸、秦皇島西部盧龍、保定東北部定興和東南部地區(qū)、定州市以及河北省南部的部分地區(qū)，受影響產(chǎn)量在4.1～8.0 t/hm2；圍場中西部和邯鄲市中部冀南新區(qū)在8.1～12.0 t/hm2；衡水阜城縣風(fēng)險最大，受影響產(chǎn)量為16.1 t/hm2。30年一遇致災(zāi)強度下，大部分蘋果種植區(qū)的受影響產(chǎn)量增至2.1～4.0 t/hm2，圍場西部在12.1 t/hm2以上，甚至高達16.9 t/hm2。50年一遇，受影響產(chǎn)量和范圍進一步擴大，其中，圍場自西北向東南呈逐級遞減，受影響產(chǎn)量最大為19.9 t/hm2，其次為阜城縣。

圖6 河北省蘋果低溫凍害風(fēng)險分布

4 結(jié)論與討論

（1）基于災(zāi)情信息確定的河北省蘋果開花—幼果期低溫致災(zāi)閾值為-1.4 ℃，總體準確度為83.95%。

（2）基于優(yōu)勢分析法確定最低氣溫和持續(xù)小時對受凍率的影響權(quán)重分別為76.9%和23.1%。致災(zāi)強度與致災(zāi)因子危險性的空間分布大體一致，自西北向東南逐漸降低。

（3）蘋果低溫凍害風(fēng)險高值區(qū)主要分布在暴露度水平較高的地區(qū)。張家口東部和承德中西部的致災(zāi)危險性較高，但單產(chǎn)水平較低，損失風(fēng)險較低。10年一遇致災(zāi)強度下，南部蘋果受影響產(chǎn)量在2.1～4.0 t/hm2，衡水阜城風(fēng)險最大，為16.1 t/hm2。受致災(zāi)因子和承災(zāi)體綜合作用，圍場西部風(fēng)險自西北向東南逐級遞減，50年一遇下，受影響產(chǎn)量為全省最大。

從致災(zāi)因子角度來看，燕山以北的河北省東北部的致災(zāi)強度和危險性處于中等級及以上，尤其張家口東北部處于低溫凍害高和較高危險性，因此蘋果不宜種植在海拔1200～1500 m的壩上地區(qū)。許俊卿等[30]根據(jù)適宜指標對比，指出承德隆化、興隆以及唐山灤州、樂亭等地區(qū)比較適宜種植，均分布在致災(zāi)危險性的低等級和較低等級地區(qū)。對于蘋果種植較多的承德和張家口的部分縣區(qū)，可考慮多栽種抗逆性強的品種或增加樹體營養(yǎng)，以提高生命活性和抗性。

蘋果低溫凍害高風(fēng)險主要分布在單產(chǎn)水平較高的地區(qū)，這與暴露水平密切相關(guān)。對于張家口東部和承德中西部，蘋果種植面積雖大，但產(chǎn)量和種植密度相對較低，單產(chǎn)及暴露度較低，即使其低溫凍害致災(zāi)強度及危險性較高，受低溫凍害影響的風(fēng)險也會較低。隨著致災(zāi)強度年遇水平的增加，承德圍場的受影響產(chǎn)量超過衡水阜城，成為全省最高，這種時空分異特征是致災(zāi)因子危險性、承災(zāi)體暴露度和脆弱性共同作用的結(jié)果。

本文確定的河北省蘋果開花—幼果期低溫凍害致災(zāi)閾值為-1.4 ℃，符合多位學(xué)者[31-33]給出的開花期受凍臨界低溫為-2.0 ℃，幼果受凍臨界低溫為-1.0 ℃和-2 ℃～0 ℃為輕度凍害等結(jié)論。通過氣溫特征分析，南部開花期和中部開花—幼果期有多個年份的日最低氣溫低于-1.4 ℃，花期的氣溫相比幼果期較低，降溫頻次較多，其危險性相對較高，更需關(guān)注寒潮天氣和霜凍災(zāi)害的預(yù)報預(yù)警信息，及時采取灌水、噴防凍劑、撒石灰等防御措施，最大程度減輕不利影響。

本研究基于風(fēng)險系統(tǒng)理論，通過建立脆弱性曲線，對河北省開花—幼果期的蘋果低溫凍害的定量風(fēng)險進行了初步探索，在實際業(yè)務(wù)中，可與氣溫預(yù)報相結(jié)合，提供區(qū)域性的風(fēng)險預(yù)警服務(wù)。脆弱性分析和風(fēng)險評價結(jié)果可與種植區(qū)特色農(nóng)業(yè)氣象災(zāi)害保險結(jié)合，用于建立氣象指數(shù)保險或政策性保險制度，科學(xué)厘定費率，化解蘋果種植風(fēng)險。今后在構(gòu)建蘋果低溫凍害脆弱性曲線時，可以針對不同品種的蘋果結(jié)合實驗?zāi)M來進行脆弱性研究。應(yīng)進一步考慮地形地貌導(dǎo)致的孕災(zāi)環(huán)境差異，挖掘并積累更多災(zāi)情，針對四個蘋果種植分區(qū)來考慮脆弱性，以完善評價模型。

[1] 王景紅,柏秦鳳,梁軼,等.2013年陜西蘋果花期凍害氣象條件分析及受凍指標研究[J].果樹學(xué)報,2015,32(1):100-107

[2] 楊小利.甘肅平?jīng)鍪刑O果花期凍害風(fēng)險分析及區(qū)劃[J].干旱地區(qū)農(nóng)業(yè)研究,2014,32(4):231-235

[3] 劉英勝.果樹花期凍害發(fā)生原因與預(yù)防技術(shù)[J].河北果樹,2010(3):31-32

[4] 孫小諾,許俊卿,王蓉蓉,等.河北蘋果花期低溫冷害預(yù)報模型研究[J].陜西農(nóng)業(yè)科學(xué),2020,66(1):71-74

[5] Diffenbaugh NS, Burke M. Global warming has increased global economic inequality [J]. Preceedings of the National Academy of Sciences, 2019,116(20):9808-9813

[6] Augspurger CK. Reconstructing patterns of temperature, phenology, and frost damage over 124 years: Spring damage risk is increasing [J]. Ecology, 2013,94(1):41-50

[7] 楊小利,江廣勝.隴東黃土高原典型站蘋果生長對氣候變化的響應(yīng)[J].中國農(nóng)業(yè)氣象,2010,31(1):74-77

[8] 秦大河,張建云,閃淳昌,等.中國極端天氣氣候事件和災(zāi)害風(fēng)險管理與適應(yīng)國家評估報告[M].北京:科學(xué)出版社,2015

[9] 史培軍.中國自然災(zāi)害風(fēng)險地圖集[M].北京:科學(xué)出版社,2011

[10] IPCC. Climate change 2022: impacts, adaptation, and vulnerability [M/OL]. Cambridge: Cambridge University Press, 2022[2022-03-31]. http://www.ipcc.ch/report/ar6/wg2/

[11] 章國材.自然災(zāi)害風(fēng)險評估與區(qū)劃原理和方法[M].北京:氣象出版社,2014

[12] 李婷,陳笑娟,張靜,等.基于脆弱性曲線的河北省日光溫室大風(fēng)災(zāi)害風(fēng)險評估[J].湖北農(nóng)業(yè)科學(xué),2020,59(4):63-69

[13] 李婷,孫玉龍,陳笑娟,等.基于經(jīng)濟損失的河北省雪災(zāi)風(fēng)險定量評估[J].災(zāi)害學(xué),2018,33(4):72-77

[14] 高歌,王長源,吳先華,等.山洪溝暴雨洪澇災(zāi)害情境下居民室內(nèi)財產(chǎn)的脆弱性曲線研究——以2016年9月福建省長 泰縣馬洋溪山洪災(zāi)害為例[J].災(zāi)害學(xué),2018,33(3):126-131

[15] 王丹丹,潘東華,賈慧聰,等.基于歷史災(zāi)情的貴州雪災(zāi)脆弱性分析[J].災(zāi)害學(xué),2018,33(3):24-28

[16] 王景紅,劉璐,高峰,等.陜西富士系蘋果花期霜凍災(zāi)害氣象指標的修訂[J].中國農(nóng)業(yè)氣象,2015,36(1):50-56

[17] 張山清,普宗朝,李新建,等.氣候變化對新疆蘋果種植氣候適宜性的影響[J].中國農(nóng)業(yè)資源與區(qū)劃,2018,39(8):255-264

[18] 李健,劉映寧,李美榮,等.陜西果樹花期低溫凍害特征及防御對策[J].氣象科技,2008,36(3):318-322

[19] 中國氣象局.QX/T392-2017富士系蘋果花期凍害等級[S].北京:中國氣象出版社,2017

[20] 屈振江,周廣勝,魏欽平.蘋果花期凍害氣象指標和風(fēng)險評估[J].應(yīng)用氣象學(xué)報,2016,27(4):385-395

[21] 周欣昕,吳艷蘭,李夢雅,等.基于特征分離機制的深度學(xué)習(xí)植被自動提取方法[J].地球信息科學(xué)學(xué)報,2021,23(9):1675-1689

[22] 徐敏,崔玉龍,袁文華,等.基于GIS和信息量法的伊犁黃土滑坡易發(fā)性評價[J].河南城建學(xué)院學(xué)報,2021,30(4):25-31

[23] Azen R, Budescu DV. The dominance analysis approach for comparing predictors in multiple regression [J]. Psychological methods, 2003,8(2):129-148

[24] 孫玉龍,王麗榮,張琪,等.基于優(yōu)勢分析法的河北省暴雨對房屋致災(zāi)閾值研究[J].氣候與環(huán)境研究,2021,26(3):289-298

[25] 中國氣象局.QX/T405-2017雷電災(zāi)害風(fēng)險區(qū)劃技術(shù)指南[S].北京:中國氣象出版社,2017

[26] 榮艷淑,王文,王鵬,等.淮河流域極端降水特征及不同重現(xiàn)期降水量估計[J].河海大學(xué)學(xué)報(自然科學(xué)版),2012,40(1):1-8

[27] 許瀚卿,譚金凱,李夢雅,等.基于Copula函數(shù)的沿海城市雨潮復(fù)合災(zāi)害風(fēng)險研究[J].自然災(zāi)害學(xué)報,2022,31(1):40-48

[28] 郭昱.權(quán)重確定方法綜述[J].農(nóng)村經(jīng)濟與科技,2018,29(8):252-253

[29] 孫魯云,王力.基于層次分析法-熵權(quán)法的中國棉花質(zhì)量綜合評價與分析[J].江蘇農(nóng)業(yè)學(xué)報,2022,38(3):642-649

[30] 許俊卿,孫小諾,王蓉蓉,等.河北省蘋果主要種植縣氣象條件分析[J].陜西農(nóng)業(yè)科學(xué),2019,65(9):90-93

[31] 李美榮,李星敏,柏秦鳳,等.蘋果極端氣象災(zāi)害氣溫極值的分布及重現(xiàn)期預(yù)測[J].干旱地區(qū)農(nóng)業(yè)研究,2012,30(3):257-261

[32] 汪景彥.紅富士蘋果高產(chǎn)栽培[M].北京:金盾出版社,1993

[33] 柏秦鳳,霍治國,王景紅,等.中國主要果樹氣象災(zāi)害指標研究進展[J].果樹學(xué)報,2019,36(9):1229-1243

Risk Assessment on the Freezing Injury of Apple in Hebei Province Based on the Fragility Curve

LI Ting1, SUN Yu-long1, CHEN Xiao-juan1,2*, XU Yi-fu1

1.,050021,2.,050021,

Based on the temperature data of 142 stations in Hebei Province from 1981 to 2020, freezing injury intensity and hazard danger index during the apple florescence-young fruit period were constructed by Dominance analysis and Entropy method. The injury intensity of different probabilities was analyzed. Based on the risk system theory, ‘intensity index - reduction rate of yield’ curves were established to study the freezing injury risk. The results showed that the overall accuracy of temperature between -2.0°C and -1.1°C was above 80%, and the threshold was determined to be -1.4℃. The injury intensity and hazard gradually decreased from northwest to southeast, and the high hazard danger was distributed in the northwest of Zhangjiakou, while other areas had low and lower levels. The high-risk was mainly concentrated in the counties with higher unit yield. In 10-year return period, the risk of Fucheng County was greatest, with an affected yield of 16.1 t/hm2. In 30-year return period, the affected yield in most planting areas increased to 2.1～4.0 t/hm2. In 50-year return period, Weichang County showed a gradual decline from northwest to southeast, and the highest affected yield was 19.9 t/hm2. Due to the combined effect of hazard, exposure and disaster-bearing body vulnerability, the risk of Weichang was higher. The risk of eastern Zhangjiakou and center-western Chengde was lower because of the exposure. In the areas with more planting in Chengde and Zhangjiakou, strong resistance varieties should be considered.

Apple; freezing injury; risk assessment

S661.1/X820.4

A

1000-2324(2022)05-0742-08

2022-02-13

2022-03-25

河北省氣象局科研開發(fā)項目(19ky24);河北省氣象局省級財政投資業(yè)務(wù)建設(shè)項目(202201-01)

李婷(1989-),女,碩士,工程師,從事氣象災(zāi)害風(fēng)險評估工作. E-mail:liting29320@sina.com

通訊作者:Author for correspondence. E-mail:juan2044335@163.com