肖森元，蘇 軍，楊 廣*，劉 兵，張 迪，仇澤豪，羅雅潔，黃晶晶

(1.石河子大學(xué)水利建筑工程學(xué)院/寒旱區(qū)生態(tài)水利工程兵團重點實驗室，新疆 石河子 832000;2.第八師石河子市水利工程管理服務(wù)中心，新疆 石河子 832000)

氣候變化和人類活動改變了陸地水文循環(huán)的過程，影響了水文和水資源系統(tǒng)的結(jié)構(gòu)和功能，并對水資源的開發(fā)和利用構(gòu)成重大挑戰(zhàn)[1-2]。近幾十年來，氣溫變暖加速了全球水文循環(huán)，從而影響了水循環(huán)過程并改變了盆地降水量與水循環(huán)之間的時空格局[3-4]。人類活動則進一步促進了降水和徑流的變化，導(dǎo)致流域水資源的時空分布發(fā)生了變化[5-6]，水資源供需平衡遭到破壞，干旱問題日益突出[7]。因此揭示氣候變化和人類活動影響下流域徑流和干旱演變特征成為當(dāng)前亟待解決的重要問題[8]。

目前，國內(nèi)外學(xué)者在氣候變化和人類活動對徑流和干旱的影響方面做了大量研究，并取得豐碩的成果。目前應(yīng)用較為廣泛的主要有數(shù)理統(tǒng)計法、彈性系數(shù)法及水文模型法。數(shù)理統(tǒng)計法以徑流和氣候因子間的線性關(guān)系為基礎(chǔ)，分離氣候變化和人類活動對徑流影響。劉可群等[9]對洪湖春季干旱分析，指出洪湖春旱增多加重主要受氣候變化和人類活動影響，而人類活動產(chǎn)生的影響更顯著。Alifujiang Y[10]等采用雙累積曲線對伊塞克湖盆地徑流量變化進行等量分析，發(fā)現(xiàn)氣候變化和人類活動對徑流減少的貢獻率分別為64.7%和35.3%，氣候變化是導(dǎo)致伊塞克湖盆地徑流減少的主要原因。彈性系數(shù)法主要基于水熱耦合平衡方程，其中基于Budyko假設(shè)的分析方法應(yīng)用最為廣泛。薛聯(lián)青等[11]基于Budyko 假設(shè)，對塔里木河流域徑流變化對氣候變化和人類活動的響應(yīng)進行了敏感性分析。Zheng等[12]基于Budyko假設(shè)的互補關(guān)系方法探討氣候變化和人類活動對黃河流域徑流變化的貢獻率，結(jié)果表明黃河流域徑流減少主要是由于氣候變化，其貢獻為60%～70%。水文模型是分析流域水文過程的工具，可實現(xiàn)不同氣候情景下的徑流模擬。 Seguis等[13]應(yīng)用分布式水文模型對尼日爾的薩赫勒小流域徑流進行模擬研究發(fā)現(xiàn),人類活動對徑流的影響比氣候因素對徑流的影響重要。Lehner等[14]通過采用WaterGAP模型模擬未來水流不同流態(tài)認為，未來歐洲北部到東北部在全球變化的影響下易發(fā)生洪水，而歐洲南部和東南部則易發(fā)生干旱。

瑪納斯河流域是典型內(nèi)陸河流域，地處中國西北干旱區(qū)，水文生態(tài)平衡更加脆弱敏感。近年來，隨著水資源的日益消耗，水問題變得非常尖銳，出現(xiàn)了嚴重的生態(tài)環(huán)境問題，如地下水位下降、農(nóng)田次生鹽漬化、植被大面積退化等。目前，在瑪納斯河流域展開的氣候變化和人類活動研究主要集中在對水文水資源產(chǎn)生的影響進行定量區(qū)分方面[15]，并未涉及水文干旱演變，而如何加強人類活動對瑪納斯河流域水文變化的正向引導(dǎo)作用，實現(xiàn)區(qū)域水資源可持續(xù)發(fā)展，亟待對流域水文及干旱演變規(guī)律進行研究。基于此，本研究通過Mann-Kendall非參數(shù)檢驗法、雙積累曲線法和SRI指數(shù)法對氣候變化和人類活動影響下瑪納斯河流域徑流和水文干旱進行分析，并對比分析了氣候變化和人類活動影響下流域干旱狀況，為流域水土資源合理利用提供基礎(chǔ)依據(jù)。

1 數(shù)據(jù)來源與研究方法

1.1 研究區(qū)域與數(shù)據(jù)來源

瑪納斯河流域位于新疆天山北麓、準(zhǔn)噶爾盆地南緣，地理位置為43°27′～45°21′N，85°01′～86°32′E(圖1)。流域地勢東南高西北低，自東向西分布有塔西河、瑪納斯河、寧家河、金溝河以及巴音溝河，其中瑪納斯河是流域內(nèi)流量最大的河流，年均徑流量為12.9×108m3，流域年均降雨量為115～200 mm，屬于冰雪融水和降雨混合補給型河流。由于流域地處新疆腹地，形成氣溫差異大，蒸發(fā)量大的溫帶荒漠氣候，流域年均氣溫為4.7～5.7℃，年均蒸發(fā)量為1 500～2 100 mm。

本研究數(shù)據(jù)來源于瑪納斯河流域肯斯瓦特水文站與石河子和烏蘭烏蘇氣象站提供的1956—2020年實測的逐月徑流和逐月降雨數(shù)據(jù)。

圖1 研究區(qū)位置

1.2 研究方法

a)徑流變化趨勢及變異分析。Mann-Kendall是基于秩的非參數(shù)檢驗，它對異常值不敏感，并且不要求數(shù)據(jù)呈正態(tài)分布，因此廣泛應(yīng)用于降水、溫度及徑流等氣象水文序列的變化趨勢檢驗[16]。在分析時間序列的趨勢和突變點時，首先計算樣本統(tǒng)計量UF與UB，確定時間序列的變化趨勢。當(dāng)UF超過顯著性置信度臨界線(±1.96)時，該序列在統(tǒng)計上存在顯著的增加/減少趨勢，反之則不表現(xiàn)出顯著的趨勢。若UF和UB 2條曲線的交點在臨界線內(nèi)，則交點對應(yīng)的時刻被視為時間序列的突變點。具體步驟見文獻[17-18]。

b)徑流變化歸因分析。雙累積曲線法是研究2個序列間變化關(guān)系和判別其是否具有一致性的方法，由于該方法簡單、直觀、有效而在水文氣象要素分析中得到廣泛應(yīng)用。首先，將2個相同長度的序列進行累加，并在直角坐標(biāo)系中繪制關(guān)系線，依據(jù)關(guān)系線的斜率是否發(fā)生變化出現(xiàn)拐點，判斷其是否發(fā)生異常，而出現(xiàn)的拐點對應(yīng)的時間則為發(fā)生異常的年份。以異常年份為時間點，對長度序列進行比較，若相對于變化前，變化后序列出現(xiàn)偏離，則表示變化后序列可能受外界影響出現(xiàn)顯著性變化，否則表示兩序列一致性較好，未出現(xiàn)顯著性變化。而后，將第一個突變點之前有實測資料的時段作為自然期，之后作為變化期，根據(jù)自然期內(nèi)降水-徑流累積關(guān)系，反推變化期受氣候變化影響徑流量，實現(xiàn)氣候變化和人類活動對徑流影響的分離，從而探究人類活動和氣候變化對突變期相對于基準(zhǔn)期徑流量變化的貢獻率，對徑流量變化進行歸因分析。具體步驟見文獻[19-21]。

c)水文干旱分析。標(biāo)準(zhǔn)化徑流指數(shù)(SRI)是持續(xù)時間內(nèi)累積的水文徑流百分比的單位標(biāo)準(zhǔn)正態(tài)偏差，流域內(nèi)不同時間尺度下的徑流情況和干旱程度可以通過SRI的大小反映[22-23]。其計算方法如下。

首先獲取徑流的時間序列，并對時間序列所代表的樣本進行概率分布擬合：

(1)

然后對其進行累加，得到累積概率計算公式見式(2)：

(2)

最后將累積概率分布函數(shù)轉(zhuǎn)換為標(biāo)準(zhǔn)正態(tài)分布函數(shù)：

(3)

由于標(biāo)準(zhǔn)化徑流指數(shù)(SRI)的計算方法與標(biāo)準(zhǔn)化降水指數(shù)(SPI)類似，故可參照國家氣候中心頒布的氣象干旱等級分類標(biāo)準(zhǔn)[24]確定SRI指數(shù)的干旱等級。具體等級劃分見表1。

表1 標(biāo)準(zhǔn)化徑流指數(shù)干旱等級劃分

2 結(jié)果與分析

2.1 徑流變化趨勢和變異分析

根據(jù)1956—2020年實測徑流和降水資料對瑪納斯河流域年徑流變化進行分析可知，從整體上看，研究時段內(nèi)瑪納斯河流域年徑流深呈不規(guī)則周期波動(圖2a)，其中1999年徑流深為歷年最大值，達到了435.46 mm，而1992年僅有201.78 mm，為歷年最小值。年降水量在1955—1965 年波動幅度較大且有先增加后減少的趨勢,1966—1977 年呈先波動增加后減少，在1977 年達到最小值，從1978 年起降水量呈波動增加的變化趨勢。從年際上看，65 a來瑪納斯河流域年徑流深呈現(xiàn)出比較明顯的線性增加趨勢，尤其是1995年以后遞增幅度最大。從年代上來看(圖2b)，流域年徑流深呈現(xiàn)出“先增后減再增”的變化趨勢。除20世紀(jì)60年代到70年代徑流深呈現(xiàn)出減少趨勢外，其他年代徑流深均呈現(xiàn)出不同程度遞增狀態(tài)，尤其是20世紀(jì)80年代到90年代增加最為明顯，變化幅度高達17.27%。進入21世紀(jì)后，流域年徑流變化趨勢有所減緩，10 a變化幅度約為2%。在這一時期，降雨量則呈現(xiàn)“先減后增再平穩(wěn)”的變化趨勢。

a)流域徑流

b)降雨量

利用Mann-Kendall非參數(shù)統(tǒng)計檢驗法對瑪納斯河流域年徑流分析可知(圖3)，在0.05顯著性水平下，流域年徑流在1956—2004年均未出現(xiàn)顯著性變化，在2005—2020年統(tǒng)計檢驗指標(biāo)超過臨界線，表明在此時間段內(nèi)年徑流呈顯著性增加趨勢。UF和UB 2條曲線相交的時間節(jié)點是1995年，且交點位于顯著性置信度臨界線內(nèi)，表明流域年徑流發(fā)生突變的時間是1995年。1995年恰是瑪納斯河流域推行節(jié)水灌溉技術(shù)的時期，流域規(guī)模化節(jié)水措施作為重要的人類活動，在迅速推動綠洲化進程的同時，加劇了流域徑流的年際變化趨勢。

圖3 年徑流量Mann-Kendall檢驗

2.2 徑流變化影響因素分析

2.2.1氣候變化和人類活動對徑流的影響

氣候變化和人類活動是影響地表徑流的2個重要因素，為識別不同驅(qū)動因素的貢獻，本研究利用降水-徑流雙累積曲線來分析量化氣候變化和人類活動對瑪納斯河流域徑流的貢獻比例(圖4)。

圖4 瑪納斯河流域年降水-徑流雙累積曲線

由圖4可知，降水-徑流雙累積曲線在1995年出現(xiàn)拐點，說明1995年為其確定的突變年份，這與利用Mann-Kendall非參數(shù)統(tǒng)計檢驗結(jié)論相一致。因此，依據(jù)判定的突變年份，以1995年為時間節(jié)點，將年徑流序列劃分為1956—1995年和1996—2014年2個子序列。選取1956—1995年作為自然期，認為在該時期徑流受自然環(huán)境變化影響較大，受到人類活動影響較小，選取1996—2020年作為變化期，認為在該時期徑流受氣候變化和人類活動的雙重影響。經(jīng)計算基準(zhǔn)期和變化期2個子序列的年徑流均值可知，基準(zhǔn)期流域年徑流均值為250.93 mm，變化期流域年徑流均值較基準(zhǔn)期增加62.43 mm，約增加24.88%，可見，流域徑流在變化期增加較為明顯。

2.2.2氣候變化和人類活動對徑流的貢獻

為定量分析氣候變化和人類活動對流域徑流的影響程度，以自然期月累積降水-徑流為原始數(shù)據(jù)，構(gòu)建月降水-徑流雙累積關(guān)系模型(∑R=0.77633·∑p-264.62,R2=0.99)，反推變化期受氣候變化影響徑流量，實現(xiàn)氣候變化和人類活動對徑流影響的分離，而后通過計算實測值和擬合值差值占變化值的百分比，確定其貢獻率，結(jié)果見表2。

由表2可知，氣候變化對月徑流量變化的影響量為2.95 mm，貢獻率為56.73%，人類活動對月徑流量變化的影響量為2.25 mm，貢獻率為43.37%，表明氣候變化是導(dǎo)致瑪納斯河流域徑流量變化的主要因素。

表2 氣候變化和人類活動對徑流變化的貢獻率

2.3 水文干旱演變特征分析

在進行水文干旱分析時，采用與SPI指數(shù)相類似的方法，所以SRI指數(shù)和SPI指數(shù)一樣具有多時間尺度特征，可進行不同時間尺度下的干旱特征分析?？紤]到研究區(qū)主要發(fā)展綠洲農(nóng)業(yè)，干旱事件的發(fā)生又多以季節(jié)性干旱為主，因此，在上述徑流分析的基礎(chǔ)上，利用3個月尺度SRI指數(shù)對氣候變化和人類活動影響下流域水文干旱演變進行分析。

為了便于進行不同影響條件下季節(jié)性干旱變化分析，根據(jù)研究區(qū)地理和氣候特征將12個月依次劃分為春季(3—5月)、夏季(6—8月)、秋季(9—11月)和冬季(12月至次年2月)[25]。并繪制不同影響條件下不同季節(jié)SRI變化過程，見圖5。

a)春季 b)夏季

c)夏季 d)夏季

春季，實測SRI和受氣候變化影響下的SRI在春季的變化趨勢存在較大差異，實測SRI主要呈現(xiàn)出輕度干旱化趨勢，受氣候變化影響下的SRI則呈現(xiàn)出無旱化趨勢。比較春季SRI值可知，前者共出現(xiàn)輕旱事件23次，在1996—2007年和2009—2020年出現(xiàn)連續(xù)春旱，后者僅2012年出現(xiàn)1次輕旱。說明在人類活動影響下，流域春季出現(xiàn)持續(xù)干旱化趨勢。這主要是由于瑪納斯河屬季節(jié)性河流，春季上游來水量少，而農(nóng)田播種又需要消耗大量的水資源，故使得流域春季呈現(xiàn)出干旱化趨勢。

夏季，實測SRI和受氣候變化影響下的SRI在夏季的變化趨勢大致相同，總體上均呈下降趨勢。比較夏季SRI值可以發(fā)現(xiàn)，在所取的研究時段內(nèi)，兩者均無干旱事件發(fā)生，但后者SRI值波動幅度明顯大于前者。分析原因可知，膜下滴灌技術(shù)的實施降低了作物的棵間蒸發(fā)，增加了土壤的濕度，使得流域夏季干旱化趨勢得到有效控制。

秋季，實測SRI和受氣候變化影響下的SRI變化趨勢大致相似，主要呈現(xiàn)出無旱化趨勢。比較秋季SRI值可以發(fā)現(xiàn)，前者在所取的研究時段內(nèi)，無干旱事件發(fā)生；后者處于輕旱以上共出現(xiàn)8次，分別在1996—1998年、2005—2006年、2008年、2013年以及2019年出現(xiàn)秋旱，說明人類活動減少了流域干旱事件的發(fā)生。這主要得益于流域綠洲區(qū)人工綠洲面積的增長，改善了流域的生態(tài)環(huán)境，從而使得秋季無旱化趨勢更加穩(wěn)定。

冬季，實測SRI和受氣候變化影響下的SRI值均在-0.5以下，干旱程度明顯。進一步分析比較冬季SRI值可知，前者處于中旱以上共出現(xiàn)3次，是在1996—1997年和2014年；后者處于中旱以上共出現(xiàn)18次，其中在1996—1997年、2001年、2011年和2019年出現(xiàn)了重度干旱。經(jīng)分析可知，冰川融雪水對瑪納斯河徑流有著重要的補給作用，而人類活動加速了冰川消融，提高了瑪納斯河冬季的徑流量，從而在短期內(nèi)改善了流域冬季的干旱程度。

對發(fā)生的季節(jié)性干旱事件進行頻次和頻率統(tǒng)計可知(圖6)，春季實測SRI以輕旱為主，無旱以及輕旱發(fā)生的頻次和頻率分別為2次、23次和8%和92%；春季氣候變化影響下SRI僅發(fā)生1次輕旱，其余均為無旱，發(fā)生輕旱和無旱的頻率分別為4%和96%。夏季實測SRI和氣候變化影響下SRI均無干旱事件發(fā)生，統(tǒng)計的25次均處于無旱。秋季實測SRI無干旱事件發(fā)生；秋季氣候變化影響下SRI無重旱事件發(fā)生，無旱、輕旱和中旱發(fā)生的頻次和頻率分別為17次、7次、1次和68%、28%、4%。冬季實測SRI發(fā)生3次中旱和21次輕旱，發(fā)生中旱和輕旱的頻率為12.5%和87.5%；冬季氣候變化影響下SRI發(fā)生6次輕旱、13次中旱和5次重旱，發(fā)生的頻率為25%、54.2%和20.8%。

對比分析可知，實測SRI易在春季發(fā)生輕度干旱事件，而氣候變化影響下SRI則更易在冬季發(fā)生中旱及其以上干旱事件。可見人類活動一方面使流域春季出現(xiàn)干旱化趨勢，需要人類警惕；另一方面緩解了流域冬季的干旱程度，對于改善季節(jié)性干旱起到推動作用。

a)實測

b)受氣候變化影響

3 討論

流域徑流的形成和變化不僅與降水的變化有關(guān)，而且還與人類活動引起的地表條件的變化以及人類活動的直接影響有關(guān)[26]。本研究從氣候變化和人類活動的角度出發(fā)，定量分析了這2種影響因素對流域徑流的影響程度，并通過對比分析探討了在這2種影響因素下流域水文干旱演變情況。流域徑流量的增加受氣候變化和人類活動的雙重影響，相對于人類活動所產(chǎn)生的影響，氣候變化對流域流量變化的影響量為2.95 mm，貢獻率為56.73%，而人類活動對徑流量變化的影響量為2.24 mm，貢獻率為43.37%，可見，氣候變化對流域流量變化的影響更加顯著。這與陳伏龍等[15]、高培等[27]的研究結(jié)果保持一致。同時需要指出的是，雖然降水-徑流雙累積關(guān)系模型引入了累積量以減少年際波動對測量數(shù)據(jù)的影響且年與累積之間的線性擬合具有較高的相關(guān)系數(shù)(R2=0.99)，但本研究僅將降水視為氣候因素，實際上應(yīng)包括日照，溫度，風(fēng)速和其他因素[28]，因此對氣候因素和人類活動的分解缺乏徹底性，定量研究仍存在不確定性。

Mann-Kendall非參數(shù)統(tǒng)計檢驗和降水-徑流雙累積曲線表明瑪納斯河流域的徑流量在1995年發(fā)生突變，這與流域?qū)嵭泄?jié)水灌溉的時間相符，流域內(nèi)人口和耕地面積大幅增加(圖7)，人類活動對徑流的影響加劇。另外瑪納斯河屬冰川融雪及降雨混合補給型河流，由于人類對自然資源不合理的開發(fā)利用，導(dǎo)致全球氣溫持續(xù)升高，山區(qū)冰川消融加劇。研究[29]發(fā)現(xiàn)瑪納斯河流域已有 32.5%的冰川面積消失，雖然短期內(nèi)冰川融雪水對徑流補給量的增加能夠降低流域干旱事件的發(fā)生，但從長遠來看冰川融水量先增后減的拐點一旦出現(xiàn)，以降水和冰川融水為主要補給的瑪納斯河徑流量會出現(xiàn)陡降，在枯水季節(jié)或年份可能出現(xiàn)區(qū)域性缺水危機。同時由于人類活動的影響，現(xiàn)階段瑪納斯河流域春季呈現(xiàn)出持續(xù)干旱化趨勢，針對這一現(xiàn)象，瑪納斯河流域應(yīng)合理控制開墾面積，結(jié)合研究區(qū)實際情況建議調(diào)整農(nóng)作物種植結(jié)構(gòu)，減少高耗水農(nóng)作物種植面積，發(fā)展低耗水農(nóng)作物。更為重要的是，應(yīng)大力培養(yǎng)高效農(nóng)業(yè)技術(shù)管理人員，盡管研究區(qū)高效節(jié)水技術(shù)已經(jīng)發(fā)展的相當(dāng)成熟，但相應(yīng)的技術(shù)管理水平還有待進一步提高，因此，培養(yǎng)高效農(nóng)業(yè)技術(shù)管理人員是接下來發(fā)展的首要任務(wù)。

注：橙色實虛線分別代表1975—1995、1996—2020年的平均值。

4 結(jié)論

a)近65 a來瑪納斯河流域徑流量總體上呈增加趨勢，Mann-Kendall突變檢驗法和降水-徑流雙累積曲線法識別出流域年徑流在1995年左右發(fā)生突變，2005—2020年徑流量呈顯著性增加趨勢。

b)氣候變化是引起瑪納斯河流域徑流增加的主要原因，對流域月徑流量增加的貢獻率為56.73%；人類活動主要通過流域開墾及農(nóng)業(yè)灌溉等方式間接影響徑流變化，對流域月徑流量增加的貢獻率為43.37%。

c)通過利用SRI對變化期徑流進行季尺度分析發(fā)現(xiàn)，相對于受單一氣候變化的影響，人類活動降低了流域干旱事件發(fā)生的頻率，短期內(nèi)改善了流域的干旱狀況；但人類活動同樣打破了原有的水文生態(tài)平衡，使得流域春季出現(xiàn)了干旱化趨勢。