王鶯，王健順，張強

（1. 中國氣象局蘭州干旱氣象研究所，甘肅省干旱氣候變化與減災(zāi)重點實驗室，中國氣象局干旱氣候變化與減災(zāi)重點開放實驗室，甘肅蘭州 730020；2. 蘭州大學(xué)大氣科學(xué)學(xué)院，甘肅 蘭州 730020）

草地通常指以草本植物占優(yōu)勢的植物群落［1］。作為重要的生態(tài)系統(tǒng)類型，草地在防風(fēng)、固沙、保土、調(diào)節(jié)氣候、凈化空氣、涵養(yǎng)水源等方面發(fā)揮著重要的生態(tài)功能。同時，草地資源也是草地農(nóng)業(yè)的重要基礎(chǔ)，提供飼草飼料支撐著畜牧業(yè)的生產(chǎn)，在改善居民膳食結(jié)構(gòu)和營養(yǎng)水平、促進農(nóng)業(yè)農(nóng)村經(jīng)濟發(fā)展和農(nóng)牧民增收方面發(fā)揮著重要作用。

中國的天然草原面積約4 億hm2，占國土面積的42%，總面積位居世界第二。中國也是畜牧業(yè)大國之一，截至2019 年，大牲畜飼養(yǎng)規(guī)模達(dá)到9877.42 萬頭。中國的草原大約有2/3 位于干旱、半干旱地區(qū)，主要依靠雨水涵養(yǎng)，是陸地生態(tài)系統(tǒng)中最為脆弱的開放性系統(tǒng)，干旱災(zāi)害成為制約中國草原生產(chǎn)力的最重要因素［2］。20 世紀(jì)70 年代以來，由于大氣環(huán)流系統(tǒng)從對流層到平流層發(fā)生了明顯的年代際轉(zhuǎn)折，中國北方干旱化趨勢增強，大范圍干旱連年發(fā)生，與此同時，南方干旱也出現(xiàn)了明顯的增加和加重趨勢。中國農(nóng)作物干旱年均受災(zāi)面積和成災(zāi)面積分別達(dá)2.1×107和8.9×106hm2，干旱造成的直接經(jīng)濟損失高達(dá)440 億元·年-1［3］。氣候變化背景下，中國未來年平均氣溫上升，整體將呈干旱化趨勢［4-5］。這將導(dǎo)致更多的干旱事件發(fā)生，并可能對草地生態(tài)系統(tǒng)造成不可逆的影響，從而為中國畜牧業(yè)的穩(wěn)定發(fā)展帶來極大挑戰(zhàn)。因此，如何進行草原干旱災(zāi)害的防災(zāi)減災(zāi)工作就成了社會關(guān)注的熱點問題。

隨著從被動抗旱向主動防旱的防災(zāi)減災(zāi)新理念的提出，干旱災(zāi)害由危機管理向風(fēng)險管理轉(zhuǎn)變，而干旱災(zāi)害風(fēng)險評估就是風(fēng)險管理的核心內(nèi)容［6］。IPCC（the intergovernmental panel on climate change）認(rèn)為災(zāi)害風(fēng)險是脆弱性、危險性和暴露性共同作用的結(jié)果［7］。張繼權(quán)等［8］認(rèn)為干旱災(zāi)害風(fēng)險系統(tǒng)是由致災(zāi)因子危險性和旱災(zāi)脆弱性（承災(zāi)體暴露性、災(zāi)損敏感性和抗旱能力）相互聯(lián)系、相互作用下形成的復(fù)雜系統(tǒng)。目前，對于草原干旱的研究主要是基于干旱指數(shù)來定量描述干旱程度，其中包括氣象干旱指數(shù)，例如降水距平百分率、標(biāo)準(zhǔn)化降水指數(shù)、標(biāo)準(zhǔn)化降水蒸散指數(shù)、帕默爾干旱指數(shù)等［9-10］，以及遙感干旱指數(shù)，例如植被供水指數(shù)（vegetation supply water index，VSWI）、植被狀態(tài)指數(shù)（vegetation condition index，VCI）、溫度植被干旱指數(shù)（temperature vegetation dryness index，TVDI）等［11-12］。這些成果屬于干旱致災(zāi)因子研究，是干旱災(zāi)害形成的直接影響因素，對于認(rèn)識草原干旱的起止時間、發(fā)生頻率、出現(xiàn)強度、受災(zāi)面積等具有重要意義。作為一個復(fù)雜系統(tǒng)，干旱災(zāi)害風(fēng)險研究還需要充分考慮暴露度和脆弱性在干旱成災(zāi)中的作用。關(guān)于脆弱性，研究人員主要基于敏感性和適應(yīng)能力，構(gòu)建評價指標(biāo)體系，分析脆弱性特征［13-14］。草原干旱災(zāi)害風(fēng)險是干旱危險性、暴露性和脆弱性的耦合。李文龍等［15］基于標(biāo)準(zhǔn)化降水蒸散指數(shù)（standardized precipitation evapotranspiration index，SPEI）與植被指數(shù)，評價了高寒草地的干旱生態(tài)風(fēng)險。張存厚等［16］通過建立干旱風(fēng)險評價指標(biāo)和評價模型，評價了內(nèi)蒙古草原干旱災(zāi)害綜合風(fēng)險。

從以上分析可以看出，目前對于草原干旱風(fēng)險的研究大多著眼于致災(zāi)因子，而缺少了對承災(zāi)體和孕災(zāi)環(huán)境構(gòu)成的復(fù)雜系統(tǒng)的影響評估研究，特別是針對中國草原的干旱風(fēng)險評估研究還較薄弱?；谝陨显?，本研究通過干旱風(fēng)險分析，從致災(zāi)因子危險性、承災(zāi)體暴露性和敏感性，以及孕災(zāi)環(huán)境脆弱性4 個方面構(gòu)建了草原干旱風(fēng)險評價指標(biāo)框架，對中國草原干旱災(zāi)害進行風(fēng)險評估和區(qū)劃，以期為實現(xiàn)草原的科學(xué)抗旱和主動防旱提供基礎(chǔ)數(shù)據(jù)和理論支撐。

1 材料與方法

1.1 研究區(qū)概況

中國位于亞歐大陸東部，太平洋西部，海拔落差大，氣候類型復(fù)雜多樣。西北部深居內(nèi)陸，海拔較高，氣候干旱；東南部瀕臨海洋，海拔較低，氣候溫濕。草原是中國面積最大的陸地生態(tài)系統(tǒng)，約占國土面積的42%，其空間分布是氣候與地貌等長期共同作用的結(jié)果。主要分為草甸草原、典型草原、荒漠草原、高寒草原和草叢等。草甸草原主要分布在松遼平原和內(nèi)蒙古高原的東部邊緣。典型草原主要分布在內(nèi)蒙古、東北西南部、黃土高原中西部和阿爾泰山、天山以及祁連山的部分區(qū)域?；哪菰饕植荚趦?nèi)蒙古中部、黃土高原北部以及祁連山和天山的低山帶。高寒草原主要分布在青藏高原、帕米爾高原及祁連山和天山的高海拔處。草叢主要分布在森林退化后的區(qū)域。

1.2 資料及其來源

氣象數(shù)據(jù)主要包括1978-2020 年全國836 個氣象站的日降水量、氣溫、風(fēng)速、日照時數(shù)、相對濕度等觀測資料，資料來源于中國氣象局的國家氣象科學(xué)數(shù)據(jù)共享服務(wù)平臺（http：//data.cma.cn/）。該數(shù)據(jù)經(jīng)過嚴(yán)格的質(zhì)量控制［17］。

遙感資料數(shù)據(jù)來自美國宇航局中分辨率成像光譜輻射計（national aeronautics and space administration/moderate resolution imaging spectroradiometer，NASA/MODIS）陸地產(chǎn)品組按照統(tǒng)一算法開發(fā)的2000-2020 年的月最大合成植被指數(shù)產(chǎn)品MODIS13A3（1000 m×1000 m）和2000-2020 年的年最大合成凈初級生產(chǎn)力（net primary productivity，NPP）產(chǎn)品MODIS17A3（1000 m×1000 m），下載網(wǎng)址為https：//search. earthdata. nasa.gov/search。該資料被廣泛應(yīng)用于中國的植被研究［18］。

中國草地類型數(shù)據(jù)來自1996 年出版的1∶100 萬草地資源圖集［19］。該數(shù)據(jù)根據(jù)植被-生境多因素分類法，將中國天然草地分為18 類。土壤田間持水量和土壤萎蔫系數(shù)空間數(shù)據(jù)是由中國科學(xué)院南京土壤研究所根據(jù)1994年完成的第二次中國土壤調(diào)查數(shù)據(jù)制作的1∶100 萬中國土壤空間數(shù)據(jù)庫提供。

1.3 研究方法

通過對中國草原干旱災(zāi)害風(fēng)險的分析與識別，從危險性、暴露性、敏感性和脆弱性入手，構(gòu)建草原干旱災(zāi)害風(fēng)險評價指標(biāo)體系和風(fēng)險評價模型，獲得中國草原干旱災(zāi)害風(fēng)險空間分布特征。具體技術(shù)路線見圖1。

圖1 中國草原干旱災(zāi)害風(fēng)險評價技術(shù)路線Fig.1 Technical route of grassland drought disaster risk assessment in China

1.3.1干旱災(zāi)害致災(zāi)因子危險性 就目前各類干旱致災(zāi)指標(biāo)來看，SPEI 能夠較好地表征氣象干旱致災(zāi)程度，其負(fù)值越大意味著氣象條件越干旱。SPEI 是由Vicente-Serrano 等［20］首次提出的，基于水分平衡原理，通過對降水量和蒸散量之間的差值序列累積概率值進行正態(tài)標(biāo)準(zhǔn)化構(gòu)成的指數(shù)。它用降水量和潛在蒸散量差值的平均偏離程度來綜合表征局地干濕狀況。SPEI 具備多時間尺度特征，有助于識別不同的干旱類型和干旱對不同系統(tǒng)的影響，同時也綜合考慮了溫度和降水對干旱的敏感性，適用于干旱災(zāi)害風(fēng)險研究［21-23］。SPEI 的計算步驟見文獻(xiàn)［20］。

已有研究［15，24］發(fā)現(xiàn)，6 個月時間尺度的SPEI（SPEI6）能夠較好地表征草地水分的虧缺或盈余，常用于草地干旱問題研究，所以這里用SPEI6 作為草原干旱災(zāi)害的致災(zāi)因子指標(biāo)?？筛鶕?jù)SPEI 的正態(tài)分布曲線，獲得干旱等級劃分閾值［25］（表1）。

從理論上講，致災(zāi)因子危險性一般由致災(zāi)因子強度與它的出現(xiàn)概率共同決定。這里將草原干旱致災(zāi)因子危險性定義為：

式中：H是干旱致災(zāi)因子危險性，單位為%；IC是干旱致災(zāi)因子強度等級值，為無因次數(shù)；ICx是干旱致災(zāi)因子特強等級值，為無因次值；i是干旱致災(zāi)因子等級（表1）；PCi是等級i的干旱致災(zāi)強度的出現(xiàn)概率，單位為%，可用信息分配法計算得到。

表1 標(biāo)準(zhǔn)化降水蒸散指數(shù)旱澇等級劃分Table 1 Classification standard of drought based on SPEI

信息分配法是在傳統(tǒng)直方圖法計算樣本概率基礎(chǔ)上的優(yōu)化方法，即將直方圖邊界模糊化處理，使得過渡邊界的信息也得以利用，從而增加小樣本概率估算的合理性［26］，其具體計算流程如下：

首先，令X為研究區(qū)在過去n年內(nèi)的實際觀測值樣本集合，在本研究中X可代表SPEI 樣本集合，即：

式 中：x1、x2、x3……xm是SPEI 的樣本。然后，根 據(jù) 樣本集合X中的最大值和最小值，選取適當(dāng)?shù)牟介L（△），得到觀測樣本的控制點空間（U）：

式中：u1、u2、u3……un是控制點樣本。然后，可以將每個單值觀測樣本xi所攜帶的信息擴散給U中的所有點：

式中：h是信息擴散系數(shù)，其解析表達(dá)式如下［27］：

相應(yīng)的模糊子集的隸屬函數(shù)為：

可把式（7）中μxi(uj)稱為樣本xi的歸一化信息分布。令：

式中：Q為控制點uj獲得的信息總量?？傻茫?/p>

式中：p(uj)為樣本落在uj處的頻率值，可作為概率的估計值。

1.3.2承災(zāi)體暴露性 承災(zāi)體暴露性反映了系統(tǒng)暴露在致災(zāi)因子影響下可能遭受損失的價值或數(shù)量。一個區(qū)域暴露于干旱危險因素的價值密度越高，可能遭受的潛在損失越大。NPP 是植被在單位時間和單位面積由光合作用所產(chǎn)生的有機物總量中扣除自養(yǎng)呼吸后的剩余部分，可以代表中國草原的生產(chǎn)能力［28-29］。因此選擇NPP作為承災(zāi)體暴露性指標(biāo)。

1.3.3承災(zāi)體敏感性 承災(zāi)體敏感性指承災(zāi)體對于環(huán)境變化的敏感程度。在本研究中選擇增強型植被指數(shù)（enhanced vegetation index，EVI）變異系數(shù)和草地類型來表示。

植被指數(shù)是反映草地長勢和營養(yǎng)信息的重要參數(shù)。很多研究人員通過構(gòu)建植被指數(shù)與草地生物量的關(guān)系模型，進而去評估草地生產(chǎn)力［30-31］。其中EVI 通過加入藍(lán)色波段來增強植被信號，進而矯正土壤背景和氣溶膠散射對植被的影響，避免了歸一化植被指數(shù)（normalized difference vegetation index，NDVI）在植被高生物量區(qū)易飽和的缺點［32］。因此選擇EVI 作為評價草地生長狀況的評價因子。草地EVI 的變異系數(shù)指某一區(qū)域EVI 的波動大小，體現(xiàn)了草地對外來干擾的敏感性。EVI 波動大的地區(qū)草地的干旱敏感性較高，EVI 波動小的地區(qū)草地的干旱敏感性則較低。具體計算公式如下：

式中：CVEVI是EVI變異系數(shù)；SDEVI是EVI的標(biāo)準(zhǔn)偏差是EVI的多年平均值；i是統(tǒng)計單元。根據(jù)植被-生境學(xué)分類法，將具有相同水熱氣候帶特征和植被特征，具有獨特地帶性，或具有廣域性分布的隱域性特征的草地分為一類，可將中國天然草地分為18 類。溫性草甸草原類主要位于溫帶半濕潤地區(qū)，中旱生多年生叢生禾草及根莖禾草占優(yōu)勢，?；煊写罅恐猩蚝抵猩参?。溫性草原類主要分布在溫帶半干旱氣候帶，旱生多年生叢生禾草占優(yōu)勢，混生一定數(shù)量的旱生小灌木。溫性荒漠草原類發(fā)育在溫帶干旱地區(qū)，強旱生多年生草本植物占優(yōu)勢，也有部分旱生小半灌木。高寒草原類主要分布在高山和青藏高原寒冷干旱的氣候區(qū)，植被以寒旱生多年生禾草為主，混生墊狀、匍匐狀植物和高寒藻類。溫性荒漠類分布于溫帶極干旱與嚴(yán)重缺水的生境條件下，植被以超旱生小灌木和小半灌木占優(yōu)勢。高寒草甸類主要分布在高原（高山）亞寒帶和寒帶地區(qū)寒冷濕潤的氣候條件下，植被以耐寒多年生中生草本植被為主。灌草叢類主要是由于森林、灌叢被砍伐，生境趨于干旱化所形成的次生類型，以中生或旱中生多年生草本為建群種。由于草原植被的分布對長期干旱分布具有一定的選擇性，因此根據(jù)草原不同類的空間分布特征和優(yōu)勢種的抗旱特性，結(jié)合牧草干鮮比［33］，通過專家打分法為18 類天然草地和改良草地進行干旱敏感性賦值（表2）。值越高，干旱敏感性越低。

表2 不同草原類型的干旱敏感性Table 2 Drought sensitivity of different grassland types

承災(zāi)體敏感性評估模型為：

式中：S為承災(zāi)體敏感性；XSi為第i種指標(biāo)的敏感性指數(shù)；WSi為第i種指標(biāo)的敏感性權(quán)重；n為指標(biāo)個數(shù)。

1.3.4孕災(zāi)環(huán)境脆弱性 孕災(zāi)環(huán)境脆弱性表征了環(huán)境條件對災(zāi)害傳遞的脆弱性程度。主要包括土壤環(huán)境和地形地貌等自然因素。

研究區(qū)的土壤環(huán)境用0～10 cm、10～20 cm 和20～30 cm 剖面的土壤有效水含量來表示。土壤有效水含量是田間持水量與土壤凋萎系數(shù)之間的土壤含水量，表征了土壤中能被草地吸收利用的水量。

研究區(qū)的地形地貌用坡度和坡向來表示。土壤水分累積入滲率隨坡度的增加而減少［34］。坡度小的區(qū)域，土壤容重小，單位面積降水量多，降水入滲快，地表徑流少，土壤的保水能力強，有利于草地生物量的形成［35］。而大坡度還會造成草地根系表面化，影響植物對水分的吸收利用。將坡度分為0°～5°、5°～10°、10°～15°、15°～20°、20°～25°、25°～40°和＞40°共7 類，構(gòu)建坡度指數(shù)。坡向?qū)τ谌照諘r數(shù)和太陽輻射強度有較大的影響。對于北半球來說，陽坡的輻射收入最多，其次為半陽坡，再次為半陰坡，最少為陰坡。將坡向分為陰坡（315°～45°）、陽坡（135°～225°）、半陰坡（90°～135°，225°～270°）、半陽坡（45°～90°，270°～315°）和平地（0°）共5 類，構(gòu)建坡向指數(shù)。

孕災(zāi)環(huán)境脆弱性評估模型為：

式中：V為脆弱性；XVi為第i種指標(biāo)的脆弱性指數(shù)；WVi為第i種指標(biāo)的脆弱性權(quán)重；n為指標(biāo)個數(shù)。

1.3.5干旱風(fēng)險評價指標(biāo)標(biāo)準(zhǔn)化處理 草原干旱風(fēng)險評價指標(biāo)體系中涉及的多種指標(biāo)存在量綱不統(tǒng)一的問題。為了解決這個問題，使各指標(biāo)間存在可比性，就必須對其做歸一化處理，使各指標(biāo)值位于0～1。歸一化公式如下：

式中：Rij為j區(qū)第i個指標(biāo)的標(biāo)準(zhǔn)化值；Xij為j區(qū)第i個指標(biāo)值；maxjXij和minjXij分別為j區(qū)第i個指標(biāo)中的最大值和最小值。

1.3.6指標(biāo)權(quán)重的確定 為了增加指標(biāo)權(quán)重的客觀性，選擇熵權(quán)法確定各指標(biāo)權(quán)重。該方法的基本原理是若某個指標(biāo)的信息熵越小，則說明指標(biāo)值的變異程度越大，提供的信息量越大，在綜合評價中所占比重越大，其權(quán)重值也就越大。反之，亦然。具體計算方法見文獻(xiàn)［36］。

1.3.7草原干旱風(fēng)險評價模型的建立 構(gòu)成風(fēng)險的4 項要素中，危險性、暴露性、敏感性、脆弱性與風(fēng)險生成的作用方向是一致的，因此，將草原干旱風(fēng)險評價模型定義為：

式中：F為草原干旱風(fēng)險；H、E、S和V分別為危險性、暴露性、敏感性和脆弱性；WH、WE、WS和WV分別為危險性、暴露性、敏感性和脆弱性對應(yīng)的權(quán)重。

2 結(jié)果與分析

2.1 致災(zāi)因子危險性

由公式（1）可得中國草原主要生長季5-9 月的干旱危險性。中國草原干旱的高及次高危險區(qū)主要位于內(nèi)蒙古高原、大興安嶺、寧夏以及甘肅的河西走廊地區(qū)和河?xùn)|地區(qū)（圖2），柴達(dá)木和塔里木盆地邊緣，青藏高原西部以及云貴高原。低危險區(qū)主要位于青藏高原東部和祁連山區(qū)。通過區(qū)域統(tǒng)計發(fā)現(xiàn)，高及次高危險區(qū)占草原總面積的39.4%。

圖2 草原干旱的危險性空間分布Fig.2 The distribution of the drought hazard for grassland

2.2 承災(zāi)體暴露性

中國草地NPP 具有從東南向西北逐漸減少的趨勢。對于南方的草山草坡來說，草地NPP 在500 g C·m-2以上，華南地區(qū)水熱條件好，草地NPP 可達(dá)800 g C·m-2以 上。對 于 青 藏 高 原，草 地NPP 在200 g C·m-2左右。新疆北部、內(nèi)蒙古東部和甘肅南部的草地NPP 在350 g C·m-2以上。荒漠草原帶的草地NPP 在100 g C·m-2以下。以單位面積草地NPP 代表承災(zāi)體暴露性（圖3），可以看出高及次高暴露性區(qū)域主要位于中國的華南和華東地區(qū)，低及次低暴露區(qū)主要位于內(nèi)蒙古、山西、中國西北地區(qū)以及青藏高原。通過區(qū)域統(tǒng)計發(fā)現(xiàn)，高及次高暴露區(qū)占草地總面積的7.2%。

圖3 承災(zāi)體暴露性空間分布Fig.3 The distribution of the drought exposure for grassland

2.3 承災(zāi)體敏感性

從草地EVI 變異系數(shù)來看，高變異區(qū)主要位于內(nèi)蒙古高原、甘肅河西走廊、寧夏地區(qū)、準(zhǔn)格爾和塔里木盆地邊緣、青藏高原南部、川西地區(qū)和云南北部（圖4）。低變異區(qū)位于內(nèi)蒙古西部、青藏高原腹地及中國南方地區(qū)。通過分析不同草地類型的生境，發(fā)現(xiàn)青藏高原東部和南部、天山山脈以及祁連山中段的主要草地類型為高寒草甸，生境寒冷而潮濕，植被具有耐寒的特性；內(nèi)蒙古北部及東北地區(qū)、塔里木和柴達(dá)木盆地邊緣的主要草地類型是山地和低地草甸類，生境水濕條件較好；云南、四川盆地及長江以南地區(qū)的主要草地類型是熱性草叢類和灌草叢類，山西、河北、河南的主要草地類型是暖性草叢和灌草叢類，生境濕潤度較大，主要是森林、灌叢退化后形成的。以上草地類型對干旱的響應(yīng)敏感，屬于高敏感區(qū)。位于干旱區(qū)的內(nèi)蒙古中部和西部、新疆北部、青藏高原北邊緣、西藏北部、寧夏、甘肅河西走廊等地的主要草地類型是溫性和高寒荒漠、溫性和高寒荒漠化草原、溫性草原化荒漠及溫性草原，生境干旱缺水，建群種對干旱的適應(yīng)能力強，屬于較低敏感區(qū)。

圖4 承災(zāi)體敏感性指標(biāo)及空間分布Fig.4 Sensitivity indices for grass and their spatial distribution

由公式（12）得，草原干旱敏感性較高的區(qū)域主要位于內(nèi)蒙古東南部、東北地區(qū)、青藏高原南部和東部邊坡、塔里木盆地北邊緣以及西南地區(qū)。敏感性較低的區(qū)域主要位于內(nèi)蒙古西部、新疆北部和青藏高原西北部。通過區(qū)域統(tǒng)計發(fā)現(xiàn)，高及次高敏感區(qū)占草原總面積的24.0%。

2.4 孕災(zāi)環(huán)境脆弱性

從土壤有效水含量指數(shù)來看，脆弱性較高的區(qū)域主要位于內(nèi)蒙古、新疆北部、青藏高原、云貴高原和四川東部地區(qū)（圖5）。這是由于土壤有效水含量越高，土壤中潛在的可被植物吸收利用的水分越多。從坡度指數(shù)可以看出，脆弱性較高的區(qū)域主要位于青藏高原邊緣區(qū)域及青藏高原南部、川西高原、天山山脈和祁連山山脈。從坡向指數(shù)來看，高脆弱區(qū)主要位于陽坡和半陽坡。這是因為陽坡、半陽坡的日照時間長、太陽輻射強，水分條件相對較差，陰坡和半陰坡正好相反，水分條件相對較好。而草地主要分布在降水較少的半濕潤、半干旱和干旱地區(qū)，使得水分成為草原干旱孕災(zāi)環(huán)境脆弱性的主導(dǎo)因素。

根據(jù)公式（13），獲得中國草原干旱孕災(zāi)環(huán)境脆弱性，孕災(zāi)環(huán)境脆弱性較高的區(qū)域主要分布在內(nèi)蒙古中部和東部、甘肅河西走廊、新疆北部、青藏高原大部分地區(qū)、云南和川東地區(qū)（圖5）。脆弱性較低的區(qū)域主要位于內(nèi)蒙古西北部、青藏高原北部和東部邊坡。通過區(qū)域統(tǒng)計發(fā)現(xiàn)，高及次高脆弱區(qū)占草原總面積的43.1%。

圖5 孕災(zāi)環(huán)境脆弱性指標(biāo)及空間分布Fig.5 Vulnerability indices for grass and their spatial distribution

2.5 草原干旱風(fēng)險

用熵權(quán)法獲得四大因子在干旱風(fēng)險形成中的權(quán)重值，危險性、暴露性、敏感性和脆弱性的權(quán)重分別為0.30、0.21、0.24 和0.25，通過公式（15）在ArcGIS 軟件中得到中國草原干旱風(fēng)險空間分布圖，再通過自然斷點法［37］將干旱風(fēng)險F劃分為高風(fēng)險區(qū)、次高風(fēng)險區(qū)、中風(fēng)險區(qū)、次低風(fēng)險區(qū)和低風(fēng)險區(qū)，對應(yīng)的F分別為＜0.23，0.23～0.40，0.40～0.54，0.54～0.69，＞0.69。中國草原干旱高及次高風(fēng)險區(qū)主要位于內(nèi)蒙古東北部的呼倫貝爾草原和錫林郭勒草原、東北西部科爾沁草原、寧夏北部、塔里木盆地北邊緣、青藏高原南部、云貴高原、河南和山東地區(qū)（圖6）。較低風(fēng)險區(qū)主要位于內(nèi)蒙古西部、甘肅河西走廊、新疆北部和青藏高原東北部。通過區(qū)域統(tǒng)計發(fā)現(xiàn)，高及次高風(fēng)險區(qū)占草原總面積的13.9%。

圖6 中國草原干旱風(fēng)險空間分布Fig.6 Spatial distribution of grassland drought risk in China

3 討論

從研究結(jié)論可以看出，各地草原干旱高風(fēng)險性的主導(dǎo)因素是不同的。內(nèi)蒙古東北部和新疆塔里木盆地北邊緣主要是由高的致災(zāi)危險性、承災(zāi)體敏感性和脆弱性決定的；青藏高原西南部主要是由高的敏感性和脆弱性決定的；云貴高原主要是由高的危險性、暴露性、敏感性和脆弱性決定的；南方地區(qū)主要是由高的暴露性決定的。雖然內(nèi)蒙古西部的干旱危險性很高，但由于該區(qū)域的暴露性、敏感性和脆弱性較低，因此屬于干旱的低風(fēng)險區(qū)。因此，各地區(qū)應(yīng)因地制宜地制定防旱抗旱措施，以便最大限度地降低草原干旱風(fēng)險。例如降低致災(zāi)危險性，實施人工影響天氣；通過科學(xué)的草原管理技術(shù)，提高草地的抗旱性，并酌情發(fā)展人工高產(chǎn)飼草料基地，降低承災(zāi)體敏感性；加大草原畜牧業(yè)的基礎(chǔ)設(shè)施建設(shè)，例如棚圈；增強氣象干旱及其風(fēng)險的監(jiān)測、預(yù)報和預(yù)警技術(shù)，提高服務(wù)的時效性；通過發(fā)展規(guī)?；投嘣慕?jīng)營，改善牧民的收入結(jié)構(gòu)，提高牧民對草原干旱的防旱抗旱能力；建立草原生態(tài)紅線和生態(tài)補充機制，提高草原生態(tài)系統(tǒng)的干旱防御能力；在個別地區(qū)還可以實施生態(tài)移民，進行草原生態(tài)環(huán)境保護。

本研究從不同方面展示了中國草原干旱的危險性、暴露性、敏感性和脆弱性特征，并繪制了草原干旱風(fēng)險區(qū)劃圖。該結(jié)果可以為中國草原干旱風(fēng)險的管理提供一定的理論依據(jù)，但依然存在諸多不足。草原干旱風(fēng)險是一個復(fù)雜的系統(tǒng)，該系統(tǒng)由若干個子系統(tǒng)組成，而各子系統(tǒng)之間又存在著隨機性和動態(tài)性的關(guān)聯(lián)，如何構(gòu)建科學(xué)完整的風(fēng)險評價指標(biāo)體系，并定量化各指標(biāo)之間的關(guān)系就顯得非常重要。在構(gòu)建草原干旱風(fēng)險評價指標(biāo)體系前應(yīng)進行充分的干旱調(diào)研和風(fēng)險分析，不僅要考慮數(shù)據(jù)的代表性和連續(xù)性，還要考慮數(shù)據(jù)的可獲取性和時空分辨率。未來應(yīng)結(jié)合大量實測數(shù)據(jù)驗證評價結(jié)果，不斷優(yōu)化評價指標(biāo)和評價模型。并且，草地生長是一個動態(tài)的復(fù)雜過程，干旱對不同生長階段的影響不同，在今后的研究中也應(yīng)關(guān)注不同時期的草原干旱風(fēng)險特征。

4 結(jié)論

基于中國草原干旱風(fēng)險形成機理，從危險性、暴露性、敏感性和脆弱性四方面入手，選擇具有代表性的7 個指標(biāo)，建立了草原干旱風(fēng)險評估模型，將中國草原干旱風(fēng)險分為5 級并進行評估與區(qū)劃，得到以下結(jié)論：1）高及次高危險區(qū)主要位于內(nèi)蒙古地區(qū)、河西走廊、柴達(dá)木和塔里木盆地邊緣、青藏高原西部以及云南地區(qū)。2）高及次高暴露區(qū)主要位于中國的西南和華南地區(qū)。3）敏感性較高的區(qū)域主要位于內(nèi)蒙古東南部、東北地區(qū)、青藏高原南部和東部邊坡、塔里木盆地北邊緣以及西南地區(qū)。4）脆弱性較高的區(qū)域主要分布在內(nèi)蒙古中部和東部、甘肅河西走廊、新疆北部、青藏高原大部分地區(qū)、云南和川東地區(qū)。5）干旱高及次高風(fēng)險區(qū)主要位于內(nèi)蒙古東北部、東北西部、寧夏北部、塔里木盆地北邊緣、青藏高原南部、云貴高原、河南和山東地區(qū)。從草地生態(tài)系統(tǒng)服務(wù)總價值來看，極端干旱影響下其變化趨勢為草甸草原＞典型草原＞荒漠草原［38］；從草地水分利用效率來看，青藏高原西北部和內(nèi)蒙古北部最低［39］。這也從另一方面驗證了評估結(jié)果的可靠性。