劉振男 周靖楠 龔宇 徐敏

摘要：

構(gòu)建有效干旱指數(shù)是識別、評價實際旱情的關(guān)鍵，以貴州省為研究區(qū)域，基于月尺度降水、氣溫、土壤含水量數(shù)據(jù)，以GB/T 20481-2017《氣象干旱等級》中綜合干旱指數(shù)（CI）的構(gòu)建方法為參考，采用主成分分析對上述多源數(shù)據(jù)集進行融合計算，提出了一種物理機制清晰、應(yīng)用簡便的綜合干旱指數(shù)（PTS）。評價了該指數(shù)在貴州省的適用性，檢驗了該指數(shù)表征旱情時空演變?nèi)^程特征的能力，分析了該指數(shù)與其他常用干旱指數(shù)間的性能差異。結(jié)果表明：在干旱歷時與強度方面，新構(gòu)建的綜合干旱指數(shù)能夠準確地識別出貴州省全部6次典型旱情，同時其刻畫的干旱發(fā)生-發(fā)展-消退全過程時空演變特征與歷史資料基本一致。該干旱指數(shù)不但克服了典型氣象干旱指數(shù)出現(xiàn)的高估旱情等級問題，而且對旱情態(tài)勢走向評價更加準確。

關(guān) 鍵 詞：

干旱指數(shù)； 旱情評價； 主成分分析； 貴州省

中圖法分類號： TV11

文獻標志碼： A

DOI：10.16232/j.cnki.1001-4179.2023.08.008

0 引 言

近年來，隨著全球氣候變暖和人類活動加劇，中國傳統(tǒng)濕潤區(qū)貴州省干旱事件頻發(fā)，給當?shù)亟?jīng)濟造成了巨大損失 ［1-2］。準確刻畫貴州省干旱全過程特征，不僅有助于指導抗旱減災(zāi)工作的具體部署，對深入剖析喀斯特地區(qū)干旱形成機理也具有重要研究價值。目前，對貴州省干旱評價的研究主要基于單要素干旱指數(shù)，如毛春艷等［3］采用降水距平百分率指數(shù)（PI）對貴州山區(qū)的干旱情勢進行了評價并指出當?shù)馗珊登閯菖c地形因子關(guān)系密切；吳建峰等［4］基于標準降水指數(shù)（SPI）分析了貴州近半世紀的氣象干旱時空演變特征；劉振男等［5］基于中國氣象局全球大氣/陸面再分析產(chǎn)品（CMA-RA）計算了貴州標準化土壤濕度指數(shù)（SSI）并驗證了該指數(shù)對研究區(qū)農(nóng)業(yè)干旱評價的適用性。由于各研究切入點與數(shù)據(jù)來源不同，僅從氣象、農(nóng)業(yè)等視角片面地對干旱情勢進行評價，無法全面客觀地反映干旱特征全貌。為此，能夠多維度識別干旱特征的指數(shù)也被引入到研究區(qū)的應(yīng)用當中。例如，慎東方等［6］選用標準化降水蒸發(fā)指數(shù)（SPEI）對貴州氣象干旱的歷時與烈度進行了識別與評價，結(jié)果表明同時考慮降水與蒸發(fā)的指數(shù)所得結(jié)果更加可靠；吳哲紅等［7］分別運用綜合氣象干旱指數(shù)（CI）對貴州安順地區(qū)的氣象干旱進行了評價，結(jié)果表明該指數(shù)的評定結(jié)果更能反映實際旱情特征；王文等［8］基于遙感監(jiān)測數(shù)據(jù)研究了標準化降水指數(shù)（SPI）、偵測干旱指數(shù)（RDI）及蒸散發(fā)脅迫指數(shù)（ESI）在云貴高原區(qū)的干旱評價效果，結(jié)果證實同時考慮實際蒸散發(fā)和潛在蒸散發(fā)的蒸

散發(fā)脅迫指數(shù)比考慮單一蒸散發(fā)因素的干旱指數(shù)在監(jiān)測干旱方面更有效；張建平等［9］基于相對濕潤指數(shù)（MI）與歸一化植被指數(shù)（NDVI），構(gòu)建了綜合干旱指數(shù)（DI）并驗證了其對西南地區(qū)的干旱評價效果優(yōu)于其他單一干旱指數(shù)。由此可見，采用多要素干旱指數(shù)對研究區(qū)的干旱評價更加可靠，但當前用于研究區(qū)的綜合干旱指數(shù)數(shù)據(jù)來源過于單一，或均采集自監(jiān)測站點，或均來源于遙感監(jiān)測，而采用多源數(shù)據(jù)耦合構(gòu)造綜合干旱指數(shù)的研究并不多見，對其評價效果的適用性也沒有充分認識。因此，本文基于監(jiān)測站點的降水、氣溫數(shù)據(jù)及全球大氣/陸面再分析產(chǎn)品的土壤含水量數(shù)據(jù)，采用主成分分析方法構(gòu)建一種適用于評價貴州省干旱的綜合干旱指數(shù)，并與SPEI指數(shù)和SSI指數(shù)進行對比，驗證了新構(gòu)建指數(shù)在貴州省干旱評價中的適用性。研究成果可為提升研究區(qū)干旱監(jiān)測預(yù)警能力提供技術(shù)參考。

1 研究區(qū)域與方法

1.1 研究區(qū)概況

貴州省（東經(jīng)103°36′～109°31′，北緯24°37′～29°13′）位于中國西南，是典型的喀斯特地貌集中區(qū)，山地與丘陵占貴州省總面積的90%以上，地勢呈西高東低之態(tài)，平均海拔為1 100 m，為濕潤地區(qū)，屬亞熱帶濕潤季風氣候，年內(nèi)氣溫差異性較小，多年平均氣溫約為15 ℃，雨量大但分布不均，降水主要集中在5～10月，多年平均降水量約為1160.6 mm，常年相對濕度在70%以上，局部氣候特征顯著，干旱、凌凍等災(zāi)害性天氣多發(fā)，給當?shù)厣鐣斐闪司薮笪：?，引起了當?shù)卣拖嚓P(guān)組織的高度關(guān)注［10］。

1.2 數(shù)據(jù)來源

貴州省37個氣象站的降水與氣溫數(shù)據(jù)收集自中國國家氣象信息中心（NMIC），以數(shù)據(jù)完整度與系列長度為標準，對各站點的數(shù)據(jù)進行了再次遴選，最終確定選用19個氣象站點的數(shù)據(jù)進行研究，上述站點空間分布均勻，具有較高的代表性與可靠性，具體氣象站點信息如圖1所示。土壤含水量數(shù)據(jù)采集自中國氣象局全球大氣/陸面再分析產(chǎn)品（CMA-RA），由于0～10 cm土層的含水量能夠有效表征農(nóng)業(yè)干旱的實際情況［11-12］，為此，本文選用的0～10 cm土壤含水量數(shù)據(jù)提取自同化產(chǎn)品“CRA40LAND”，空間分辨率為0.5°×0.5°，時間分辨率為1個月。降水、氣溫與土壤含水量數(shù)據(jù)的時間起止為2009年1月至2020年12月，上述數(shù)據(jù)資料下載自中國氣象數(shù)據(jù)服務(wù)網(wǎng)（http：∥data.cma.cn/）。

1.4 綜合干旱指數(shù)評定等級的確定

綜合干旱指數(shù)（PTS）應(yīng)用的關(guān)鍵在于確定其合理的干旱評定等級。本文參考文獻［18］的確定方法，即基于月時間尺度的SPEI指數(shù)和SSI指數(shù)對應(yīng)的不同干旱發(fā)生頻率均值（見表1），反算得到了PTS的干旱評定等級標準（見表2）。SSI指數(shù)的計算原理與SPEI指數(shù)相近，具體可參考文獻［5］。

2 結(jié)果與分析

2.1 綜合干旱指數(shù)適用性分析

為了檢驗綜合干旱指數(shù)在貴州省干旱事件評價的適用性，依據(jù)水利部發(fā)布的《中國水旱災(zāi)害公報》（http：∥www.mwr.gov.cn/sj/tjgb/zgshzhgb/），整理了貴州省2009～2020年間的實際旱情資料。據(jù)統(tǒng)計，貴州省12 a間有6 a發(fā)生了干旱，具體時間分別為：2009年、2010年、2011年、2012年、2014年與2018年。根據(jù)PTS的干旱等級評定標準，將PTS≤-0.36作為干旱發(fā)生的臨界點，將其持續(xù)時間作為干旱歷時，干旱歷時期中PTS的最小值作為干旱強度的判定條件。

表3和圖2展示的是貴州省旱情歷史資料與綜合干旱指數(shù)PTS間的比對結(jié)果?？梢钥闯?，新構(gòu)建的PTS指數(shù)較好地識別出發(fā)生在貴州省的4次干旱事件，PTS指數(shù)所顯示的干旱歷時與干旱強度與實際資料基本相符，尤其是對2009年、2010年及2011年干旱事件的評價更加準確。例如，就干旱歷時比對情況而言，歷史資料指出“自2009年9月中旬，貴州出現(xiàn)旱情，到2010年5月，旱情減緩”，而PTS指數(shù)顯示此次旱情起止時間約為2009年9月至2010年5月；就干旱強度比對情況而言，歷史資料指出：2011年9月，旱情達到高峰期，而PTS指數(shù)的最小值亦出現(xiàn)在2011年9月份，數(shù)值為-1.25，旱情等級為重旱。由此可見，PTS指數(shù)能對貴州省實際干旱情況進行準確的評價，說明其適用于研究區(qū)的干旱評價，且評價結(jié)果具有較好的可靠性。

2.2 綜合干旱指數(shù)監(jiān)測典型干旱能力分析

選用貴州省2009～2010年秋冬春三季連旱作為典型案例，檢驗綜合干旱指數(shù)PTS監(jiān)測旱情過程的能力。2009年與2010年的《中國水旱災(zāi)害公報》記錄：此次旱情開始于2009年9月中旬，貴州東北部、北部、西北部與東南部地區(qū)最早出現(xiàn)旱情，2009年12月旱情開始迅速發(fā)展，直至2010年4月末，旱情逐漸開始緩解，2010年3，4月旱情最為嚴重，部分地區(qū)干旱強度達到特旱。

基于Arcgis軟件平臺，采用反距離插值繪制出綜合干旱指數(shù)PTS監(jiān)測2009～2010年旱情時空演變過程（見圖3）。就時間尺度而言，此次旱情最早發(fā)生于2009年9月，見圖3（a），且研究區(qū)局部有5個月出現(xiàn)了特旱，即2009年9月、2009年11月、2010年1月、2010年2月及2010年3月，PTS的峰值出現(xiàn)于2010年3月（PTS≈-2.44，特旱），自2010年4月開始，旱情逐漸開始減弱。就空間尺度分析，旱情最早發(fā)生于貴州西部、北部、東北部及東南部，旱情等級均達重旱以上，見圖3（a），2009年10月旱情出現(xiàn)了大幅蔓延之勢，見圖3（b），直至2010年1月旱情開始穩(wěn)定擴散至全省，見圖3（e），2010年2月與3月旱情最為嚴重，貴州東南部、東北部及西南部地區(qū)旱情等級已到達特旱，見圖3（f）與圖3（g），旱情于2010年4月開始有所減弱，見圖3（h）。顯然，綜合干旱指數(shù)PTS反映的旱情時空演變特征與歷史資料描述情況基本一致，說明PTS具備良好的旱情監(jiān)測能力，可用于評價貴州省實際旱情的指數(shù)。

2.3 綜合干旱指數(shù)與其他干旱指數(shù)對比分析

基于月尺度的氣象干旱指數(shù)（SPEI）、農(nóng)業(yè)干旱指數(shù)（SSI）與綜合干旱指數(shù)（PTS）計算了不同干旱指數(shù)間的相關(guān)系數(shù)，結(jié)果表明PTS與SPEI、SSI間的相關(guān)系數(shù)值分別為0.682 8及0.879 2，顯然PTS與SPEI、SSI均呈正相關(guān)關(guān)系，這與PTS指數(shù)的構(gòu)建原理相符，表明了PTS在評價干旱時能同時兼顧降水、氣溫與土壤含水量的影響；同時發(fā)現(xiàn)，由于計算SPEI與SSI時的基礎(chǔ)數(shù)據(jù)無交叉，因此SPEI與SSI間的相關(guān)系數(shù)僅為0.126 1。PTS、SPEI與SSI指數(shù)監(jiān)測研究區(qū)干旱的變化過程見圖4。由圖4可知，氣象干旱指數(shù)SPEI隨著時間變化曲線波動明顯，同時對干旱過程反應(yīng)敏感，能夠較早捕捉到干旱未來發(fā)展態(tài)勢信息，但識別干旱穩(wěn)定性一般，會出現(xiàn)高估旱情等級的現(xiàn)象。農(nóng)業(yè)干旱指數(shù)SSI對實際旱情評價較為客觀，與歷史資料信息反映的情況大體相符，但有評價遲滯現(xiàn)象，即SSI反映旱情時間往往落后于實際情況。綜合干旱指數(shù)PTS的評價結(jié)果與SSI的評價結(jié)果較為接近，但PTS克服了SSI的遲滯現(xiàn)象，對干旱的反映更加及時。

為了深入剖析不同干旱指數(shù)間的性能差異，以2009～2010年貴州省典型干旱事件為例，再次對3種干旱指數(shù)的評價過程進行分析?？傮w而言，3種干旱指數(shù)均對典型旱情做出了有效識別，即干旱歷時均為8個月，干旱強度均為特旱。就評價旱情過程細節(jié)而言，SPEI顯示2009年9月干旱等級就達到了特旱，并持續(xù)6個月，直到2010年3月逐漸緩解；SSI顯示2009年9月起旱情為重旱水平，2010年2月與3月旱情達特旱；PTS同樣顯示2010年2月與3月旱情為特旱水平，但在2009年9～12月，旱情有從輕旱向重旱過渡階段，這與歷史資料所描述的情況更加相符。就旱情等級評價細節(jié)而言，歷史資料指出“本輪旱情在2010年3月與4月份最為嚴重”，PTS顯示3月份旱情為特旱、4月份旱情為重旱，SPEI顯示3月份旱情為輕旱、4月份旱情為無旱，SSI顯示3月份旱情為特旱、4月份旱情為輕旱。顯然，PTS的評價結(jié)果較SPEI，SSI更加貼合實際情況，即對4月份的旱情走勢評價更為準確。

3 討 論

本文基于降水、氣溫以及土壤含水量數(shù)據(jù)，采用主成分分析法構(gòu)建了一種適用于貴州省的綜合干旱指數(shù)PTS，以《中國水旱災(zāi)害公報》的記載資料為背景，檢驗了PTS評價實際旱情的能力，且探明了其與SPEI、SSI指數(shù)在識別旱情方面的差異性，研究成果可為貴州省干旱防治工作提供一定的技術(shù)參考。

雖然研究取得些許成果，但仍有不足之處：① 致旱要素眾多，諸如降水、氣溫、徑流、日照時數(shù)、地貌特征等都會對干旱產(chǎn)生一定影響［19］，而本文僅考慮了降水、氣溫與土壤含水量的耦合作用，后續(xù)可在此研究基礎(chǔ)上，適時增加其他相關(guān)影響要素，進而提高綜合干旱指數(shù)的普適性與可靠性，為準確刻畫研究區(qū)干旱全貌特征提供有力工具；② 當前干旱指數(shù)種類眾多，本文僅對PTS與SPEI、SSI指數(shù)的性能差異性進行了對比分析，未將PTS與其他不同類型干旱指數(shù)進行對比分析，今后可作進一步的討論、驗證；③ 旱情嚴重期，對干旱實時評價的時間往往以日、候為尺度［20］，本文受限于數(shù)據(jù)資料時間分辨率的客觀原因，未以日、候的時間尺度對PTS進行計算分析，下一步可采用降尺度等方法提高數(shù)據(jù)的時間分辨率，從而滿足政府與相關(guān)部門的實際業(yè)務(wù)需求。

4 結(jié) 論

本文收集了貴州省月尺度降水、氣溫與土壤含水量數(shù)據(jù)，基于主成分分析法構(gòu)建了綜合干旱指數(shù)PTS，檢驗了PTS指數(shù)在貴州省干旱評價中的適用性與可靠性，得到結(jié)論如下：

（1） 綜合干旱指數(shù)PTS能夠很好地識別貴州省的實際旱情，是一種研究貴州省干旱的有效指數(shù)。

（2） 綜合干旱指數(shù)PTS具備良好的旱情監(jiān)測能力，其反映的旱情時空演變特征與歷史資料描述情況基本一致。

（3） 綜合干旱指數(shù)PTS、SPEI與SSI指數(shù)在評價干旱情勢方面能力大體相當，但在旱情細節(jié)評價方面，PTS較SPEI、SSI指數(shù)更有優(yōu)勢。

參考文獻：

［1］ 劉振男，杜堯，韓幸燁，等.基于遺傳算法優(yōu)化極限學習機模型的干旱預(yù)測：以云貴高原為例［J］.人民長江，2020，51（8）：13-18.

［2］ LI Q，DU Y，LIU Z，et al.Drought prediction in the Yunnan-Guizhou Plateau of China by coupling the estimation of distribution algorithm and the extreme learning machine［J］.Nat Hazards，2022，113：1635-1661.

［3］ 毛春艷，戴麗，楊廣斌，等.1960-2016年貴州喀斯特山區(qū)干旱時空動態(tài)分析［J］.水資源與水工程學報，2021，32（3）：64-72，79.

［4］ 吳建峰，張鳳太，盧海芬，等.基于標準化降水指數(shù)的貴州省近54年干旱時空特征分析［J］.科學技術(shù)與工程，2018，18（15）：207-214.

［5］ 劉振男，焦衛(wèi)國，王笑宇，等.基于標準化土壤濕度指數(shù)的貴州省農(nóng)業(yè)干旱評價［J］.中國農(nóng)村水利水電，2022（11）：1-5，12.

［6］ 慎東方，商崇菊，方小宇，等.貴州省干旱歷時和干旱烈度的時空特征分析［J］.干旱區(qū)資源與環(huán)境，2016，30（7）：138-143.

［7］ 吳哲紅，詹沛剛，陳貞宏，等.基于CI指數(shù)和K指數(shù)的近40a貴州安順區(qū)域干旱演變特征［J］.冰川凍土，2013，35（4）：1044-1055.

［8］ 王文，黃瑾，崔巍.云貴高原區(qū)干旱遙感監(jiān)測中各干旱指數(shù)的應(yīng)用對比［J］.農(nóng)業(yè)工程學報，2018，34（19）：131-139，309.

［9］ 張建平，劉宗元，王靖，等.西南地區(qū)綜合干旱監(jiān)測模型構(gòu)建與驗證［J］.農(nóng)業(yè)工程學報，2017，33（5）：102-107.

［10］ 劉振男.中長期干旱預(yù)測研究：以云貴高原為例［D］.南京：河海大學，2019.

［11］ 牛乾坤，劉瀏，程湫雅，等.基于多源遙感數(shù)據(jù)的河套灌區(qū)干旱時空演變特征［J］.干旱地區(qū)農(nóng)業(yè)研究，2020，38（4）：266-277.

［12］ SHEFFIELD J，GOTETI G，WEN F，et al.A simulated soil moisture based drought analysis for the United States［J］.Journal of Geophysical Research，2004，109：1-19.

［13］ 吳志勇，程丹丹，何海，等.綜合干旱指數(shù)研究進展［J］.水資源保護，2021，37（1）：36-45.

［14］ 袁星，王鈺淼，張苗，等.關(guān)于驟旱研究的一些思考［J］.大氣科學學報，2020，43（6）：1086-1095.

［15］ 粟曉玲，梁曉萱，吳海江，等.表征生態(tài)狀況的綜合干旱指數(shù)構(gòu)建及干旱風險分析［J］.水資源保護，2023，39（2）：50-58，100.

［16］ 郭盛明.綜合干旱指數(shù)構(gòu)建及其適用性研究［D］.楊凌：西北農(nóng)林科技大學，2021.

［17］ 常文娟，梁忠民，馬海波.基于主成分分析的干旱綜合指標構(gòu)建及其應(yīng)用［J］.水文，2017，37（1）：33-38，82.

［18］ 馬俊超，杭慶豐，李瓊芳，等.基于VIC模型的灤河流域綜合干旱指數(shù)的構(gòu)建與應(yīng)用［J］.水資源與水工程學報，2015，26（2）：79-84.

［19］ 馮寶飛，邱輝，紀國良.2022年夏季長江流域氣象干旱特征及成因初探［J］.人民長江，2022，53（12）：6-15.

［20］ 虞美秀.綜合干旱指數(shù)的構(gòu)建與應(yīng)用［D］.南京：河海大學，2013.

（編輯：江 文）

Abstract：

The construction of an effective drought index is the key to identifying and evaluating the actual drought situation.Taking Guizhou Province as the study area，a comprehensive drought index with a clear physical mechanism and easy application was developed based on monthly precipitation，temperature and soil moisture content according to the method in GB/T 20481-2017.The above multi-source data sets were fused and calculated by Principal Component Analysis.The applicability of the index in the study area was evaluated.The ability of the index to represent the whole process of spatio-temporal evolution of drought was tested.The performance difference between this index and other commonly used drought indexes was analyzed.The results show that the newly constructed comprehensive drought index can accurately identify all 6 typical droughts in the study area in terms of drought duration and intensity.Meanwhile，the spatio-temporal characteristics of the whole process of drought occurrence，development and subside are basically consistent with the historical data.This drought index not only overcomes the problem of overestimation of drought grade that exists in typical meteorological drought index but also evaluates the drought trend more accurately.The research results can provide a technical reference for drought relief work in Guizhou Province.

Key words：

drought index；drought assessment；principal component analysis；Guizhou Province