陳紅光 ，孟凡浩* ，薩楚拉 ，羅 敏 ，王牧蘭

(1.內(nèi)蒙古師范大學(xué)地理科學(xué)學(xué)院，內(nèi)蒙古 呼和浩特 010022；2.內(nèi)蒙古自治區(qū)遙感與地理信息系統(tǒng)重點實驗室，內(nèi)蒙古 呼和浩特 010022；3.內(nèi)蒙古自治區(qū)蒙古高原草原災(zāi)害與生態(tài)安全重點實驗室，內(nèi)蒙古 呼和浩特 010022)

以全球變暖為主的背景下，極端降水頻發(fā)、徑流銳減以及植被退化等問題日益突出，對生態(tài)系統(tǒng)極為脆弱的干旱、半干旱區(qū)草原內(nèi)陸河流域的影響不言而喻[1]。歸一化植被指數(shù)(Normalized Difference Vegetation Index，NDVI)是能夠反映地表植被覆蓋和生長狀況的有效指標(biāo)，通常用來表征植被活動的強弱[2]。近些年，在氣候變化和人類活動的影響下，全球植被呈顯著的增長趨勢。然而，植被變化不僅能影響氣候，還可以通過蒸騰蒸發(fā)，故而影響流域水文循環(huán)過程[2]。水資源緊缺及草場退化對草原流域生態(tài)系統(tǒng)安全及可持續(xù)發(fā)展會產(chǎn)生重大影響。因此，了解草原區(qū)內(nèi)陸河流域氣候、水文變化特征，并揭示干旱、半干旱區(qū)內(nèi)陸河流域植被及水資源對氣候變化的響應(yīng)程度，能夠為生態(tài)保護和生態(tài)水文變化分析提供理論依據(jù)。

近年來，許多中國學(xué)者在部分典型大流域開展了關(guān)于植被及徑流變化特征的研究[3-4]。 如李帥等[5]、劉酌希等[6]、安晨等[7]分別定量分析了氣候變化及人類活動對流域徑流的影響程度，發(fā)現(xiàn)氣候變化是主控因子。上述研究涉及內(nèi)容雖廣泛，但對流域徑流量演變趨勢分析的研究僅分析時空變化，且對其影響因素探討局限在降水氣溫以及土地利用變化等，未考慮氣候因素、生態(tài)環(huán)境演變與徑流之間的關(guān)系[8]。并且已有研究大多集中在大型流域及濕潤區(qū)，很少有學(xué)者關(guān)注北方缺乏資料的干旱、半干旱區(qū)內(nèi)陸河流[9]。鑒于此，本研究采用SWAT水文模型，結(jié)合數(shù)理統(tǒng)計法，研究變化背景下植被演變特征及其對徑流的影響，以期為流域水資源開發(fā)與生態(tài)保護提供參考與借鑒。

烏拉蓋河流域位于錫林郭勒盟草原(中國大草原之一)，是內(nèi)蒙古最大的內(nèi)陸河流域，也是內(nèi)蒙古自治區(qū)重要的畜牧業(yè)基地[10]。烏拉蓋河流域地處于氣候變化敏感區(qū)，氣候變化對該區(qū)植被生態(tài)、徑流量有顯著的影響。草場退化及水資源變化嚴(yán)重影響流域生態(tài)健康及人民的生產(chǎn)生活。然而，目前對烏拉蓋河流域的氣候和植被變化以及徑流演變趨勢的研究鮮有報道。因此，本研究以烏拉蓋河流域為研究區(qū)，分析近40年氣候變化背景下烏拉蓋河流域植被響應(yīng)特征及其對徑流量的影響。研究結(jié)果能夠為干旱、半干旱區(qū)草原內(nèi)陸河流域生態(tài)保護和生態(tài)水文變化分析提供理論依據(jù)。

1 材料與方法

1.1 研究區(qū)概況

烏拉蓋河流域位于內(nèi)蒙古錫林郭勒盟東部，地處116°20′～119°59′E，44°02′～46°42′N，是內(nèi)蒙古最大的內(nèi)陸河流域，流域面積為38 549.25 km2，平均海拔為990 m，相對高程差約為1 100 m(圖1)。烏拉蓋河發(fā)源于大興安嶺西側(cè)的寶格達山，橫貫內(nèi)蒙古烏珠穆沁草原，流域內(nèi)形成大片濕地和湖泊，最終注入烏珠穆沁盆地最低處索林淖爾洼地[11]。烏拉蓋河流域氣候類型為干旱、半干旱大陸性氣候，具有降水量少、蒸發(fā)量多的氣候特征。根據(jù)1982—2020年研究區(qū)附近氣象資料統(tǒng)計，年平均降水量由東向西遞減，約為250～400 mm，年平均溫度約為-0.9 ℃。研究區(qū)東北部為草甸草原植被，西南部為典型草原，中部為兩者過渡地帶。為了明晰氣候變化對烏拉蓋河流域不同河段水文及植被的影響，以烏拉蓋水庫和胡稍廟為結(jié)點[10]，將烏拉蓋河流域分為上游、中游和下游3段。

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

為構(gòu)建SWAT水文模型，收集了研究區(qū)的數(shù)字高程模型、5個時期土地利用數(shù)據(jù)(1980/1990/2000/2010/2020年)、土壤數(shù)據(jù)以及氣象站點數(shù)據(jù)(逐日降水、最高/低氣溫、風(fēng)速、相對濕度)。實測徑流數(shù)據(jù)由1981—2003年奴乃廟水文站和2004—2012年烏拉蓋水庫入庫水文站合成的逐月徑流量資料。模型驗證的積雪數(shù)據(jù)和蒸散發(fā)數(shù)據(jù)分別選擇2000—2012年的MOD10A1和MOD16A2產(chǎn)品，時空分辨率為8 d/500 m。歸一化植被指數(shù)NDVI為1982—2020年中國5 km分辨率逐月NDVI數(shù)據(jù)集[12]，根據(jù)研究區(qū)域空間范圍經(jīng)進一步處理后得到年尺度的NDVI數(shù)據(jù)。

1.3 研究方法

1.3.1SWAT模型建模及精度評定

由美國農(nóng)業(yè)部和農(nóng)業(yè)研究局開發(fā)的SWAT水文模型，能較好的模擬變化環(huán)境下的水文循環(huán)過程[7,13]。因此，本研究以站點氣象數(shù)據(jù)以及遙感數(shù)據(jù)為基礎(chǔ)，構(gòu)建SWAT模型分析氣候變化以及植被變化背景下的烏拉蓋河流域不同河段蒸散發(fā)及水文循環(huán)過程。除徑流外，本研究選取MODIS積雪和蒸散發(fā)等數(shù)據(jù)對模型模擬結(jié)果進行多目標(biāo)驗證，并采用納什效率系數(shù)NSE、百分比偏差PBIAS、判定系數(shù)R2對模型模擬精度進行評價。評價指標(biāo)公式如下：

(1)

(2)

R2=

(3)

式中Qiobs——實測值，m3/s；Qisim——模擬值，m3/s；Qmean——實測均值，m3/s；n——實測數(shù)據(jù)量[14]。

1.3.2線性趨勢法

線性趨勢法是對研究時段內(nèi)的時間序列資料表現(xiàn)為接近于一條直線，表現(xiàn)為近似直線的上升和下降時采用的一種預(yù)測方法。其優(yōu)點是可以定性分析趨勢變化和檢驗序列的顯著性。本研究以線性趨勢法為主，以多年滑動平均法為輔助，分析烏拉蓋河流域氣象水文植被因素的多尺度變化特征。

1.3.3克里金插值法

克里金法原理是利用變量為基礎(chǔ)，對未知樣點進行線性無偏、最優(yōu)化估計[15]。被廣泛應(yīng)用于各類觀測的空間插值，如地質(zhì)學(xué)、環(huán)境科學(xué)以及近地面氣溫、降水、風(fēng)場等的單點觀測中。本研究將此方法應(yīng)用于站點氣象數(shù)據(jù)從而得到氣溫降水空間分布圖。

1.3.4M-K趨勢檢驗法

Mann-Kendall趨勢檢驗法(簡稱M-K法)是由Mann和Kendall提出并改進的統(tǒng)計方法。M-K檢驗法的優(yōu)點為對樣本數(shù)據(jù)無要求，其受異常值的干擾較小，計算步驟簡單，是目前長時間序列分析中應(yīng)用較多的一種趨勢檢驗方法[16]。當(dāng)M-K檢驗突變點時，檢驗變量S計算公式如下：

(4)

(5)

(6)

秩序Sk是第i時刻數(shù)值大于j時刻數(shù)值個數(shù)的累計值。通過分析曲線UBK和UFK可以表明要素變化趨勢以及發(fā)生的突變年份[17]。若UBK和UFK的交點位于臨界線內(nèi)，該交點為突變點。

2 結(jié)果與分析

2.1 SWAT模型校準(zhǔn)與驗證

根據(jù)站點實測徑流數(shù)據(jù)與模擬徑流數(shù)據(jù)對比來看，徑流模擬率定期及驗證期評價指標(biāo)NSE、R2和PBLAS分別為0.62和0.71、0.62和0.72、18.8%和8.5%。此外，通過對比MODIS數(shù)據(jù)的積雪面積、蒸散發(fā)與模型模擬的積雪面積、蒸散發(fā)，結(jié)果發(fā)現(xiàn)，模擬值與遙感數(shù)據(jù)擬合較好，且兩者的相關(guān)系數(shù)分均高于0.8(圖2)，可用于下一步的研究。

2.2 烏拉蓋河流域氣候變化特征

2.2.1烏拉蓋河流域降水和溫度變化特征

烏拉蓋河流域降水量和氣溫的趨勢變化特征見圖3a、3b。烏拉蓋河流域降水量以8.9 mm/10a的速率呈不顯著減少，氣溫以0.43 ℃/10a的速率呈顯著增加(p<0.01)，說明近40 a烏拉蓋河流域呈暖干化趨勢。此外，降水量變化特征顯示90年代波動較大，此后10 a降水量一致處于較低狀態(tài)(圖3a)。氣溫呈顯著的上升趨勢，這與全球氣候變化趨勢一致[8]。通過分析M-K檢驗結(jié)果可以發(fā)現(xiàn)，降水及氣溫的UF和UB 2條曲線在1994年出現(xiàn)交點，且通過顯著性檢驗(圖3c、3d)，表明在1994年烏拉蓋河流域降水和氣溫發(fā)生突變。

圖4a、4c分別給出了烏拉蓋河流域不同河段降水和氣溫變化趨勢。不難發(fā)現(xiàn)，近40 a流域上、下游降水量均表現(xiàn)為減少趨勢，而中游與其變化趨勢相反(圖4a)。由圖4b可以看出近40 a溫度在不同河段表現(xiàn)為較一致的增加趨勢，并且增溫幅度相差較小。烏拉蓋河流域降水及氣溫分布具有明顯的空間異質(zhì)性。流域降水量由東向西逐漸增加。在流域上游區(qū)年均降水量達到399～457 mm(圖4c)。而溫度由流域西南向東北遞減，整個流域中上游溫度最低(圖4d)。

2.2.2烏拉蓋河流域蒸散變化特征

蒸散發(fā)是水量平衡與能量平衡的重要組成部分，在水分運動過程中占有重要地位，對全球及區(qū)域水循環(huán)和能量循環(huán)過程中有重要意義。圖5顯示了由烏拉蓋河流域的SWAT模型中提取的實際蒸散發(fā)年際變化特征?？梢钥闯?，整個流域年均蒸散發(fā)約為200 mm，且呈不顯著的上升趨勢。從不同河段來看，流域上游蒸散量最小，最低值只有122 mm，且以2.89 mm/10a的速率減少。流域上游蒸散呈減少趨勢，而中、下游蒸散呈增加趨勢，增速分別為1.15、3.68 mm/10a。各河段氣溫和降水變化速率以及不同的植被類型均可導(dǎo)致蒸散量變化趨勢的不同[8]。

2.3 流域氣候變化對植被的影響

處于干旱、半干旱區(qū)的烏拉蓋河流域主要以草地、林地和未利用地為主，并且對氣候變化的響應(yīng)十分敏感。分析暖干化背景下該流域植被變化對區(qū)域生態(tài)系統(tǒng)可持續(xù)發(fā)展尤為重要。近40 a烏拉蓋河流域年均NDVI指數(shù)(歸一化化植被指數(shù))約為0.29，最大值和最小值分別出現(xiàn)在2018、2007年(圖6a)?？傮w上，1982—2020年年均NDVI以0.005/10a的速率呈不顯著的上升趨勢，表明烏拉蓋河流域植被恢復(fù)較明顯。從不同河段來看，變化速率相差不大，但平均NDVI指數(shù)卻不盡相同。上游NDVI在0.32～0.44波動，而下游在0.22～0.32變化。由烏拉蓋河流域的NDVI指數(shù)的空間分布(圖6b)發(fā)現(xiàn)，流域上、中游NDVI指數(shù)較高，流域下游南部地區(qū)，也就是烏拉蓋戈壁NDVI指數(shù)相對較低。

1982—2020年烏拉蓋河流域NDVI值數(shù)與降水、蒸散發(fā)、溫度均成正相關(guān)。不同河段相關(guān)系數(shù)顯示溫度及蒸散發(fā)對植被影響一致，而降水量具有明顯的區(qū)域差異，流域上游降水與植被成正相關(guān)，而中、下游成負相關(guān)(表1)。整個流域降水與植被變化的相關(guān)性最強，溫度呈最小(圖7)。而3個河段中流域蒸散發(fā)與NDVI值數(shù)的相關(guān)性呈最強。

表1 烏拉蓋河流域及各河段NDVI值數(shù)與降水量、溫度及蒸散的相關(guān)系數(shù)

注：*為通過0.1信度的顯著檢驗；**為通過0.05信度的顯著檢驗；***為通過0.01信度的顯著檢驗。

不難發(fā)現(xiàn)，流域氣候呈暖干化趨勢下，流域植被恢復(fù)較明顯。但植被變化一方面受氣候變化因素的影響，另一方面由于過度開墾以及過度放牧等人為因素造成流域內(nèi)草場退化也有較大影響[8]，且近年來烏拉蓋河流域植被好轉(zhuǎn)與多年是實施的一系列生態(tài)保護措施，包括“退耕還林還草”“京津風(fēng)沙源治理”“未封禁牧”[18]等政策的影響也不容忽視。

2.4 流域氣候及植被變化對徑流的影響

根據(jù)SWAT模擬結(jié)果可以看出，烏拉蓋河流域年均徑流深以11.9 mm/10a的速率呈顯著減少趨勢(p<0.05)。其中，徑流深的最大值出現(xiàn)在1998年，可達183.24 mm，最小值出現(xiàn)在2017年，為12.08 mm(圖8a)。且不同河段徑流變化趨勢一致，徑流值最大值均出現(xiàn)在1998年。通過分析M-K檢驗結(jié)果可以發(fā)現(xiàn)，在2000年烏拉蓋河流域徑流深發(fā)生明顯突變(圖8b)。表明2000年之前烏拉蓋河流域徑流主要受氣候變化因素的影響，自突變年以后徑流不僅受氣候變化影響，還受人類活動的影響。而下墊面變化是人類活動的主要表現(xiàn)形式之一。通過對流域近40 a土地利用以及NDVI對比分析，烏拉蓋河流域植被總體呈改善趨勢，無植被覆蓋度區(qū)及水域面積逐漸減小，低覆蓋度植被、中覆蓋度植被和建設(shè)用地面積增加。植被覆蓋度的提高，一方面增加了土壤入滲，另一方面增加了植被蒸發(fā)，從而使主要以降水為主的徑流減少[8]。

由烏拉蓋河流域多年平均徑流深的空間分布(圖9)看出，流域中部徑流較大，東部徑流較小，并且由東向西呈先增大后減小趨勢，中游徑流較小，流域下游徑流較大，主要分布在烏拉蓋濕地處。

為了明晰氣候暖干化以及植被好轉(zhuǎn)的背景下烏拉蓋河流域徑流深的演變特征，本研究基于近40 a內(nèi)流域氣候和土地利用數(shù)據(jù)，應(yīng)用SWAT模型分別模擬了4個不同時期(1980—1990年、1990—2000年、2000—2010年、2010—2020年)的徑流量，4個時段的產(chǎn)流見表2。如表所示，在流域氣候暖干化以及植被恢復(fù)背景下，自1980—1990年到1990—2000年為止總產(chǎn)水量增加，到2000年時減少，而后又增加。自2000年后總產(chǎn)水量的急速減少，不僅因為流域內(nèi)氣候呈暖干化趨勢，還因為植被恢復(fù)而導(dǎo)致吸收更多的水分。除此之外，與2003年在烏拉蓋河上游建立的水庫，截取了上游水資源有較大的關(guān)系[18]。而自2010—2020年總產(chǎn)水量的增加主要因為近幾年降水較豐富，當(dāng)然，也離不開自2000年開始實施的水土保持措施的效應(yīng)[19]。

表2 烏拉蓋河流域不同時段產(chǎn)流比 單位：mm

3 結(jié)論

基于1981—2020年烏拉蓋河流域的水文氣候要素，結(jié)合SWAT水文模型，分析了近40 a烏拉蓋河流域氣候、草地生態(tài)和徑流深的變化特征，揭示了氣候變化對干旱、半干旱區(qū)植被的影響以及變化環(huán)境對內(nèi)陸河流域徑流量的影響程度。主要結(jié)論如下。

近40 a流域氣候呈暖干化趨勢明顯，降水減少(-8.9 mm/10a)，溫度顯著升高(0.43 ℃/10a)，蒸散發(fā)增加(3.3 mm/10a)。在氣候變暖背景下，流域的植被覆蓋整體呈好轉(zhuǎn)趨勢，自1982年以來流域NDVI以0.005/10a的速率呈不顯著增長。1982—2020年流域NDVI指數(shù)與氣候因子(降水、溫度、蒸散發(fā))均成正相關(guān)，且與降水量的相關(guān)性最強。在流域氣候暖干化趨勢以及草地植被恢復(fù)的背景下，流域近40 a年均徑流深呈顯著減少趨勢(-11.9 mm/10a)，并在2000年發(fā)生突變，自突變之后總徑流量減少了69.20%。流域內(nèi)降水量的減少以及氣溫的升高使徑流深減少，而植被的恢復(fù)加劇了這一變化。