2.2.3景觀指數(shù)計(jì)算利用Fragstats軟件,從景觀類(lèi)型與斑塊水平的兩個(gè)角度計(jì)算景觀指數(shù),包括斑塊數(shù)量、斑塊形狀指數(shù)、斑塊密度、Shannon多樣性指數(shù)和Shannon均勻度指數(shù)。各個(gè)指數(shù)所代表的生態(tài)學(xué)意義及計(jì)算方法參見(jiàn)文獻(xiàn)[16]。景觀指數(shù)能夠定量地反映土地利用的結(jié)構(gòu)與空間配置,是景觀格局的高度濃縮。各景觀指數(shù)的定義及生態(tài)學(xué)意義如表1所示。

表1 景觀指數(shù)的方程及生態(tài)學(xué)意義
3 結(jié)果與分析
3.1 流域氣候變化特征
3.1.1氣溫與降水的變化趨勢(shì)額爾齊斯河流域多年平均溫度為3.26℃,從圖1A中可以看出,研究區(qū)近40 a的平均溫度整體呈上升趨勢(shì),增溫速率為0.37℃/10 a,且通過(guò)了Daniel顯著性檢驗(yàn)(p=0.0049),這與全球[1]、中亞[17]及新疆北部[18]氣溫不斷升高的變化趨勢(shì)一致。研究區(qū)氣溫在20世紀(jì)90年代中期前主要呈負(fù)距平,而之后呈正距平,特別是2000年以后升溫幅度明顯,是近40 a升溫幅度最大且持續(xù)時(shí)間最長(zhǎng)的一次升溫(圖1B)。

圖1 流域溫度序列及離差
研究區(qū)多年平均降水量為311 mm,屬于中亞干旱半干旱區(qū)的相對(duì)濕潤(rùn)區(qū)。從圖2A中可以看出,近40 a研究區(qū)降水的增加速率為9.77 mm/10 a,但未通過(guò)Daniel顯著性檢驗(yàn)。有研究表明,中亞地區(qū)各個(gè)國(guó)家的降水變化趨勢(shì)有所差異[19],本研究區(qū)降水的增加趨勢(shì)不明顯,但上游區(qū)降水量的顯著增加與中國(guó)西北干旱區(qū)[20]及新疆高山地區(qū)[21]的增濕趨勢(shì)一致。研究區(qū)降水在2000年以前呈正、負(fù)距平交替變化,總體平穩(wěn),其中最大低降水出現(xiàn)在1975年,2000年后降水以正距平為主,是近40 a降水增幅最大且持續(xù)時(shí)間最長(zhǎng)的一段時(shí)期(圖2B)。
綜上所述,研究區(qū)近40 年氣溫和降水均呈波動(dòng)上升變化趨勢(shì),水熱狀況趨于改善。溫度上升趨勢(shì)顯著(p=0.0049),降水增加趨勢(shì)不顯著(p=0.0663)。從時(shí)間尺度看,降水變化滯后于溫度4~5 年。
3.1.2氣溫與降水未來(lái)變化趨勢(shì)采用DFA法分析上游區(qū)和全流域氣溫、降水的未來(lái)變化趨勢(shì),對(duì)比子區(qū)域與整個(gè)流域未來(lái)氣候變化的差異。研究發(fā)現(xiàn),上游區(qū)氣溫和降水的DFA值分別為0.588,0.730,均介于0.5~1(圖3A—3B)。表明上游區(qū)未來(lái)氣溫和降水都將繼續(xù)保持增長(zhǎng)的趨勢(shì),而且溫度的持續(xù)上升相對(duì)于降水具有更高的可靠性。

圖2 流域降水序列及離差
對(duì)于流域整體,雖然氣溫和降水呈多年增長(zhǎng)趨勢(shì)(0.37℃/10 a,9.77 mm/10 a),但其DFA值分別為1.537,1.786,均大于1.5,屬布朗振蕩,無(wú)法判斷整個(gè)流域未來(lái)氣溫和降水變化的可能方向(圖3C—3D)。這種上游區(qū)未來(lái)增長(zhǎng)趨勢(shì)明確但全流域不明確的現(xiàn)象是由于上游區(qū)地形及地物結(jié)構(gòu)相對(duì)單一,而全流域環(huán)境復(fù)雜造成的。因此有必要分子系統(tǒng)對(duì)研究區(qū)氣候與生態(tài)的動(dòng)態(tài)變化作深入研究。

圖3 上游區(qū)域和整個(gè)流域DFA對(duì)比分析
3.2 MOD系統(tǒng)動(dòng)態(tài)變化分析
3.2.1MOD子系統(tǒng)劃分水熱條件對(duì)MOD各系統(tǒng)具有限制作用,綠洲與荒漠的相互演變過(guò)程就是以水熱為代表的環(huán)境條件逐步改善的過(guò)程[22]。在不考慮水文因子的情況下,對(duì)整個(gè)流域多年平均溫度和降水進(jìn)行全區(qū)域插值,分別劃分3級(jí),并將各區(qū)域值累積相乘后共構(gòu)成6級(jí),得到整個(gè)研究區(qū)的水熱分級(jí)圖,值越高的區(qū)域指示水熱狀況越好。在分級(jí)基礎(chǔ)上,再沿著干流從上游到下游選取3塊有代表性的子區(qū)域A(齋桑泊),B(塞米伊)和C(巴普洛達(dá)爾),分別代表山地系統(tǒng)、綠洲系統(tǒng)和荒漠系統(tǒng)(表2)。結(jié)合3個(gè)子區(qū)域近40 a內(nèi)的氣候和遙感影像數(shù)據(jù),分析研究MOD系統(tǒng)各組分變化特征及其與氣候變化和人類(lèi)活動(dòng)之間的關(guān)系。

表2 目標(biāo)子區(qū)域信息
3.2.2MOD子系統(tǒng)動(dòng)態(tài)變化分析A區(qū)作為典型的山地系統(tǒng),分析得出,A區(qū)年平均溫度為1.67 ℃,年降水量為292.5 mm,近40 a線性增溫和增濕率分別為0.39℃/10 a和4.8 mm/10 a。區(qū)域內(nèi),由于齋桑泊的存在,水體充足,達(dá)全區(qū)面積的15%;但裸地面積較大,達(dá)30%以上,廣泛分布于湖邊與山地;綠地(耕地和自然植被)主要分布于支流的階地及湖岸高地,面積約占全區(qū)面積的20%,自1979—1999年綠地面積有所增加,而1999年后有所減少;裸地面積由35.27%上升到47.19%且持續(xù)上升;水體面積穩(wěn)步小幅上升,從14.6%增加到了16.7% (圖4A)。
B區(qū)作為典型的綠洲系統(tǒng),年平均溫度為4.08℃,線性增溫率為0.16℃/10 a,年降水量為309.5 mm,無(wú)明顯升降趨勢(shì)。山前沖積平原水熱充足,該區(qū)域迎風(fēng)坡山地有大量的耕地和天然植被,植被覆蓋度高達(dá)40%~50%。其中,針葉林集中于額爾齊斯河北岸的西伯利亞低地,荒地和裸地主要分布在河流南部的哈薩克丘陵。從圖像分析也可知荒地和裸地占該區(qū)的20%左右。1979—1989年,水體面積由1979年的0.8%增加到了1989年的1.5%,綠地面積也增加了6.9%。1989—1999年沼澤和裸地面積明顯減少,分別減小了5.3%和3%,而荒地面積卻增加了7.7%。1999—2013年荒地和裸地?cái)U(kuò)大,達(dá)到了全區(qū)面積的46.7%,而綠地面積減小了14.3% (圖4B)。
C區(qū)作為典型的荒漠系統(tǒng),年平均溫度為2.87 ℃,線性增溫率為0.15℃/10 a,年降水量為326.8 mm,線性增濕率6.14 mm/10 a。該區(qū)是徑流耗散區(qū),水熱狀況一般,多年植被覆蓋度基本維持在50%左右。區(qū)域南部作為過(guò)渡緩沖帶的荒漠面積較大,達(dá)到了20%左右。區(qū)域內(nèi)有許多小型的水域,部分已干涸或鹽堿化。1979—1989年綠地面積基本維持不變,裸地以及鹽堿地面積增加,增幅分別為0.41%和3.46%。1989—1999年大量荒地退化成裸地,荒地隨之減少了18.8%,其他地物面積幾乎不變;1999—2013年綠地面積擴(kuò)大了10.99%,荒地面積,沼澤和水體都有小幅的減少(圖4C)。
3.2.3景觀格局動(dòng)態(tài)變化分析A,B,C這3個(gè)子區(qū)域中,B和C區(qū)受人類(lèi)活動(dòng)的影響較大,因此其景觀格局的變化較大。利用Fragstats軟件,從景觀類(lèi)型與斑塊水平兩個(gè)角度計(jì)算景觀指數(shù)。


圖4 1979-2013年A,B和C區(qū)土地利用類(lèi)型變化
從表3中可知,B區(qū)和C區(qū)的斑塊數(shù)量從1979—1999年持續(xù)增加,1999—2013年略有減少,可見(jiàn)區(qū)域破碎性增加。B區(qū)景觀Shannon多樣性指數(shù)和Shannon均勻度指數(shù)先降低再升高,說(shuō)明該區(qū)域的景觀異質(zhì)性程度呈現(xiàn)先下降再上升的趨勢(shì),景觀類(lèi)型趨于多樣化發(fā)展。C區(qū)景觀Shannon多樣性指數(shù)和Shannon均勻度指數(shù)則先上升再下降,說(shuō)明該區(qū)域景觀異質(zhì)性程度呈現(xiàn)先上升再下降的趨勢(shì),景觀類(lèi)型越豐富,景觀所含的信息量和信息的不定性越大。從年份對(duì)比可以看出1999年B區(qū)和C區(qū)景觀形狀結(jié)構(gòu)最為復(fù)雜。而1999年以后,斑塊形狀指數(shù)減小,斑塊復(fù)雜性降低,說(shuō)明區(qū)域的景觀結(jié)構(gòu)與形狀區(qū)域規(guī)則化、聚集化。
3.2.4MOD系統(tǒng)對(duì)氣候變化和人類(lèi)活動(dòng)的響應(yīng)由景觀變化分析可知1999年前后各子系統(tǒng)都發(fā)生了較其他時(shí)期更為明顯的變化。區(qū)分MOD各子系統(tǒng),進(jìn)一步討論其受氣候變化和人類(lèi)活動(dòng)的影響。
(1) 山地系統(tǒng)。由流域氣候的暖濕化趨勢(shì)以及A區(qū)近年水體和沼澤面積的變化(圖4和表3)可以推斷,地表水資源得到了暫時(shí)性增加。主要是因?yàn)闅鉁厣邔?dǎo)致的高山冰雪融水增加[17]。A區(qū)氣溫和降水的增長(zhǎng)率分別為0.39℃/10 a和4.8 mm/10 a,而綠地面積與平均溫度和降水的相關(guān)系數(shù)分別為0.53和0.98,氣溫和降水的增加使得A區(qū)綠地面積在1979—1999年增加,但2013年有所減少。人類(lèi)活動(dòng)使得該區(qū)自然植被面積有所減少,而耕地面積增加;不可用地面積(裸地和鹽堿—沙地)與平均溫度和降水的相關(guān)系數(shù)分別為0.64和-0.84,不可用地面積逐年增加。由進(jìn)一步的遙感解譯可知,區(qū)域內(nèi)自然植被主要是山地林帶,1989年林地面積的增加,一定程度上抑制了不可用地的增加,1999—2013年,隨著林地面積大幅減少,不可用地的面積增加迅速(表4)。

表4 A區(qū)域地物景觀百分比 %
(2) 綠洲系統(tǒng)。1979—1989年,隨著B(niǎo)區(qū)水庫(kù)的建立,水體面積顯著增加,至1999年后維持相對(duì)穩(wěn)定的規(guī)模;沼澤面積有一定的起伏,總體呈下降趨勢(shì);綠地面積同平均溫度和降水的相關(guān)性都不高,分別為0.48和-0.36。伴隨著區(qū)域水體面積的顯著增加,引水灌溉和開(kāi)墾力度加大,耕地面積隨之增加。而自然植被受限于氣候和水文條件,1999—2013年退縮明顯。荒地面積自1989年后持續(xù)增加。不可用地面積與平均溫度的相關(guān)性較好,相關(guān)系數(shù)為-0.80,與降水的相關(guān)性較差,相關(guān)系數(shù)為-0.1。1979—1999年受綠地面積擴(kuò)大和氣溫上升的影響,不可用地面積逐漸縮減。相應(yīng)的,1999—2013年綠地規(guī)模退縮,不可用地面積擴(kuò)大(表5)。

表5 B區(qū)域地位景觀百分比 %
(3) 荒漠系統(tǒng)。1979—1989年C區(qū)水體面積有所增加,之后小幅減少。綠地主要分布在干流附近及河流北岸,同平均溫度和降水的相關(guān)系數(shù)分別為0.65和0.95,與降水的相關(guān)性較大。1979—1989年植被面積變化不大,主要是其他地向耕地的轉(zhuǎn)化;1989—2013年隨著人類(lèi)活動(dòng)影響,水熱狀況的改善,耕地和自然植被面積擴(kuò)大。不可用地面積與平均溫度和降水的相關(guān)系數(shù)分別為0.88和0.61,與溫度的相關(guān)性較大。區(qū)域內(nèi)溫度升高,且伴隨著綠地向外擴(kuò)張,導(dǎo)致1989年以來(lái)荒地面積減少,但不可用地面積增加(表6)。

表6 C區(qū)域地物景觀百分比 %
4 討 論
山地高海拔冰川和積雪融水是額爾齊斯河的主要補(bǔ)給源,氣溫升高使得冰雪融水增多,伴隨流域內(nèi)降水量的增加,使得流域內(nèi)植被覆蓋面積增加[23]。但流域內(nèi)由于人為開(kāi)發(fā),生態(tài)環(huán)境進(jìn)一步變得復(fù)雜,對(duì)整個(gè)生態(tài)系統(tǒng)的影響也更加明顯。山地系統(tǒng)中,非農(nóng)業(yè)用地主要為山區(qū)森林,山區(qū)森林面積的增加有助于抑制不可用地面積的增加??梢?jiàn),森林對(duì)穩(wěn)定系統(tǒng)組分和調(diào)節(jié)氣候和水資源供給是非常重要的。
綠洲作為整個(gè)MOD系統(tǒng)中最為活躍的組分,水熱狀況較為良好,同時(shí)是人類(lèi)活動(dòng)的中心,受到自然和人為因素的共同影響與制約。綠洲內(nèi)溫度、降水和植被覆蓋變化導(dǎo)致蒸散發(fā)率、反照率、地表粗糙率等綜合效應(yīng)的不同,也會(huì)對(duì)自身起到調(diào)節(jié)作用[24]。同時(shí)人類(lèi)生產(chǎn)生活方式的變化導(dǎo)致綠洲內(nèi)景觀變化更為劇烈,生態(tài)環(huán)境的脆弱性程度加大[25]。
荒漠作為綠洲系統(tǒng)外圍的屏障和儲(chǔ)備[26],緩沖了干旱半干旱區(qū)惡劣氣候?qū)τ诰G洲的侵害,為綠洲的擴(kuò)展和演替提供了基礎(chǔ)。除氣候變化對(duì)其造成的長(zhǎng)時(shí)間尺度的影響外,目前荒漠系統(tǒng)受人類(lèi)的干預(yù)也變得越來(lái)越多,生態(tài)環(huán)境變得十分敏感。人類(lèi)活動(dòng)的持續(xù)增強(qiáng),植被覆蓋面積將進(jìn)一步減少而裸地面積增加,最終導(dǎo)致更大面積的鹽堿化和沙化,使得區(qū)域的脆弱性增加[27],并減弱未來(lái)農(nóng)業(yè)生產(chǎn)和生態(tài)的適應(yīng)性。例如20世紀(jì)著名的咸海危機(jī)[28]和塔里木河生態(tài)水文危機(jī)[29]被認(rèn)為是由于過(guò)多的人類(lèi)活動(dòng)造成的。
本文僅對(duì)MOD系統(tǒng)的土地類(lèi)型進(jìn)行了基本的分類(lèi)分析,但人工綠地、森林、灌木和草地沒(méi)有被進(jìn)一步分類(lèi)。因此,人類(lèi)和自然因素對(duì)MOD系統(tǒng)的影響需要進(jìn)一步更精確的分類(lèi)和定量分析。
5 結(jié) 論
(1) 研究區(qū)整體變得暖濕,水熱狀況趨于良性發(fā)展,但溫度和降水并不協(xié)同變化,加上地形地貌的差異,各個(gè)子系統(tǒng)對(duì)氣候變化和人類(lèi)活動(dòng)的響應(yīng)存在差異。1999年是研究區(qū)景觀格局的轉(zhuǎn)折點(diǎn),景觀結(jié)構(gòu)、形狀和數(shù)量都在一定程度上有所增加,這一轉(zhuǎn)變與人類(lèi)活動(dòng)的影響密不可分。
(2) A區(qū)山地系統(tǒng)也受到人類(lèi)活動(dòng)的影響,但影響較小。其年均氣溫最小,線性升溫率卻最高。近40 a來(lái),在氣候變暖的影響下,水體面積小幅增長(zhǎng),植被面積先增加后減少,主要表現(xiàn)為耕地面積擠占自然植被面積而使得自然植被面積的減少。裸地面積在1999年前幾乎無(wú)變化而之后有所增加。研究區(qū)土地利用類(lèi)型向不可利用地類(lèi)轉(zhuǎn)換的面積越來(lái)越大。
(3) B區(qū)綠洲系統(tǒng)地處山前沖擊平原,受氣候和人類(lèi)的共同影響,近40 a來(lái)該區(qū)氣溫有所升高,但降水量無(wú)明顯變化。該區(qū)植被覆蓋度高,但天然植被面積自1999年后降低,而耕地面積大幅增加,可知人類(lèi)耕作活動(dòng)對(duì)這一地區(qū)土地利用改變十分明顯。
(4) C區(qū)增濕率在3個(gè)區(qū)域最高,多年植被覆蓋度基本維持在50%左右,但結(jié)構(gòu)單一,區(qū)內(nèi)鹽堿化較嚴(yán)重,受綠洲擴(kuò)張的影響,荒地面積自1979年以來(lái)顯著減小。
參考文獻(xiàn):
[1]Stocker T F, Qin D, Plattner G K, et al. The physical science basis: contribution of working group to the fifth assessment report of the intergovernmental panel on climate change[M]∥IPCC. IPCC 2013. Cambridge, UK:Cambridge University Press, 2013.
[2]Evans D D, Thames J L. Water in desert ecosystems[J]. Journal of Range Management, 1981,35(4):543.
[3]Baird A, Wilby R. Eco-hydrology:plants and water in terrestrial and aquatic environments[J]. Journal of Ecology, 1999,88(6):1095-1096.
[4]Hrkal Z, Gadalia A, Rigaudiere P. Will the river Irtysh survive the year 2030? Impact of long-term unsuitable land use and water management of the upper stretch of the river catchment(North Kazakhstan)[J]. Environmental Geology, 2006,50(5):717-723.
[5]龍愛(ài)華,鄧銘江,李湘權(quán),等.哈薩克斯坦水資源及其開(kāi)發(fā)利用[J].地球科學(xué)進(jìn)展,2010,25(12):1357-1366.
[6]李捷,夏自強(qiáng),郭利丹,等.額爾齊斯河流域氣候特征及變化趨勢(shì)分析[J].河海大學(xué)學(xué)報(bào):自然科學(xué)版,2008,36(3):311-315.
[7]雷雨,龍愛(ài)華,鄧銘江,等.1926—2009年額爾齊斯河流域中游地區(qū)氣候變化及其對(duì)水資源的影響分析[J].冰川凍土,2012,34(4):912-919.
[8]Beurs D, Kirsten M, Henebry G M. Land surface phenology, climatic variation, and institutional change:analyzing agricultural land cover change in Kazakhstan[J]. Remote Sensing of Environment, 2004,89(4):497-509.
[9]沈金祥,陳曦,楊遼,等.額爾齊斯河—齋桑湖流域近20年來(lái)土地利用/土地覆被時(shí)空演變[J].干旱區(qū)地理,2010,33(2):189-195.
[10]Chuluun T, Ojima D. Land use change and carbon cycle in arid and semi-arid lands of East and Central Asia[J]. Science in China Series C:Life Sciences-English Edition, 2002,45(S1):48-54.
[11]王讓會(huì),張慧芝,趙振勇,等. MODS耦合關(guān)系的界面過(guò)程及其特征[J].干旱區(qū)地理,2004,27(3):388-392.
[12]Jing X H. Landscape pattern and fragmentation in Beitun of Irtysh River Basin, Xinjiang[J]. Scientia Silvae Sinicae, 2008,44(3):21-28.
[13]Wilks D S. Statistical methods in the atmospheric sciences[M]. Oxford: Academic Press, 2011.
[14]Peng C K, Havlin S, Stanley H E, et al. Quantification of scaling exponents and crossover phenomena in nonstationary heartbeat time series[J]. Chaos An Interdisciplinary Journal of Nonlinear Science, 1995,5(1):82-87.
[15]高沈瞳,徐長(zhǎng)春.額爾齊斯河流域上游地區(qū)近50 a來(lái)氣溫和降水變化的DFA分析[J].冰川凍土,2014,36(3):706-716.
[16]余新曉.景觀生態(tài)學(xué)[M].北京:高等教育出版利,2006.
[17]Chen F, Wang J, Jin L, et al. Rapid warming in mid-latitude central Asia for the past 100 years[J]. Frontiers of Earth Science, 2009,3(1):42-50.
[18]馬柱國(guó),符淙斌,任小波,等.中國(guó)北方年極端溫度的變化趨勢(shì)與區(qū)域增暖的聯(lián)系[J].地理學(xué)報(bào),2003,58(9):11-19.
[19]Stulina G, Eshchanov O. Climate change impacts on hydrology and environment in the Pre-Aral region[J]. Quaternary International, 2013,311(454):87-96.
[20]Wang S J, Zhang M J, Sun M P, et al. Changes in precipitation extremes in alpine areas of the Chinese Tianshan Mountains, central Asia, 1961—2011[J]. Quaternary International, 2013,311:97-107.
[21]Wang H J, Chen Y N, Chen Z S. Spatial distribution and temporal trends of mean precipitation and extremes in the arid region, northwest of China, during 1960—2010[J]. Hydrological Processes, 2013,27(12):1807-1818.
[22]潘曉玲,馬映軍,顧峰雪.中國(guó)西部干旱區(qū)生態(tài)環(huán)境演變與調(diào)控研究進(jìn)展與展望[J].地球科學(xué)進(jìn)展,2003,18(1):50-57.
[23]沈永平,蘇宏超,王國(guó)亞,等.新疆冰川、積雪對(duì)氣候變化的響應(yīng)(I):水文效應(yīng)[J].冰川凍土,2013,35(3):513-527.
[24]Chen Y N, Xu C C, Chen Y P, et al. Response of glacial-lake outburst floods to climate change in the Yarkant River basin on northern slope of Karakoram Mountains, China[J]. Quaternary International, 2010,226(1/2):75-81.
[25]Forster P, Ramaswamy V, Artaxo P, et al. Changes in atmospheric constituents and in radiative forcing[M]. Cambridge, Uk: Cambridge University Press, 2007.
[26]Su Y Z, Zhao W Z, Su P X, et al. Ecological effects of desertification control and desertified land reclamation in an oasis-desert ecotone in an arid region:a case study in Hexi Corridor, northwest China[J]. Ecological Engineering, 2007,29(2):117-124.
[27]Lioubimtseva E, Henebry G M. Climate and environmental change in arid Central Asia: impacts, vulnerability, and adaptations[J]. Journal of Arid Environments, 2009,73(11):963-977.
[28]Mainguet M, Létolle R. Human-made Desertification in the Aral Sea Basin:Planning and Management Failures[J].1998,41(3):129-142.
[29]Chen Y N, Pang Z H, Chen Y P, et al. Response of riparian vegetation to water-table changes in the lower reaches of Tarim River, Xinjiang Uygur, China[J]. Hydrogeology Journal, 2008,16(7):1371-1379.