黑河流域中游地區(qū)生態(tài)環(huán)境變化特征及驅(qū)動力

周 沙,黃躍飛,王光謙

(清華大學(xué)水沙科學(xué)與水利水電工程國家重點實驗室,北京 100084)

黑河流域中游地區(qū)生態(tài)環(huán)境變化特征及驅(qū)動力

周 沙,黃躍飛*,王光謙

(清華大學(xué)水沙科學(xué)與水利水電工程國家重點實驗室,北京 100084)

多角度分析黑河流域中游地區(qū)植被覆蓋度時空變化特征,并建立基于演變過程的生態(tài)系統(tǒng)灰色關(guān)聯(lián)度模型,分析生態(tài)環(huán)境變化的驅(qū)動因子.研究表明：1)1999～2008年,平原旱地、低覆蓋草地、有林地年最大歸一化差異植被指數(shù)(NDVI)增幅較大,達 0.1～0.2,疏林地2004年后保持高速增長,年最大NDVI增幅0.208,增長了77.6%; 2)張臨高盆地年最大植被覆蓋度線性擬合年均增幅0.0063,生長季節(jié)平均植被覆蓋度在小幅波動中呈現(xiàn)穩(wěn)定增長趨勢,擬合優(yōu)度達 0.74;甘州區(qū)生態(tài)治理成效最顯著,年增長幅度集中在0.03～0.3;臨澤縣和高臺縣以增長為主,但局部過渡帶地區(qū)下降幅度達 0.1～0.3; 3)植被 7～10月覆蓋度呈現(xiàn)明顯增長趨勢,峰值從 6～7月延遲到7～8月,2007年達0.39;植被覆蓋度分級結(jié)構(gòu)呈現(xiàn)優(yōu)化趨勢,極低覆蓋度植被逐漸轉(zhuǎn)化為低覆蓋度植被,2007年相比2000年降低25%以上,高覆蓋度植被1998～2008年間增長約16%; 4)根據(jù)3種植被覆蓋度變化與各驅(qū)動因子關(guān)聯(lián)分析,氣象水文因子主要包括降雨量、蒸發(fā)量、徑流量,最大關(guān)聯(lián)度分別為 0.91、-0.83、0.76,社會經(jīng)濟因子主要包括農(nóng)作物播種面積、第一產(chǎn)業(yè)產(chǎn)值、農(nóng)業(yè)科技水平,最大關(guān)聯(lián)度分別為-0.81、0.78、0.81.

黑河流域；植被覆蓋度；時空變化特征；灰色關(guān)聯(lián)度

20世紀 60年代以來,隨著人口的增長和經(jīng)濟快速發(fā)展,黑河流域水土資源開發(fā)過度,中游地區(qū)水資源短缺、水環(huán)境惡化、土壤鹽堿化嚴重.為了遏制黑河流域生態(tài)環(huán)境不斷惡化的趨勢,21世紀初,國家實施了黑河流域近期治理工程.2001年 8月,甘肅省黑河流域近期治理工程項目陸續(xù)實施,加強黑河干流水資源統(tǒng)一管理,在中游地區(qū)提高用水效率,建設(shè)節(jié)水型社會,改善生態(tài)環(huán)境.2010年10月,項目通過水利部總體竣工驗收,實現(xiàn)了國務(wù)院批準的《黑河流域近期治理規(guī)劃》目標(biāo).

灰色關(guān)聯(lián)分析是根據(jù)比較因子與參考因子之間發(fā)展趨勢的相似或相異程度,來判斷因子間關(guān)聯(lián)程度的一種方法[1].目前測算灰色關(guān)聯(lián)度的模型很多,如鄧氏關(guān)聯(lián)度[2]、斜率關(guān)聯(lián)度[3]、廣義灰色關(guān)聯(lián)度[4]、灰色綜合關(guān)聯(lián)度[5]等,但這些模型都存在一個共同的問題,不能反映正、負相關(guān)關(guān)系,且存在正負積分相互抵消而出現(xiàn)明顯不合理的結(jié)果.近年來,也有許多學(xué)者提出改進的灰色關(guān)聯(lián)度模型[6-8],但都存在一定的局限性,難以統(tǒng)一評價氣候水文、社會經(jīng)濟等多方面因素對于生態(tài)環(huán)境變化的影響大小.

本文以張臨高盆地為研究區(qū),分析黑河流域中游地區(qū)近期治理項目對改善生態(tài)環(huán)境的作用效果,分別從不同植被覆蓋類型區(qū)NDVI變化特征、植被覆蓋度逐年變化特征、空間分布特征、時空統(tǒng)計特征4方面,多角度分析黑河流域中游地區(qū)生態(tài)環(huán)境的變化特征.基于研究區(qū)的實際情況,在曹明霞[9]的基礎(chǔ)上,提出了基于演變過程的生態(tài)系統(tǒng)灰色關(guān)聯(lián)度模型,適用于研究干旱區(qū)生態(tài)環(huán)境變化與各影響因子間的關(guān)聯(lián)度,分析中游地區(qū)生態(tài)環(huán)境變化的驅(qū)動力.

1 研究區(qū)概況

黑河流域中游張臨高盆地土地資源豐富,是我國西北地區(qū)重要的商品糧生產(chǎn)基地,也是黑河徑流的主要利用區(qū).本文研究區(qū)包括張掖地區(qū)的甘州區(qū)、臨澤縣和高臺縣(圖1).研究區(qū)屬溫帶大陸性干旱氣候,降雨少而集中于夏季,年均降水量50～250mm,日照充足,蒸發(fā)強烈,年均蒸發(fā)量2000～3500mm,晝夜溫差很大[10].20世紀90年代以來,張臨高盆地人口逐年增加,從1990年71.05萬人增加到2010年82.75萬人,社會經(jīng)濟發(fā)展十分迅速,地區(qū)生產(chǎn)總值從1990年8.78億元增加到2010年147.01億元[11].

圖1 研究區(qū)示意Fig.1 The location map of the study area

2 數(shù)據(jù)與方法

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

本次研究中,行政邊界、氣候、水文、植被等數(shù)據(jù)來源于國家自然科學(xué)基金委員會“中國西部環(huán)境與生態(tài)科學(xué)數(shù)據(jù)中心” (http∶//westdc. westgis.ac.cn),包括數(shù)據(jù)集有“黑河流域行政邊界數(shù)據(jù)集”、“黑河流域逐日水文觀測數(shù)據(jù)集”、“黑河流域逐日氣象觀測數(shù)據(jù)集”[12]、“黑河流域長時間序列SPOT_Vegetation植被數(shù)據(jù)集(NDVI)”[13].其中,歸一化差異植被指數(shù)(NDVI)時間序列為1998年4月～2008年7月.社會經(jīng)濟數(shù)據(jù)來源于《甘肅省統(tǒng)計年鑒》(1998-2008年)[14].

2.2 研究方法

基于NDVI時間序列數(shù)據(jù),估算有植被覆蓋區(qū)域的植被覆蓋度,計算植被覆蓋度變化趨勢,研究生態(tài)環(huán)境變化特征,并建立生態(tài)系統(tǒng)灰色關(guān)聯(lián)度模型,分析生態(tài)環(huán)境變化的驅(qū)動力.

2.2.1 植被覆蓋度計算 植被覆蓋度是指植被(包括葉、莖、枝)在地面的垂直投影面積占統(tǒng)計區(qū)總面積的百分比,常用于植被變化、生態(tài)環(huán)境研究[15].NDVI 是植被生長狀態(tài)及植被覆蓋度的最佳指示因子,其范圍在-1～1之間變動[16].選擇研究區(qū)內(nèi)多年平均年最大NDVI大于0.1的象元,作為有植被覆蓋的區(qū)域[17].

NDVI與植被覆蓋分布呈線性相關(guān),根據(jù)象元二分模型原理,每個象元的NDVI值可以表示為有植被覆蓋部分的 NDVI與無植被覆蓋部分的NDVI線性組合的形式[18].因此,計算植被覆蓋度(VFC)的公式可表示為∶

式中,N DVIsoil為無植被覆蓋象元的 NDVI,理論上接近 0,但受眾多因素影響,變化范圍一般在-0.1～0.2[19];N DVIveg表示純植被象元的NDVI,隨著植被類型和植被的時空分布而變化.本文采用李苗苗[20]的研究方法計算兩者的值.

2.2.2 植被覆蓋度變化特征計算 計算研究區(qū)逐個像元1999～2007年各年的年平均、年最大、生長季節(jié)平均的VFC,采用一元線性回歸法計算每個有植被覆蓋網(wǎng)格9年年平均、年最大、生長季節(jié)平均植被覆蓋度的總體變化趨勢,用SLOPE線性趨勢線估計其變化幅度,計算公式為[21]∶

式中,VFC可代表1999～2007逐年最大、年平均、年生長季節(jié)平均植被覆蓋度,SLOPE表示1999～2007年年平均、年最大、生長季節(jié)平均植被覆蓋度的增長值,分別得到總增長幅度及增長百分率為[22]∶

2.2.3 基于演變過程的生態(tài)系統(tǒng)灰色關(guān)聯(lián)度模型 由于干旱區(qū)生態(tài)環(huán)境變化受到自然因素和人類活動的共同影響,自然因素如氣候變化的波動范圍較小,而社會經(jīng)濟發(fā)展一般呈現(xiàn)穩(wěn)定增長,變化幅度較大.為了協(xié)調(diào)各驅(qū)動因子變化趨勢對生態(tài)環(huán)境變化的影響,本文基于生態(tài)環(huán)境因子與驅(qū)動因子在研究期內(nèi)變化的相對斜率,建立基于演變過程的生態(tài)系統(tǒng)灰色關(guān)聯(lián)度模型.模型中定義了完全正相關(guān)(灰色關(guān)聯(lián)度為1)和完全負相關(guān)(灰色關(guān)聯(lián)度為-1)情況的內(nèi)涵,即參考序列和比較序列存在正、負相關(guān)關(guān)系(擬合優(yōu)度為1),能很好地反映因子間的正負相關(guān)關(guān)系.模型建立各驅(qū)動因子與生態(tài)環(huán)境變化之間的相關(guān)程度,從而可以進行生態(tài)環(huán)境變化的驅(qū)動力分析.模型如下∶

說明∶序列X0和序列Xi的灰色關(guān)聯(lián)度滿足∶當(dāng)且僅當(dāng)Xi(t)=為常數(shù)),

3 結(jié)果與討論

3.1 不同植被覆蓋類型區(qū)NDVI變化特征

圖2中,除戈壁區(qū),各類型區(qū)年最大NDVI均呈現(xiàn)明顯的增長趨勢,且峰谷的波動與鶯落峽徑流量存在一定的對應(yīng)關(guān)系,如2003、2007年,多數(shù)類型區(qū)NDVI為波峰年.相比1999年(波谷),平原旱地、低覆蓋草地、有林地均有大幅增長,增長幅度達0.1～0.2.疏林地NDVI在2004年后保持高速增長,4年間增長幅度0.208,增長了77.6%.總體來看,2004年以前,各類型區(qū) NDVI波動較大,2004年后,在小幅波動中呈現(xiàn)穩(wěn)定增長趨勢,這與黑河干流水量統(tǒng)一調(diào)度管理是分不開的.

圖2 各植被類型區(qū)年最大NDVI逐年變化Fig.2 The yearly changes of annual maximum NDVI in different vegetation types

分析不同植被覆蓋類型區(qū)多年平均 NDVI逐月變化(圖 3),可以看出,張臨高盆地植被的主要生長季節(jié)為5～10月,各類型區(qū)生長期長度略有差異,但變化趨勢與鶯落峽逐月徑流量基本一致.有林地NDVI峰值最高,其次是高覆蓋草地、灌木林地、中覆蓋草地、平原旱地、低覆蓋度草地、疏林地,戈壁和鹽堿地NDVI全年較低,幾乎無植被覆蓋.峰值代表植被全年生長最好時的情況,在一定程度上說明植被覆蓋狀況的好壞,也是生長開始衰退的時間,除戈壁和鹽堿地,其余各類型區(qū)植被到達峰值的時間均為7月.其中,灌木林地非生長季節(jié)(11～4月)NDVI最高,0.2左右,幾乎是其它植被類型的2～3倍,可以看出灌木林對改善生態(tài)環(huán)境的意義重大.平原旱地全年NDVI處于中等水平,且在7月收獲季節(jié)一過,NDVI迅速降低,對生態(tài)環(huán)境不利.

圖3 各植被類型區(qū)年最大NDVI逐月變化Fig.3 The monthly changes of annual maximum NDVI in different vegetation types

3.2 植被覆蓋度逐年變化特征

圖4 植被覆蓋度逐年變化Fig.4 The yearly changes of vegetation coverage

分析研究區(qū)內(nèi)年最大、生長季節(jié)平均、年平均植被覆蓋度逐年變化情況(圖 4).其中,年最大植被覆蓋度增長速度最快,線性擬合年均增幅0.0063,波動也最大;生長季節(jié)平均植被覆蓋度在小幅波動中呈現(xiàn)穩(wěn)定增長趨勢,擬合優(yōu)度達0.74.年平均植被覆蓋度略有增長,波動幅度較小.總體來看,經(jīng)過一系列生態(tài)治理,張臨高盆地植被覆蓋度增長趨勢明顯,生態(tài)環(huán)境逐漸得到改善.

3.3 植被覆蓋度空間分布特征

采用一元線性回歸分析方法計算研究區(qū)各象元1999～2007年各年年平均、年最大、生長季節(jié)平均植被覆蓋度的變化趨勢線,根據(jù)趨勢線斜率計算植被覆蓋度增長幅度和增長百分率,并繪制空間分布圖.其中,生長季節(jié)平均植被覆蓋度變化幅度和變化百分率空間分布如圖4所示,未計算區(qū)域為多年平均年最大NDVI值小于0.1的區(qū)域,作為荒漠區(qū)處理.

由圖5可見,生態(tài)治理成效最顯著的是張掖市甘州區(qū),植被覆蓋度增加明顯,年增長幅度集中在 0.03～0.3,其后依次是臨澤縣和高臺縣.從植被覆蓋度變化情況來看,高臺縣和臨澤縣北部被荒漠包圍的過渡帶地區(qū),以及沿黑河分布的部分耕地地區(qū),局部年下降幅度達 0.1～0.3,年降低百分率為 20%～60%.隨著過渡帶地區(qū)由于可利用的徑流性水資源減少,對降雨性水資源依賴程度大,使得植被生長的峰值情況受到降雨量的限制.

圖5 張臨高盆地1999～2007生長季節(jié)平均植被覆蓋度變化幅度及變化百分率Fig.5 The changes in amplitude and percentage changes of growth season average vegetation cover in the study area (1999～2007)

3.4 植被覆蓋度時空統(tǒng)計特征

比較 3種植被覆蓋度的變化度(圖 6),年最大、生長季節(jié)平均植被覆蓋度的好轉(zhuǎn)趨勢較明顯,分別有70%、65%以上分布在0.03～0.3,而年平均植被覆蓋度變化 80%以上分布在 0～0.1,因為植被覆蓋度的變化主要體現(xiàn)在生長季節(jié).從變化百分率來看,年最大植被覆蓋度和生長季節(jié)平均植被覆蓋度的變化分布基本一致,兩者的增長百分率分布在20%～60%的比例占50%以上,而年平均植被覆蓋度的好轉(zhuǎn)趨勢相對較差些.

從生長季節(jié)植被覆蓋度變化規(guī)律來看(圖7),各年植被覆蓋度到達峰值的時間并不一致,2002年之前,植被覆蓋度偏低,峰值到來時間較早,6～7月植被覆蓋度較高,8月植被覆蓋度已處于明顯下降階段.2002年及之后,植被覆蓋度較高,峰值到達時間較晚,7～8月植被覆蓋度較高,在0.35左右,2003年和2006年峰值出現(xiàn)在8月,2007年峰值達到 0.3924.隨著張臨高盆地灌溉條件的改善,7～10月植被覆蓋度呈現(xiàn)明顯增長趨勢,同時生長季節(jié)時間長度有增加趨勢,而植被覆蓋度峰值提高、生長期延長,是區(qū)域生態(tài)環(huán)境改善的標(biāo)志.

圖6 植被覆蓋度變化幅度、變化百分率分級百分比統(tǒng)計Fig.6 The statistical results of vegetation cover change classification

利用年最大植被覆蓋度(植被覆蓋度峰值數(shù)據(jù),采用6級表示法(0～15%、15%～30%、30%～45%、45%～60%、60～75%、75%～100%)進行植被覆蓋度分級統(tǒng)計,研究各級覆蓋度植被面積百分比逐年變化情況.根據(jù)分級情況,將植被覆蓋度分為極低、低、中低、中、高中、高植被覆蓋度區(qū).由圖8可以看出,植被覆蓋度變化主要發(fā)生在低覆蓋度和極低覆蓋度區(qū),兩者波動幅度較大,互為增減關(guān)系,累積面積比占50%以上.高覆蓋度植被呈現(xiàn)明顯增長趨勢,10年間增長了約 16%.中高、中、中低覆蓋度比例在波動中略有下降,三者比例基本相同,累積比例約 34%.整體而言,植被覆蓋度結(jié)構(gòu)呈現(xiàn)優(yōu)化趨勢,主要體現(xiàn)在高覆蓋度植被比例的增長和極低覆蓋度植被逐漸轉(zhuǎn)化為低覆蓋度植被,可見張掖市節(jié)水灌溉、退耕還林還草等一系列工程措施取得了不錯的成效.

圖7 生長季節(jié)植被覆蓋度逐年變化Fig.7 The changes of the vegetation cover from May to October(1998～2007)

圖8 年最大植被覆蓋度分級統(tǒng)計逐年變化Fig.8 The statistical results of yearly maximum vegetation cover classification (1998～2008)

3.5 關(guān)聯(lián)度計算

應(yīng)用基于演變過程的生態(tài)系統(tǒng)灰色關(guān)聯(lián)度模型,計算 3種植被覆蓋度與各氣候水文、社會經(jīng)濟因子之間的關(guān)聯(lián)度,結(jié)果如表1所示.

從年最大、年平均、生長季節(jié)平均植被覆蓋度與各影響因子關(guān)聯(lián)度分析可知,3種植被覆蓋度與各影響因子關(guān)聯(lián)度基本一致.

針對關(guān)聯(lián)度 0.5以上的影響因子,根據(jù)生長季節(jié)平均植被覆蓋度進行排序,正相關(guān)影響因子∶降水量>農(nóng)村用電量>化肥施用折純量>第一產(chǎn)業(yè)產(chǎn)值>人均生產(chǎn)總值>鶯落峽徑流量>農(nóng)業(yè)人口;負相關(guān)影響因素∶農(nóng)作物播種面積>蒸發(fā)量>第二產(chǎn)業(yè)產(chǎn)值.

表1 年最大、年平均、生長季節(jié)平均植被覆蓋度與各影響因子關(guān)聯(lián)度計算結(jié)果Table 1 The calculation results of correlation between vegetation coverage and driving factors

圖9 年最大、年平均、生長季節(jié)平均植被覆蓋度與各影響因子關(guān)聯(lián)度分析Fig.9 The calculation results of correlation between vegetation coverage and driving factors

分析各水文氣象因素對植被覆蓋度的影響,蒸發(fā)量和降雨量影響較大,其中,蒸發(fā)量與年平均植被覆蓋度呈顯著負相關(guān),關(guān)聯(lián)度達-0.8264;降雨量與生長季節(jié)植被覆蓋度顯著正相關(guān),關(guān)聯(lián)度達0.9071,這與干旱區(qū)的氣候特點密不可分,因為植被覆蓋度變化主要發(fā)生在以降雨性水資源補給為主的極低、低植被覆蓋區(qū),且植被覆蓋度變化集中在生長季節(jié).實施水量統(tǒng)一調(diào)度工程后,徑流量的利用效率提高,鶯落峽來水量與植被覆蓋度的關(guān)聯(lián)性也較強.

分析各社會經(jīng)濟因素對植被覆蓋度的影響,第一產(chǎn)業(yè)總值、人均生產(chǎn)總值與年最大植被覆蓋度呈現(xiàn)顯著正相關(guān),均為0.75以上,而農(nóng)作物播種面積與生長季節(jié)植被覆蓋度負相關(guān)達 0.8以上.2000年之后,張掖市節(jié)水型社會建設(shè),不斷優(yōu)化水資源利用效率,開展退耕還林還草工程,農(nóng)作物播種面積下降,林牧業(yè)得到發(fā)展,高覆蓋度植被面積增加約 16%,在第一產(chǎn)業(yè)進步的同時促進了生態(tài)環(huán)境的整體改善.代表農(nóng)業(yè)科技進步的農(nóng)村用電量和化肥施用折純量與生長季節(jié)植被覆蓋度呈顯著正相關(guān)關(guān)系,關(guān)聯(lián)度分別達 0.8121、0.7948,這進一步證明了農(nóng)業(yè)技術(shù)進步有利于耕地乃至區(qū)域植被覆蓋度的提高.而第二產(chǎn)業(yè)產(chǎn)值與植被覆蓋度呈現(xiàn)負相關(guān),說明工礦、建筑業(yè)等的發(fā)展對區(qū)域植被覆蓋度有一定的負面影響.

基于生態(tài)環(huán)境變化特征和驅(qū)動力分析,中游地區(qū)生態(tài)環(huán)境在整體改善的基礎(chǔ)上,仍存在不可忽視的問題,比如北部荒漠包圍的過渡帶地區(qū)、中部沿黑河分布的條帶狀耕地覆蓋度呈降低趨勢.因此,有必要通過水資源合理配置和調(diào)度管理,在不引起耕地區(qū)植被覆蓋度降低和農(nóng)業(yè)產(chǎn)量減少的基礎(chǔ)上,保證過渡帶植被生態(tài)用水,進一步改善中游地區(qū)生態(tài)環(huán)境.

4 結(jié)論

4.1 不同植被覆蓋類型區(qū)年最大NDVI呈現(xiàn)逐年增長趨勢.平原旱地、低覆蓋草地、有林地增幅較大,達0.1～0.2,疏林地在2004年后保持高速增長,增幅 0.208,增長了 77.6%.各類型區(qū) NDVI逐月變化趨勢與鶯落峽徑流量基本一致,生長季節(jié)集中在5～10月,峰值集中在7月.其中,灌木林地非生長季節(jié) NDVI最高,幾乎是其他類型的2～3倍,對改善生態(tài)環(huán)境意義重大,而平原旱地在7月收獲季節(jié)一過,NDVI迅速降低,對生態(tài)環(huán)境不利.

4.2 張臨高盆地植被覆蓋度增長趨勢明顯,其中,甘州區(qū)生態(tài)治理成效最顯著.研究區(qū)年最大植被覆蓋度線性擬合年均幅度 0.0063,生長季節(jié)平均植被覆蓋度在小幅波動中呈現(xiàn)穩(wěn)定增長趨勢,擬合優(yōu)度達0.74.張掖市甘州區(qū),年增長幅度集中在0.03～0.3,其后依次是臨澤縣和高臺縣.其中,高臺縣和臨澤縣北部被荒漠包圍的地區(qū),荒漠化趨勢明顯,局部年下降幅度達 0.1～0.3,年降低百分率達20%～60%.

4.3 從生長季節(jié)植被覆蓋度變化規(guī)律來看,植被7～10月覆蓋度呈現(xiàn)明顯增長趨勢,峰值從6～7月延遲到7～8月,2007年達0.3924,植被覆蓋度峰值提高、生長期延長,表明區(qū)域生態(tài)環(huán)境有所改善.從年最大植被覆蓋度分級統(tǒng)計來看,植被覆蓋度分級組成結(jié)構(gòu)整體呈現(xiàn)優(yōu)化趨勢,極低覆蓋度植被逐漸轉(zhuǎn)化為低覆蓋度植被,2007年相比2000年降低了25%以上,高覆蓋度植被比例的10年間增長約 16%,可以看出張掖市節(jié)水灌溉、退耕還林還草等一系列工程取得了不錯的成效.

4.4 根據(jù)關(guān)聯(lián)度模型分析結(jié)果,正相關(guān)驅(qū)動因子主要有降水量、代表農(nóng)業(yè)科技進步的用電量及化肥施用折純量,代表社會經(jīng)濟發(fā)展的第一產(chǎn)業(yè)產(chǎn)值、人均生產(chǎn)總值,和徑流性水資源鶯落峽徑流量等;負相關(guān)驅(qū)動因子主要有農(nóng)作物播種面積、蒸發(fā)量和第二產(chǎn)業(yè)產(chǎn)值.氣象水文驅(qū)動因子主要是與植被生長可利用水資源相關(guān)的降雨量、蒸發(fā)量、徑流量等,與3種植被覆蓋度最大關(guān)聯(lián)度分別為 0.9071、-0.8264、0.7553,社會經(jīng)濟驅(qū)動因子主要是與農(nóng)業(yè)生產(chǎn)相關(guān)的農(nóng)作物播種面積、第一產(chǎn)業(yè)產(chǎn)值、農(nóng)業(yè)科技水平等,最大關(guān)聯(lián)度分別為-0.8052、0.7810、0.8121.

[1] 張北贏,徐學(xué)選,劉文兆,等.黃土丘陵區(qū)不同土地利用的土壤水分灰色關(guān)聯(lián)度 [J]. 生態(tài)學(xué)報, 2008,28(1)：361-366.

[2] 鄧聚龍.灰色系統(tǒng)基本方法 [M]. 湖北：華中理工大學(xué)出版社, 1987.

[3] 李宏艷.關(guān)于灰色關(guān)聯(lián)度計算方法的研究 [J]. 系統(tǒng)工程與電子技術(shù), 2004,26(9)：1231-1234.

[4] 虞亞平,王冠中,李大治.廣義灰色關(guān)聯(lián)度的簡便計算方法 [J].南通大學(xué)學(xué)報(自然科學(xué)版), 2008,7(002)：85-90.

[5] 何滿喜.灰色綜合關(guān)聯(lián)度及其應(yīng)用 [J]. 內(nèi)蒙古師范大學(xué)學(xué)報(自然科學(xué)版), 1999,28(3)：170-173.

[6] 曾光明,曾北危.環(huán)境影響綜合評價的灰色關(guān)聯(lián)分析方法 [J].中國環(huán)境科學(xué), 1995,15(4)：247-251.

[7] 黨耀國,劉思峰,劉 斌,等.灰色斜率關(guān)聯(lián)度的改進 [J]. 中國工程科學(xué), 2004,6(3)：41-44.

[8] 王乃舉,周濤發(fā).礦業(yè)城市環(huán)境經(jīng)濟系統(tǒng)耦合評價—以安徽銅陵市為例 [J]. 中國環(huán)境科學(xué), 2012,32(7)：1339-1344.

[9] 曹明霞.灰色關(guān)聯(lián)度模型正負性問題的研究及其改進 [J]. 系統(tǒng)工程與電子技術(shù), 2008,30(6)：1086-1088.

[10] 姜朋輝,趙銳鋒,趙海莉,等.1975年以來黑河中游地區(qū)土地利用/覆被變化時空演變 [J]. 生態(tài)與農(nóng)村環(huán)境學(xué)報, 2012,28(5)：

[11] 甘肅年鑒編委會.甘肅年鑒：1991,2011 [M]. 北京：中國統(tǒng)計出版社, 1991,2011.

[12] 國家氣象信息中心,黑河流域逐日氣象觀測數(shù)據(jù)集 [Z]. 中國氣象局, 2010.

[13] Janssens Greet,馬明國.中國地區(qū)長時間序列 SPOT Vegetation植被指數(shù)數(shù)據(jù)集 [Z]. 2006.

[14] 甘肅年鑒編委會.甘肅年鑒：1999～2009 [M]. 北京：中國統(tǒng)計出版社, 1999-2009.

[15] 李 麗,童立強,李小慧.基于植被覆蓋度的石漠化遙感信息提取方法研究 [J]. 國土資源遙感, 2010,84(2)：59-62.

[16] 周兆葉,儲少林,王志偉,等.基于 NDVI的植被覆蓋度的變化分析-以甘肅省張掖市甘州區(qū)為例 [J]. 草業(yè)科學(xué), 2008,25(12)：23-29.

[17] Fletcher R S, Escobar D E, Skaria M. Evaluating airborne normalized difference vegetation index imagery for citrus orchard surveys [J]. Hort Technology, 2004,14(1)：91-94.

[18] 劉 琳,姚 波.基于 NDVI象元二分法的植被覆蓋變化監(jiān)測[J]. 農(nóng)業(yè)工程學(xué)報, 2010,26(S1)：230-234.

[19] 杜子濤,占玉林,王長耀.基于NDVI序列影像的植被覆蓋變化研究 [J]. 遙感技術(shù)與應(yīng)用, 2008,23(1)：47-50.

[20] 李苗苗.植被覆蓋度的遙感估算方法研究 [D]. 北京：中國科學(xué)院, 2003：49-55.

[21] 楊延征.基于SPOT-VGT NDVI的陜北植被覆蓋時空變化 [J].應(yīng)用生態(tài)學(xué)報, 2012,23(7)：1897-1903.

[22] 程國棟.黑河流域水-生態(tài)-經(jīng)濟系統(tǒng)綜合管理研究 [M]. 北京：科學(xué)出版社, 2009：559.

Changes in the ecological environment and there determining factors in the middle Heihe River Basin.

ZHOU Sha, HUANG Yue-fei*, WANG Guang-qian
(State Key Laboratory of Hydroscience and Engieering, Tsinghua University, Beijing 100084, China). China Environmental Science, 2014,34(3)：766～773

This paper analyzed several aspects of the spatio-temporal variations in the vegetation cover in the middle Heihe River Basin, and developed a grey correlation model of ecosystem to identify the determining factors. The results demonstrated that 1) the annual maximum normalized difference vegetation index (NDVI) increased considerably in upland, grassland with low vegetation cover, and forested area, by about 0.1～0.2; while sparse woodland has maintained a rapid growth since 2004with the annual maximum NDVI increased by 0.21; 2) the annual maximum vegetation cover for the study area as a whole has increased by 0.0063per annum, the change is most significant in Ganzhou District, where the annual growth was mostly around 0.03～0.3, however, for some oasis-desert transition zones in Linze and Gaotai, the vegetation cover decreased by 0.1～0.3; 3) the vegetation cover from July to October showed significant growth trends, the peak increased and delayed about one month, reaching 0.39in 2007; the hierarchical structure of the annual maximum vegetation cover optimized, the lowest vegetation cover gradually transformed into low vegetation cover, decreased by 25% from 2000to 2007, and the highest vegetation cover increased approximately 16% from 1998to 2008; 4) the results of correlation analysis showed that the main meteorological and hydrological factors were rainfall, evaporation and runoff ,with the maximum correlations of 0.91, -0.83, 0.76, respectively, and that the main socio-economic factors were agricultural acreage, primary industry and the level of agricultural science and technology, with the maximum correlation were -0.81, 0.78, 0.81, respectively.

Heihe River Basin；vegetation cover；spatiotemporal variation；grey correlation model

X171

：A

：1000-6923(2014)03-0766-08

周 沙(1991-),女,湖北孝感人,清華大學(xué)博士研究生,主要從事干旱區(qū)生態(tài)環(huán)境與水資源管理研究.發(fā)表論文2篇.

2013-07-29

國家自然科學(xué)基金(91125018);國家“十二五”科技支撐計劃重點項目(2013BAB05B03);水利部公益性行業(yè)科研專項項目(201301081)

* 責(zé)任作者, 教授, yuefeihuang@tsinghua.edu.cn