西南地區氣象干旱向水文干旱傳播的特征

石 朋，唐 漢，瞿思敏，聞 童，趙蘭蘭，李瓊芳

(1.河海大學水文水資源學院，江蘇 南京 210098； 2.中國灌溉排水發展中心，北京 100054；3.水利部水文局，北京 100053)

干旱是世界上影響范圍最廣、發生頻率最高、持續時間最長、致災率最高的一類自然災害，也是眾多自然災害中造成經濟損失最大的一類氣象災害[1-2]。據統計，在各類自然災害造成的總損失中，氣象災害造成的損失約占85%，而干旱造成的損失又占氣象災害的50%左右[3]。自中華人民共和國成立以來，我國因干旱造成的直接經濟損失占干旱年GDP的2.5%～3.5%，干旱已成為制約我國國民經濟發展的重要因素。

在全球氣候變化的大環境下，我國干旱時空分布規律也隨之發生了變化，干旱災害逐漸從北方干旱區向南方氣候濕潤區發展[4-5]，其中，西南地區干旱發生頻率和強度明顯增加，引起了學者們的普遍關注[6-7]。如2009—2010年云南、貴州等地秋冬春連旱，造成了巨大的經濟損失[8-9]。因此，全面認識西南地區干旱特征，對提高干旱災害防御能力具有十分重要的意義。

根據研究目標的不同，干旱可分為氣象干旱、水文干旱、農業干旱和社會經濟干旱等不同類型[10-11]。通常情況下，氣象干旱和水文干旱之間可以認為具有線性關系，可通過對干旱響應時間等關鍵指標的分析來研究氣象干旱向水文干旱傳播的特征[12]。以往的研究多以識別得到的氣象干旱和水文干旱事件之間的滯時來確定干旱響應時間[13-15]。例如， Ding等[13]研究了中國不同氣候區干旱事件之間的傳播特征，發現大部分地區的干旱滯時不到2月?；跍r確定干旱響應時間的評價方法可以有效地確定較長持續時間的干旱事件之間的響應時間，但對于持續時間較短的干旱過程難以確定其氣象干旱和水文干旱之間是否存在傳播關系，而通過干旱指數之間的相關關系確定干旱響應時間，能夠考慮持續時間較短的干旱過程。

以往對西南地區干旱的研究多以氣象干旱為主[16-18]，而綜合氣象干旱和水文干旱探究西南干旱的研究較少。本文利用標準化降雨蒸散發指數(standardized precipitation evapotranspiration index，SPEI)和標椎化徑流指數(standardized runoff index，SRI)序列確定干旱響應時間，并以干旱響應時間為約束條件，基于游程理論和數學函數建立了干旱傳播模型，并量化了西南地區氣象干旱特征和水文干旱特征之間的關系，以期為利用氣象數據預測西南地區水文干旱提供參考。

1 研究區概況和研究數據

1.1 研究區概況

研究區(圖1)位于我國西南的長江上游地區，包括四川、貴州、云南、重慶和廣西5個省(自治區，直轄市)，地處四川盆地和云貴高原，東鄰長江中下游平原，以濕潤半濕潤的亞熱帶季風氣候為主，多年平均降水量約為1 058 mm，標準差為58 mm，年平均氣溫約為16.25℃，受亞熱帶季風和西南季風影響，夏季炎熱多雨，冬季降水偏少，降水量時空分布極不均勻，整體呈東多西少的分布形態。研究區內共有8個水文站，控制流域分別為元江流域(2.15萬km2)、江邊街流域(3.16萬km2)、三匯流域(3.32萬km2)、允景洪流域(7.26萬km2)、武隆流域(7.61萬km2)、五通橋流域(12.09萬km2)、梧州流域(32.70萬km2)和萬縣流域(63.08萬km2)，8個流域均為山區流域，且徑流過程受積雪融化及水利工程的調節影響較小，近似天然徑流過程。研究區由降水不足的氣象干旱引起的水文干旱多發，基本可視為線性關系。

圖1 研究區及氣象、水文站點位置Fig.1 Location of study area, meteorological and hydrological stations

1.2 研究數據

1968—2017年101個氣象站逐日降水數據由中國氣象數據服務中心(data.cma.cn)提供，利用日照時間和平均氣溫兩個氣候變量計算潛在蒸散發量并按月匯總，在計算得到單站點的潛在蒸散發量后，選取各流域內或鄰近流域的氣象站點并用算術平均法計算流域面平均降水量和面平均潛在蒸散發量以計算氣象干旱指數。依據水利部信息中心提供的8個水文站1968—2017年逐月自然徑流數據計算水文干旱指數。

2 研究方法

2.1 干旱指數

SPEI和SRI能描述不同時間尺度的干旱特征，在水文氣象領域得到了廣泛應用[19-21]。本文分別用SPEI和SRI來描述氣象干旱和水文干旱。 SPEI和SRI的計算過程如下: ①根據需要選擇合適的時間尺度(如3月的時間尺度一般用來描述季節性干旱，12月的時間尺度一般用來描述年干旱); ②用概率分布函數分別擬合降水量與潛在蒸散發量的差值序列與徑流序列; ③進行等概率變換，將偏態序列進行正態標準化處理。詳細計算過程可參見文獻[22]。

基于1968—2017年的氣象資料和徑流數據，計算不同時間尺度(1～12月)的SPEI(SPEI-1～SPEI-12)和1月尺度的SRI(SRI-1)。參照GB/T 20481—2017《氣象干旱等級》規定的標準，SPEI和SRI的干旱等級劃分如表1所示。

表1 SPEI和SRI的干旱等級劃分Table 1 Classification of drought grade on SPEI and SRI

2.2 干旱響應時間

已有研究大多通過計算不同時間尺度的干旱指數之間的相關關系或人為設置不同滯時的干旱指數序列來計算干旱響應時間，均明確表明水文干旱對氣象干旱的響應存在滯后現象[12-15]。本文利用皮爾遜相關系數[23](Pearson correlation coefficient，PCC)分析不同時間尺度的氣象干旱(SPEI-1～SPEI-12)和1月尺度的水文干旱(SRI-1)之間的相關性，將最大皮爾遜相關系數對應的SPEI時間尺度作為干旱響應時間。具體計算步驟如下：①計算SPEI-1～SPEI-12；②計算SRI-1；③計算SPEI-1～SPEI-12與SRI-1之間的皮爾遜相關系數，最大皮爾遜相關系數對應的SPEI時間尺度就是該流域的干旱響應時間。

2.3 干旱敏感性

由于流域的調蓄作用，并非每一次氣象干旱都會觸發水文干旱，水文干旱的發生也并非都由氣象干旱引起，本文以Nm- h表示由氣象干旱引起的水文干旱次數，Nh表示水文干旱總次數。此外，一次水文干旱過程可能由一次氣象干旱或者多次氣象干旱引起，以Nh-m表示引起水文干旱的氣象干旱次數，Nm表示氣象干旱總次數。Sattar等[24]提出了干旱傳播效率(記錄期內Nh與Nm的比值)來衡量水文干旱對氣象干旱的敏感性，干旱傳播效率越高，說明水文干旱對氣象干旱越敏感。Guo等[25]修正了干旱傳播效率，用Nh-m與Nm的比值來表示干旱傳播率。本文基于干旱響應時間和游程理論對氣象干旱和水文干旱進行匹配后，采用干旱傳播率來表征氣象干旱向水文干旱的傳播率(本文稱為氣象干旱傳播率(Tr1))，用Nm-h與Nh的比值來表示水文干旱對氣象干旱的敏感度(本文稱為水文干旱敏感度Tr2))，即：

Tr1=Nh-m/Nm

(1)

Tr2=Nm-h/Nh

(2)

2.4 干旱傳播模型的構建

干旱傳播模型的構建包括3個步驟: ①基于游程理論識別干旱事件并進行融合與剔除，提取干旱特征[4,17,26]; ②根據確定的干旱響應時間匹配干旱事件; ③根據干旱特征(干旱歷時D和干旱烈度S)構建傳播模型。本文采用SPEI-1和SRI-1來識別干旱事件，并提取干旱特征。干旱歷時是指干旱事件持續時間，干旱烈度是干旱歷時內干旱指標值與閾值之差的累積和，取其絕對值。

本文不僅通過游程理論對干旱事件進行識別，同時采用了干旱事件的融合和剔除方法(圖2)以減少歷時短烈度小的干旱事件對統計特征產生的影響。具體步驟如下[27]：①當干旱指數(R)小于R1，初步識別為干旱事件；②對于間隔時間為1月的2個干旱事件，若間隔月份的干旱指數小于R0，則將這2個干旱事件融合為1個干旱事件，如圖2中干旱事件a和b，將其融合為干旱事件2，其融合后的干旱歷時為6月，干旱烈度為干旱事件a和干旱事件b的烈度之和；③由于西南地區多發干旱歷時較長的干旱事件，對于小干旱事件可以剔除，即對于干旱歷時僅為1月的干旱事件，若該月份的干旱指數值大于R2，該干旱事件應被剔除。如圖2中干旱事件1，其干旱指數值小于R2，未被剔除。圖2中，R0取0；R1為判斷干旱事件的閾值，取-0.5；R2為判斷中旱事件的閾值，取-1。

圖2 干旱事件識別、融合和剔除過程示意圖Fig.2 Schematic diagram of drought events identifying, pooling and excluding process

識別干旱事件后，根據確定的干旱響應時間對氣象和水文干旱事件進行匹配。對于單個水文干旱事件，所匹配的為該水文干旱事件發生月份前干旱響應時間至該水文干旱事件結束月份之間所發生的氣象干旱事件。根據匹配得到的干旱事件的干旱特征，建立線性干旱傳播模型。

3 結果與分析

3.1 干旱事件識別

計算干旱指數后，根據游程理論識別干旱事件。西南地區元江流域、江邊街流域、三匯流域、允景洪流域、武隆流域、五通橋流域、梧州流域和萬縣流域分別發生氣象干旱57次、68次、36次、60次、68次、62次、59次和66次，發生水文干旱31次、30次、18次、28次、43次、31次、43和38次。各流域融合和剔除后的干旱特征如表2所示。

表2 干旱事件融合、剔除后的干旱特征值Table 2 The drought characteristic value of drought events after pooling and excluding

圖3為各流域融合和剔除后的氣象干旱和水文干旱的干旱特征箱形圖，氣象干旱歷時和烈度均較小且較為集中，而水文干旱歷時和烈度整體較為分散且數值較大，水文干旱嚴重性強于氣象干旱。

(a) 干旱歷時

(b) 干旱烈度圖3 干旱事件剔除后的干旱特征箱形圖Fig.3 Box chart of drought characteristics of drought events after excluding

3.2 干旱響應時間

采用最大皮爾遜相關系數來定量確定各流域氣象干旱到水文干旱的響應時間。表3為各流域最大皮爾遜相關系數及對應的SPEI時間尺度，整體來看，西南地區干旱響應時間為2～7月。圖4為各流域皮爾遜相關系數箱形圖，三匯流域皮爾遜相關系數偏小，允景洪流域皮爾遜相關系數較為集中。以梧州流域分析，梧州流域最大皮爾遜相關系數對應的SPEI時間尺度為2月，SRI-1與SPEI-2相關性最好，梧州流域干旱響應時間為2月。

表3 最大皮爾遜相關系數及對應的SPEI時間尺度Table 3 The maximum PCC and the correspondent time scale of SPEI

圖4 各流域皮爾遜相關系數箱形圖Fig.4 Box chart of PCC in each basin

以降水量和潛在蒸散發量計算得到的SPEI和以徑流量計算得到的SRI之間的皮爾遜相關系數反映了流域內徑流對降雨的響應關系。三匯流域和萬縣流域中下游多為平原地形，流域下墊面調蓄作用使得其SPEI和SRI相關性較弱，而多為高原和山地地形的元江流域和梧州流域的SPEI和SRI相關性較強。此外，在干旱傳播過程中，各流域干旱響應時間是多種因素綜合影響的結果，對于氣象因子和下墊面因子存在差異的不同流域，水文干旱對氣象干旱的響應關系也存在差異，可能導致不同的干旱響應時間。

3.3 干旱事件的匹配

將干旱響應時間用于匹配每個流域的水文干旱和氣象干旱事件，元江流域、江邊街流域、三匯流域、允景洪流域、武隆流域、五通橋流域、梧州流域和萬縣流域分別匹配干旱過程23次、23次、11次、27次、38次、24次、38次和23次，其干旱特征統計見表4。

表4 干旱事件匹配后的干旱特征統計Table 4 Statistics of drought characteristics of drought events after matching

圖5為干旱事件匹配后的干旱歷時和干旱烈度的統計特征，各流域水文干旱歷時和烈度均大于氣象干旱，且氣象干旱歷時和烈度相對集中，說明各流域氣象干旱的干旱特征變化幅度均較小，水文干旱比氣象干旱嚴重且干旱特征變化幅度較大。

(a) 干旱歷時

(b) 干旱烈度圖5 干旱事件匹配后的干旱特征箱形圖Fig.5 Box chart of drought characteristics of drought events after matching

表5為各流域氣象干旱傳播率和水文干旱敏感度。從干旱傳播率的結果可以看出，武隆流域和梧州流域的氣象干旱向水文干旱傳播的概率較大，超過80%；而三匯流域氣象干旱傳播率較低，僅為39%。此外，五通橋流域和萬縣流域氣象干旱傳播率也低于江邊街流域、允景洪流域和元江流域。水文干旱敏感度分布與氣象干旱傳播率分布較為一致，允景洪流域的水文干旱對氣象干旱最敏感，三匯流域的水文干旱對氣象干旱最不敏感。

表5 氣象干旱傳播率和水文干旱敏感度Table 5 Drought translation rate of meteorolgical drougt and sensitivity of hydrological drought

3.4 氣象干旱和水文干旱特征擬合及傳播閾值

為了建立氣象干旱和水文干旱之間的關系，根據干旱匹配結果計算得到氣象干旱和水文干旱事件的干旱歷時和干旱烈度，利用線性擬合方程建立干旱傳播模型描述水文干旱與氣象干旱特征因子之間的關系。方程的基本形式為

y=ax+b

(3)

式中：y為水文干旱歷時或烈度；x為氣象干旱歷時或烈度；a、b為線性擬合參數。

選取各流域前70%(取整)干旱事件來擬合方程，得到方程及其擬合優度(R2)，選取各流域其余的干旱事件進行驗證，得到方程納什效率系數(NSE)。擬合結果如表6所示，可見擬合結果合理。

表6 各流域干旱特征擬合方程Table 6 The fitting equations of drought characteristics in each basin

在氣象干旱向水文干旱傳播的過程中，已有研究表明干旱歷時具有延長的特征[28-29]。本文定義當發生1月氣象干旱時，由此觸發的水文干旱的干旱歷時為干旱歷時傳播閾值(Dt)。對于干旱傳播的閾值，更加關注干旱烈度，因此定義當發生水文干旱時，引起此水文干旱的氣象干旱的干旱烈度為干旱烈度傳播閾值(St)。根據表6的干旱特征擬合方程，進一步分析干旱特征傳播閾值，結果如表7所示。當發生1月歷時的氣象干旱事件時，由此觸發的水文干旱歷時分別為0.34月、1.58月、2.72月、2.04月、2.11月、2.62月、1.75月和2.82月。在分別發生干旱烈度為1.67、1.26、0.06、-0.58、-0.52、0.41、0.92和-0.60(負值表示只要發生氣象干旱事件，水文干旱事件的干旱烈度大于0.5)的氣象干旱時，會觸發干旱烈度為0.5的水文干旱事件。結合表5計算結果，以梧州流域為例，如果發生持續時間為1月、干旱烈度超過0.92的氣象干旱，有81%的概率會導致干旱歷時為1.75月、干旱烈度為0.5的水文干旱。此外，較大的Dt值對應的St值均較小，即氣象干旱烈度觸發閾值較小的流域，水文干旱事件歷時更長。因此，在更可能觸發水文干旱的流域，水文干旱的預防將更加重要。

表7 各流域干旱歷時和干旱烈度傳播閾值Table 7 Drought duration and propagation threshold of drought intensity in each basin

4 結 論

a.西南地區氣象干旱事件的干旱歷時和干旱烈度均較小且較為集中，而水文干旱事件的干旱歷時和干旱烈度整體較為分散且數值較大，水文干旱嚴重性強于氣象干旱。由于各流域氣象因子和下墊面因子存在差異，水文干旱對氣象干旱的響應時間也存在明顯的差異性，西南地區干旱響應時間主要為2～7月。

b.西南地區水文干旱敏感度分布與氣象干旱傳播率分布較為一致，表明西南地區氣象干旱極易發展成為水文干旱，且大部分的水文干旱由氣象干旱觸發。元江流域、江邊街流域、三匯流域、允景洪流域、武隆流域、五通橋流域、梧州流域和萬縣流域干旱歷時傳播閾值分別為0.34月、1.58月、2.72月、2.04月、2.11月、2.62月、1.75月和2.82月，干旱烈度傳播閾值分別為1.67、1.26、0.06、-0.58、-0.52、0.41、0.92和-0.60。結合干旱特征傳播閾值(干旱歷時傳播閾值和干旱烈度傳播閾值)及氣象干旱傳播率，可定量分析西南地區水文干旱對氣象干旱的響應關系，為流域水文干旱早期預警提供參考。