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

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

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

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

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

1 材料與方法

1.1 研究區(qū)概況

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

1.2 數(shù)據(jù)來(lái)源

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

1.3 研究方法

1.3.1SWAT模型建模及精度評(píng)定

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

(1)

(2)

R2=

(3)

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

1.3.2線(xiàn)性趨勢(shì)法

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

1.3.3克里金插值法

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

1.3.4M-K趨勢(shì)檢驗(yàn)法

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

(4)

(5)

(6)

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

2 結(jié)果與分析

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

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

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

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

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

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

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

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

2.3 流域氣候變化對(duì)植被的影響

處于干旱、半干旱區(qū)的烏拉蓋河流域主要以草地、林地和未利用地為主，并且對(duì)氣候變化的響應(yīng)十分敏感。分析暖干化背景下該流域植被變化對(duì)區(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的速率呈不顯著的上升趨勢(shì)，表明烏拉蓋河流域植被恢復(fù)較明顯。從不同河段來(lái)看，變化速率相差不大，但平均NDVI指數(shù)卻不盡相同。上游NDVI在0.32～0.44波動(dòng)，而下游在0.22～0.32變化。由烏拉蓋河流域的NDVI指數(shù)的空間分布(圖6b)發(fā)現(xiàn)，流域上、中游NDVI指數(shù)較高，流域下游南部地區(qū)，也就是烏拉蓋戈壁NDVI指數(shù)相對(duì)較低。

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

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

注：*為通過(guò)0.1信度的顯著檢驗(yàn)；**為通過(guò)0.05信度的顯著檢驗(yàn)；***為通過(guò)0.01信度的顯著檢驗(yàn)。

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

2.4 流域氣候及植被變化對(duì)徑流的影響

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

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

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

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

3 結(jié)論

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

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