明清時期渭南地區旱災的多尺度變化特征及其對氣候變化的響應

郁耀闖，王長燕*，龔蘭蘭，衛旭東，馬佳俊，王林剛

(1.寶雞文理學院 地理與環境學院，陜西 寶雞 721013；2.陜西省災害監測與機理模擬實驗室，陜西 寶雞 721013)

1 引言

旱災是我國和世界上歷史時期以來頻繁發生的主要自然災害之一[1]，已對全球的農業生產、社會經濟和政治等造成了嚴重的影響和危害[2-3]。在全球變暖背景下，區域干旱災害風險呈增加趨勢[4-5]。因此，理解區域干旱災害的發生機制及其與氣候變化的關系對災害風險管理和氣候變化適應與減緩措施的制定是非常重要的。

渭南地區位于黃河中游，陜西省關中平原東部，是典型的季風邊緣區，受季風活動影響，該區氣候類型屬暖溫帶半濕潤半干旱季風氣候，區域降水空間變率較大，容易發生干旱災害。國內學者對我國西北[6]、關中[7]和寶雞地區[8]等區域旱澇災害的時空分布特征和重建進行了大量研究，得出了區域旱澇災害在時間尺度上存在著周期規律等的重要認識，但對區域干旱災害在時間尺度上的多尺度變化特征及其與氣候變化關系的研究還鮮見報道，這不利于揭示區域干旱災害的發生機制及其與氣候變化的聯系。

本研究選取我國明清小冰期(1368～1911年)作為主要研究時段，以位于季風邊緣區的渭南地區為研究案例區，搜集整理明清時期渭南地區的歷史文獻資料，利用集合經驗模態分解(EEMD)等方法，探討該時期渭南地區旱災的多尺度變化特征及其與氣候變化的關系，以期為區域干旱災害的預測和防治提供參考依據。

2 研究區域與研究方法

2.1 研究區概況

渭南位于黃河中游，地處關中平原陜西省東部，位于東經108°58′～110°35′，北緯34°13′～35°52′之間，面積1.3×104km2，2019年末總人口527.81萬人[9]。目前，轄臨渭、華州2個區，有潼關、大荔、合陽、澄城、白水、蒲城、富平7個縣，代管韓城、華陰2個省轄市(圖1)。地形上南依秦嶺，北靠黃龍山和橋山，南北高，中間低，東西開闊，呈仰瓦狀。海拔在330～2645 m，屬暖溫帶半濕潤半干旱季風氣候。轄區內各地全年平均氣溫12.7～15.6 ℃，區域降水時空分布差異較大，全年降水量390.7～592.2 mm，區域年降水變率較大，容易發生旱澇災害。

圖1 渭南政區圖

2.2 數據來源與方法

本文的干旱災害數據主要來源于《中國三千年氣象記錄總集》[10]、《中國氣象災害大典·陜西卷》[11]和渭南志等歷史文獻資料。

鑒于歷史文獻資料來源的龐雜性、記錄的主觀性和空間地點的難以定量性，本研究采用區域旱災空間分布代表區域旱災規模的方法來進行旱災的統計分析[12]。對于某個區域而言，通常選取10個具有代表性的地點進行統計分析[12]。例如，若本年度該區域只有一個點發生旱災則記為1，那么10個地點的平均值為0.1；若有2個點發生旱災，則記為2，10個地點的平均值為0.2，……，以此類推，沒有旱災記載的年份則記為0。平均值大于0.3可視為區域內旱災空間規模較大[12]。利用上述方法，對明清時期渭南地區的渭南、潼關、大荔、合陽、澄城、白水、蒲城、富平、華陰和華縣等10個縣(圖1)的旱災記錄進行統計分析。采用集合經驗模態分解 (EEMD)和Morlet小波函數的方法對該區的旱災頻次序列進行多尺度周期分解[13]。

3 結果與分析

3.1 渭南地區旱災的規模、頻次和階段性

明清時期的544 a間，渭南地區10個地點共有旱災記錄1749次，占321.5%，其中，以全區域(10個點)旱災記錄頻次的空間規模為主，10個點的旱災記錄共有154次，占區域旱災記錄總數的88.1%，5～9個點的旱災記錄有24次，占區域旱災記錄總數的10.6%，4個點以下的旱災記錄共有53次，占區域旱災記錄總數的1.3%。

由圖2可知，在14世紀中期、15世紀70年代末期至16世紀60年代、17世紀30年代末期、18世紀20年代初期至19世紀前半葉和19世紀70年代末期渭南地區旱災發生較為頻繁，多為全區域(10個點)的特大旱災，如1637～1643年的明朝崇禎大旱、1719～1723年的清朝康雍奇災和1875～1879年的光緒丁戊奇荒等歷史極端干旱事件在渭南地區均有較好記錄。

明清時期渭南地區旱災總體上表現為先下降后上升的變化趨勢(圖2)。由圖2中的11 a滑動平均曲線可以看出，渭南地區旱災具有3個明顯不同的變化階段，其中，1368～1486年和1684～1911年為該區域旱災的增加階段，增加幅度分別為146.7%和11.0%；1487～1683年為該區域旱災的減少階段，減少幅度為65%。

圖2 明清時期渭南地區旱災頻次序列

3.2 渭南地區旱災的多尺度變化特征

集合經驗模態分解(EEMD)所得各IMF分量序列主要指示各分量在不同時間尺度上的周期信號，IMF 各分量振幅的大小分別代表相應周期振蕩信號的強弱程度，振幅越大表明相應時間尺度上周期震蕩信號越強；反之，則表明相應時間尺度上周期震蕩信號越弱[14]。EEMD 分解明清時期渭南地區旱災頻次序列共得到9個IMF分量，分別表示該區旱災頻次序列在不同的年際、年代際和百年際尺度上的周期變化，具體如下。

(1)在年際尺度上，IMF1、IMF2和IMF3的平均周期分別為準3 a、準6 a和準8 a周期，方差貢獻率分別為43.6%、12.9%和8.45%，說明這3個周期信號對明清時期渭南地區旱災的影響較大，其中，IMF1 分量在1414～1474年、1480～1569年、1599～1649年、1686～1734年、1747～1814年和1846～1900年波動的振幅較大，說明該周期信號在這6個階段震蕩較為強烈，對渭南地區旱災發生的影響較大；該周期震蕩信號在1385～1413年和1649～1685年波動的振幅較小，說明該周期信號在這2個階段震蕩較弱，對渭南地區旱災發生的影響較小。IMF2 分量的振幅在1508～1568年和1853～1895年波動較大，說明該周期信號在這2個階段對該區旱災發生的影響較大；該周期震蕩信號在1389～1412年、1437～1452年、1653～1687年和1802～1831年波動的振幅較小，說明該周期信號在這4個階段震蕩較弱，對旱災發生的影響較小。IMF3 分量的振幅在1417～1435年、1479～1508年、1583～1628年、1691～1721年、1748～1785年和1836～1900年波動較大，說明該周期信號在這6個時期對旱災發生的影響較大；該周期震蕩信號在1392～1410年、1439～1454年和1647～1682年波動的振幅較小，說明該周期信號在這3個階段震蕩較弱，對旱災發生的影響較小。

(2)在年代際尺度上，IMF4、IMF5和IMF6的平均周期分別為準17 a、準28 a和準84 a周期，方差貢獻率分別為4.92%、3.21%和3.04%，說明這3個周期信號與明清時期渭南地區的干旱災害存在著某種關聯，其中，IMF4 分量在1476～1508年、1555～1602年和1735～1749年波動的振幅較大，說明該周期信號在這3個階段震蕩較為強烈，對旱災發生的影響較大；該周期震蕩信號在1397～1463年、1649～1685年和1856～1911年波動的振幅較小，說明該周期信號在這3個階段震蕩較弱，對旱災的發生影響較小。IMF5 分量的周期信號在1529～1640年波動的振幅較大，除此之外，在其他時段該周期信號振幅均較小，說明該周期信號可能只在1529～1640年間對渭南地區旱災產生影響。IMF6 分量的周期信號在1368～1446年和1635～1692年波動的振幅較大，在其他時段振幅均較小，說明該周期信號在這2個時期對旱災的發生影響較大。

(3)在百年際尺度上，IMF7和IMF8的平均周期分別為準181 a和準363 a周期，方差貢獻率分別為2.20%和0.43%，總體來看，這2個分量的周期信號振幅均較小，說明百年際尺度周期信號對明清時期渭南地區旱災的影響相對較小。

(4)RES為趨勢項，主要指示明清時期渭南地區旱災在近544 a的波動趨勢，從RES曲線的趨勢可以看出(圖3)，渭南地區在近544 a的旱災變化趨勢呈先降低后增加的趨勢，降低和增加的時間拐點出現在1685年前后。

圖3 明清時期渭南地區旱災頻次序列的集合經驗模態分解結果

3.3 明清時期渭南地區旱災對氣候變化的響應

區域旱澇災害與氣候變化具有明顯的相關關系[8,15]。干旱災害發生時，氣候通常較為寒冷干燥，降水較少；洪澇災害發生時，氣候通常較為溫暖濕潤，降水較多。葛全勝等[15]研究表明：在區域尺度上，過去2000 a以來溫度和干濕變化對降水產生了重要影響，明清時期尤為顯著。在該時期，北半球已經進入寒冷階段，但受區域位置、地理環境和季風活動強弱變化等影響，我國氣候變化呈現出了不穩定的狀態，從而導致該時期旱澇災害在全國范圍內表現出了明顯的區域差異性。14世紀后半葉至20世紀初，全球氣候變化處于“小冰期”階段，此時，我國氣候較為干冷[16-17]，并且在小尺度范圍內存在區域差異，區域尺度上氣候寒冷期降水較少，往往導致旱災發生較為頻繁[16-18]。王紹武等[19]研究表明，明清小冰期我國氣候變化響應比較強烈，并存在較大區域差異。葛全勝等[15,20]認為15世紀、17世紀和19世紀為小冰期中比較偏冷的3個時期，在這3個時期內旱災發生較為頻繁。朱士光等[16]研究認為，15世紀后期、16世紀80年代至17世紀30年代、17世紀后半葉及19世紀后半葉關中地區氣候寒冷干燥。萬紅蓮等[8]研究表明，明清時期寶雞地區旱澇災害的交替出現是對該時期關中地區氣候冷暖干濕交替變化和全球氣候變化的響應。Yu等[12]研究認為，氣候溫暖時期太陽活動較為頻繁，涇河中游地區旱澇災害多發；反之，旱澇災害發生較少。由此可見，區域旱澇災害的發生與氣候變化存在著一定關聯。

從渭南地區干旱災害頻次的11 a滑動平均曲線來看(圖2)，15世紀中葉、15世紀70年代末期至16世紀60年代、17世紀前半葉、18世紀初期至19世紀前半葉和19世紀后半葉渭南地區旱災發生次數較多，尤其是17、18世紀前半葉和19世紀后半葉該區干旱災害發生頻次較高，而17世紀中期至18世紀初期該區干旱災害發生頻率較低，這可能與該時期太陽活動的減弱有關[12]。渭南地區旱災頻次序列與我國同時期冬季溫度的變化總體上較為一致，與該時期北半球溫度的變化基本同向，但不完全同步(圖4)，這可能與渭南地區所處的地理位置、自然環境、季風活動強弱和海氣震蕩脅迫等因素有關。旱災的發生主要取決于區域降水量的多少和降水的集中程度，而降水量的多寡往往會受到區域氣溫高低的影響[16]。當氣溫較高時，區域氣候波動較大，降水減少，往往會導致旱災發生。明清時期渭南地區旱災發生的頻次序列與小冰期氣候的冷暖波動具有較好的對應關系，氣候溫暖時期，太陽活動相對較為頻繁，旱災發生相對較多，這一認識與葛全勝[15,20]和Yu等[12]研究結果較為一致，說明渭南地區明清時期旱災的發生是對全球氣候變化的響應(圖4)。

圖4 渭南地區旱災頻次序列與太陽活動[21-23]、重建溫度[24,25]、ENSO[26]和PDO[27]的比較

4 討論

太陽活動是發生在太陽大氣層和對流層中的各種復雜物理過程，也是引起地球系統氣候變化的主要能量來源[28]。太陽活動的強弱變化會引起大氣中溫度的升高或降低，導致大氣環流模式的改變，從而影響區域降水空間分布的多寡，進而導致區域旱澇災害的發生[12,29]。太陽黑子數和太陽總輻照是表征太陽活動水平強弱的2個重要參數[30]。Le Mou?l等[31]研究表明，太陽黑子活動主要通過熱帶和副熱帶太平洋地區海氣耦合交互機制來影響全球溫度和降水的變化。葛全勝等[29]研究認為，當太陽輻射處于高值期時, 我國關中平原等地降水偏少，旱災多發。Li等[32]得出了太陽黑子活動具有11和112 a顯著周期的重要認識。郝志新等[7]研究認為，關中地區960～2010年干濕變化的準10 a周期與太陽活動相對應。本研究中，集合經驗模態分解(EEMD)得到的IMF3分量中的10.4 a和11.6 a周期與太陽黑子活動的11年顯著周期較為接近，說明該區干旱災害的發生可能受到太陽活動的影響，也為太陽活動影響旱災發生提供了有力證據，該研究結果與趙新華等[32-33]的研究結果較為一致。

集合經驗模態分解(EEMD)和功率譜分析表明：明清時期渭南地區旱災頻次序列中具有準3 a、準6 a和準8 a的年際變化周期，這與厄爾尼諾(ENSO)的2～8 a周期[34]具有較好的對應關系，說明ENSO可能是影響該區明清時期旱災發生年際尺度變化的重要原因之一，同時也說明赤道太平洋海表溫度可能對該時期渭南地區的降水模式產生了重要影響，進而導致該區旱災的發生。前人關于厄爾尼諾對我國季風區降水模式的影響進行了大量的研究[35-36]，認為ENSO對全球和東亞地區年際氣候變化產生了重要的影響，本研究結果也印證了這一觀點。如圖4所示，當ENSO處于暖相位時，渭南地區的旱災發生較為頻繁，但由于目前ENSO指數重建和旱災頻次序列重建還存在著很大的不確定性，因此，ENSO指數和旱災頻次序列的定量關系還難以確定。

太平洋年代際振蕩(PDO)是北太平洋海溫年代際循環變化的一種海洋現象，通常具有20～30 a的變化周期[37]，不僅對北美乃至整個北半球大氣環流產生重要影響[37]，還可以通過調制厄爾尼諾事件來間接影響區域氣候變化[38]，進而導致旱澇災害的發生。與厄爾尼諾類似，太平洋年代際振蕩通過海洋表面溫度的高低變化影響著全球和東亞地區年代際尺度上的氣候變化。Yang等[39]研究表明，當PDO處于暖位相時，華北對應高溫、少雨和干旱時段，長江中下游多雨；反之亦然。裴琳等[40]研究表明，1600～2012年來中國東部旱澇變率與PDO序列之間具有顯著的正相關關系，PDO的冷、暖位相分別對應著華北地區的偏旱和偏澇。郝志新等[7]研究認為，PDO與關中平原的干濕變化從長時間尺度來看存在著不穩定性。本研究中，集合經驗模態分解(EEMD)表明：明清時期渭南地區旱災頻次序列IMF分量中具有的準17 a和28 a的變化周期與PDO的20～30 a周期較為接近，說明PDO可能是影響該區明清時期旱災發生年代際尺度變化的因素。如圖4所示，當PDO處于暖相位時，渭南地區旱災發生比較頻繁，這一認識與Yang等[39]的認識較為一致，與裴琳等[40]的認識相悖。綜合可知，太平洋年代際振蕩與區域干旱災害(或干濕變化)的關系還存在著很大的不確定性，二者的定量關系還有待于進一步探討。

5 結論

(1) 明清時期渭南地區旱災總體上表現為先下降后上升的變化趨勢，以全區域的特大旱災為主，大致可分為1368～1486年、1487～1683年和1684～1911年3個明顯不同的變化階段，明朝崇禎大旱、清朝康雍奇災和光緒丁戊奇荒等歷史時期重大旱災事件在渭南地區均有較好記錄。

(2)渭南地區旱災在時間尺度上存在著顯著的準3 a、準6 a和準8 a的年際周期，準11 a、準17 a、準28 a和準84 a的年代際周期和準181 a和準363 a的百年際周期，分別與厄爾尼諾的2～8 a周期、太陽活動的11 a周期和太平洋年代際振蕩的20～30 a周期相對應，說明厄爾尼諾、太陽活動和太平洋年代際振蕩可能是導致該區旱災發生的重要原因。

(3)明清時期渭南地區旱災的年際、年代際和百年際尺度周期分別可以解釋旱災方差變量的64.95%、11.17%和2.63%，說明該區旱災在年際尺度上的振蕩信號可能要強于在年代際和百年際尺度上的振蕩信號。

(4)渭南地區旱災頻次序列與我國同時期冬季溫度的變化序列總體上較為一致，與該時期的北半球溫度變化基本同向，但不同步，說明該區旱災的發生是對全球氣候變化的響應。