青海湖流域植被動(dòng)態(tài)變化驅(qū)動(dòng)力及空間粒度效應(yīng)

周美林 ,劉家宏 ,劉希勝 ,王亞琴 (1.中國水利水電科學(xué)研究院,流域水循環(huán)模擬與調(diào)控國家重點(diǎn)實(shí)驗(yàn)室,北京 100038；2.中國科學(xué)院地理科學(xué)與資源研究所,北京 100101；3.水利部數(shù)字孿生流域重點(diǎn)實(shí)驗(yàn)室,北京 100038；.青海省水文水資源測(cè)報(bào)中心,青海 西寧 810001；5.中國國土勘測(cè)規(guī)劃院,北京 100081)

植被覆蓋度變化與區(qū)域氣候條件、生態(tài)環(huán)境、土地利用和人類活動(dòng)等因素密切相關(guān)[1-3].多采用相關(guān)性分析[4]、多元線性回歸法[5]等分析氣候和人類活動(dòng)因素的影響.青海湖流域是維系青藏高原生態(tài)安全的重要屏障[6],生態(tài)系統(tǒng)脆弱.氣候、地表因素和人類活動(dòng)是青海湖流域植被覆蓋變化的主要驅(qū)動(dòng)力.氣候變化改變陸地生態(tài)系統(tǒng)格局和分布范圍[7],地形地貌和人類活動(dòng)對(duì)植被覆蓋度的影響具有顯著的空間分異[8-9],且各因素影響存在顯著的交互效應(yīng)[10-11].

已有研究[12-13]分析了2000 年后,不同降水、氣溫、高程和植被類型下青海湖流域植被覆蓋度變化特征.2000 年以前青海湖流域植被覆蓋狀況認(rèn)識(shí)不足.同時(shí),針對(duì)植被覆蓋度變化驅(qū)動(dòng)研究,現(xiàn)有研究?jī)H定性分析不同氣候帶和地形地貌上的差異[14],對(duì)影響因素相對(duì)貢獻(xiàn)的量化和交互效應(yīng)的認(rèn)識(shí)仍不清楚.地理探測(cè)器是廣泛使用的多因素交互作用模型[15-16],可識(shí)別并量化因子及其交互作用對(duì)變量空間變異的貢獻(xiàn)[17-18].空間粒度影響變量特征及相互關(guān)系[19],氣候、地形和人類活動(dòng)等因素對(duì)植被覆蓋度的影響與空間粒度密切相關(guān).現(xiàn)有研究多為單粒度分析,對(duì)尺度效應(yīng)研究較少.

本文利用Google Earth Engine 平臺(tái)以Landsat遙感影像為基礎(chǔ),分析1986～2020 年青海湖流域植被覆蓋度時(shí)空變化特征,基于Pearson相關(guān)系數(shù),多元線性回歸模型及多尺度地理探測(cè)器模型,揭示降水、氣溫、地表?xiàng)l件(地形、地貌、土壤、植被、高程、坡度、坡向、凍土、距河/青海湖距離)和人為因素對(duì)流域植被覆蓋度時(shí)空變化的驅(qū)動(dòng)機(jī)制及相對(duì)貢獻(xiàn),并探討空間粒度對(duì)驅(qū)動(dòng)因素的影響.

1 數(shù)據(jù)和方法

1.1 研究區(qū)概況

青海湖是我國最大的內(nèi)陸咸水湖,青海湖流域位于97°50′E～101°20′E,36°15′N～38°20′N,之間,流域總面積約29661km2,海拔在3196～5290m 之間;多年平均降水量為379mm,多集中在6～8 月;多年平均氣溫為-0.3℃,7月氣溫最高,1月氣溫最低,12月到次年3 月湖面結(jié)冰[20].青海湖區(qū)年蒸發(fā)量為895mm,夏季占年蒸發(fā)量的60%,冬春季蒸發(fā)量最小[21].青海湖流域是重要的凍土分布區(qū),西北部為多年永久凍土區(qū),中部及環(huán)湖區(qū)域?yàn)榧竟?jié)性凍土區(qū).主要植被類型為草甸、草原、荒漠、灌叢和林地,其中草甸占流域面積45%,草原占流域面積18%.土壤類型有草甸土、黑鈣土、栗鈣土、沼澤土和風(fēng)沙土等.青海湖流域地理位置如圖1 所示.

圖1 青海湖流域地理位置Fig.1 Geographical location of Qinghai Lake basin(QHLB)

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

本文使用數(shù)據(jù)有遙感影像、氣候、地形地貌、植被、土壤、凍土分區(qū)和人類足跡.遙感影像來自美國陸地衛(wèi)星Landsat 4TM、Landsat 5TM、Landsat 7ETM+和Landsat 8OLI 2 級(jí)產(chǎn)品,已完成輻射校正和幾何校正,時(shí)間分辨率為16d,空間分辨率為30m,計(jì)算時(shí)段為1986～2020 年.高程數(shù)據(jù)來自SRTM 數(shù)字高程模型(DEM),空間分辨率為30m.坡度和坡向基于DEM 數(shù)據(jù),利用ArcGIS 軟件空間分析工具獲得.氣象數(shù)據(jù)采用時(shí)空三極環(huán)境大數(shù)據(jù)中心發(fā)布的中國1km 分辨率降水和氣溫?cái)?shù)據(jù),該數(shù)據(jù)集是在CRU和World Clim發(fā)布的全球0.5°氣候數(shù)據(jù)基礎(chǔ)上,基于Delta 空間降尺度方案生成,并通過了496 個(gè)氣象站監(jiān)測(cè)結(jié)果的驗(yàn)證[22-23].青藏高原凍土分布[24]基于改進(jìn)的中分辨率成像光譜儀(MODIS)和地表溫度凍結(jié)和融化指數(shù)生成,空間分辨率為937m.土壤和植被數(shù)據(jù)來自中國科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心.土壤采用傳統(tǒng)的“土壤發(fā)生類”系統(tǒng),將中國土壤分為12 個(gè)土綱,61 個(gè)土類,227 個(gè)亞類,其基礎(chǔ)數(shù)據(jù)是1:100萬中國土壤圖和中國土種志.植被共分為11個(gè)植被類型組、54 個(gè)植被型,涵蓋了796 個(gè)群系和亞群系,空間分辨率為1km.

采用人類足跡評(píng)估人類活動(dòng)強(qiáng)度對(duì)青海湖流域生態(tài)系統(tǒng)的影響.人類足跡指數(shù)法[25]是常用的人類足跡評(píng)估指標(biāo),綜合考慮了人口密度、土地利用變化、交通通達(dá)性和電力基礎(chǔ)設(shè)施等四類空間數(shù)據(jù).通過設(shè)置各數(shù)據(jù)圖層的緩沖區(qū)和賦值權(quán)重,經(jīng)過疊加分析和空間歸一化,得到人類足跡數(shù)據(jù).青藏高原人類足跡數(shù)據(jù)考慮人口密度、土地利用、放牧密度、夜間燈光、鐵路和道路等指標(biāo)表征人類活動(dòng)強(qiáng)度,采用對(duì)數(shù)法對(duì)不同類型數(shù)據(jù)賦值,賦值范圍為1～10.詳細(xì)賦值過程可參考文獻(xiàn)[26].該數(shù)據(jù)集包含1990、1995、2000、2005、2010、2015 和2017 年共7 期數(shù)據(jù),空間分辨率為1km.

1.3 研究方法

1.3.1 植被覆蓋度計(jì)算 利用Landsat 影像的近紅外波段和紅色波段進(jìn)行植被歸一化指數(shù)(NDVI)計(jì)算,隨后基于中值法構(gòu)建月度NDVI[27],基于生長(zhǎng)季(4 月～10 月)數(shù)據(jù)構(gòu)建年際NDVI 數(shù)據(jù)集.利用年際NDVI 基于像元二分模型計(jì)算年際植被覆蓋度(FVC).具體計(jì)算公式如下:

式中NIR 和Red 是近紅外波段和紅色波段的反射率,NDVImin為青海湖流域?qū)?yīng)5%頻次的NDVI值,NDVImax為對(duì)應(yīng)95%頻次的NDVI 值.

基于《水土流域綜合治理技術(shù)規(guī)范》,將青海湖流域FVC 分類,其中FVC≤0.1 為裸地,0.1

1.3.2 MK 趨勢(shì)檢驗(yàn) Mann-Kendall(MK)[28]方法是一種非參數(shù)檢驗(yàn)方法,廣泛應(yīng)用于時(shí)間序列趨勢(shì)變化分析.對(duì)于時(shí)間序列X={x1, x2, …, xn},假定序列內(nèi)各變量相互獨(dú)立,構(gòu)造檢驗(yàn)統(tǒng)計(jì)量S 如下:

對(duì)于給定的置信區(qū)間α,若Z≥Z1-a/2則認(rèn)為在α置信水平上序列有明顯的上升(Z>0)或下降趨勢(shì)(Z<0),其中|Z|>2.576 說明變化趨勢(shì)極顯著(α=0.01,P<0.01),2.576>|Z|>1.96 說明變化趨勢(shì)顯著(α=0.05,0.01

1.3.3 地理探測(cè)器 地理探測(cè)器是常用的分析變量空間變異驅(qū)動(dòng)機(jī)制的統(tǒng)計(jì)方法,廣泛應(yīng)用于土地利用和植被覆蓋等驅(qū)動(dòng)機(jī)制分析[29].采用地理探測(cè)器中的因子探測(cè)器和交互探測(cè)器進(jìn)行青海湖流域FVC 空間變化驅(qū)動(dòng)機(jī)制分析.因子探測(cè)器分析驅(qū)動(dòng)因素對(duì)研究區(qū)域FVC 空間變化影響程度的大小,用q 量化.交互探測(cè)器用來探測(cè)不同因子交互作用對(duì)FVC 空間變異的影響,分為非線性增強(qiáng),雙因子增強(qiáng),非線性減弱,雙因子減弱和獨(dú)立.q 計(jì)算公式如式(2).采用t 檢驗(yàn)法對(duì)相關(guān)顯著性進(jìn)行檢驗(yàn),設(shè)置顯著性標(biāo)準(zhǔn)為α=0.05.

式中:q 表征驅(qū)動(dòng)因素對(duì)植被覆蓋度變化解釋力的大小;h 為1,2,…,L 為變量的分類(區(qū))數(shù);Nh和N 分別指分類h 和全區(qū)的分類數(shù);σ2h 和σ2分別為某分層h 和全區(qū)的方差總和.

基于R語言GD包構(gòu)建地理探測(cè)器模型,在模型構(gòu)建過程中,采用自然斷點(diǎn)法、等距離散法和等頻率法對(duì)連續(xù)變量進(jìn)行空間離散,設(shè)置類別數(shù)量為4～10.參考已有文獻(xiàn)的基礎(chǔ)上,結(jié)合實(shí)際情況,選取地表因素(高程、坡度、坡向、土壤、植被、距水系距離、距青海湖距離、凍土分區(qū)),氣候因素(生長(zhǎng)季降水、生長(zhǎng)季氣溫,降水集度、降水重心)和人類足跡指數(shù)作為FVC 空間變化的驅(qū)動(dòng)因素.

空間粒度對(duì)生態(tài)系統(tǒng)結(jié)構(gòu)和影響因素空間分布具有重要影響[30],以青海湖流域FVC 為研究對(duì)象,分析不同空間粒度下氣候、人類足跡及地表因素對(duì)FVC 解釋力的變化.考慮已有數(shù)據(jù)的空間分辯率,以1km 為間隔,分別計(jì)算1～15km 分辨率下各因素對(duì)植被覆蓋度空間變化的解釋力.以最大程度表征因素交互作用和各因素解釋力90%分位數(shù)最大為標(biāo)準(zhǔn),確定流域植被覆蓋度空間變化最佳分析粒度,并作為后續(xù)驅(qū)動(dòng)因素分析的基礎(chǔ).

1.3.4 降水結(jié)構(gòu)計(jì)算 降水重心和降水集度是降水結(jié)構(gòu)的表征.降水重心(PC)反應(yīng)某區(qū)域年內(nèi)降水集中時(shí)間.降水集度(PCI)指年內(nèi)降水特征在年際尺度上的分異.PCI 越大說明年內(nèi)降水分配越集中,PCI越小說明年內(nèi)降水分配越均勻.具體計(jì)算公式如下,式中Pi為第i 個(gè)月降水,P 為年際降水量.

采用t 檢驗(yàn)法對(duì)相關(guān)顯著性進(jìn)行檢驗(yàn),設(shè)置顯著性標(biāo)準(zhǔn)為α=0.05.

1.3.5 相關(guān)性分析 采用Pearson 相關(guān)系數(shù)評(píng)估不同時(shí)期氣候因素與FVC 變化相關(guān)性.采用t 檢驗(yàn)法對(duì)相關(guān)顯著性進(jìn)行檢驗(yàn),設(shè)置顯著性標(biāo)準(zhǔn)為α=0.05.

2 結(jié)果與討論

2.1 FVC 年際變化

由圖2 可見,1986～2020 年青海湖流域平均FVC呈波動(dòng)增加趨勢(shì),其中1986 年流域平均FVC 最大,1990 年FVC 最小.MK 突變檢驗(yàn)發(fā)現(xiàn)2001 年流域平均FVC 發(fā)生突變(圖2),表現(xiàn)為流域FVC 整體增加.比較突變前后期FVC 變化,發(fā)現(xiàn)1986～2000 年流域FVC 呈減小趨勢(shì),均值為0.33,2001～2020 年流域FVC 呈增加趨勢(shì),均值為0.37.2001 年后青海湖流域平均植被覆蓋度增加幅度為12%,流域植被狀況由退化轉(zhuǎn)為改善.

圖2 青海湖流域平均FVC 變化Fig.2 Average FVC changes in QHLB

由圖3 和表1 可見,1986～2000 年,僅3%區(qū)域FVC 發(fā)生極顯著/顯著變化,其中環(huán)青海湖北部、南部和東部FVC顯著減小,流域西北部FVC顯著增加.2001～2020 年流域FVC 極顯著/顯著變化區(qū)域占比16%,其中環(huán)青海湖西北部、東部及布哈河流域上游和下游沿河地區(qū)顯著增加(占比15%),流域北部沙柳河上游及布哈河中部零星區(qū)域顯著減小.2001 年后青海湖流域FVC 顯著減小區(qū)域減小,顯著增加區(qū)域增大.

表1 不同時(shí)期青海湖流域FVC 變化空間特征Table 1 Spatial characteristics of FVC changes in QHLB in different periods

圖3 不同時(shí)期青海湖流域FVC 空間變化Fig.3 The spatial variations of FVC in QHLB in different periods

2.2 FVC 類型空間變化

以1990 年和2020 年表征突變前后期青海湖流域FVC 類型變化(圖4).相較而言,2020 年青海湖流域中等、中高和高覆蓋度植被面積增加,表現(xiàn)為從環(huán)湖區(qū)域向流域中部延伸的特征.FVC 類型變化表現(xiàn)為裸地和低覆蓋度植被面積由55%減小為25%,中等覆蓋度植被面積由6%增加為30%,中高和高覆蓋度植被面積由1%增加為16%.同時(shí),1990～2020年湖流域FVC類型轉(zhuǎn)移關(guān)系為裸地轉(zhuǎn)為低覆蓋度植被,低覆蓋度植被轉(zhuǎn)為中低、中等覆蓋度植被,中等覆蓋度植被轉(zhuǎn)為中高覆蓋度植被.從流域FVC類型變化角度分析,FVC類型變化呈向好態(tài)勢(shì),中高、高等覆蓋度植被類型區(qū)域的FVC顯著增加,青海湖流域生態(tài)環(huán)境顯著改善.

圖4 不同時(shí)期青海湖流域植被覆蓋類型空間格局Fig.4 Spatial patterns of FVC types in QHLB in the different periods

2.3 FVC 時(shí)間變化驅(qū)動(dòng)力

2.3.1 氣候因素 如圖5所示,1986～2000年降水呈減小趨勢(shì),氣溫呈增加趨勢(shì),青海湖流域氣候?yàn)榕苫?2001～2020 年,降水增加,氣溫持續(xù)增加,氣候呈暖濕化,因此2000 年以前青海湖流域植被退化和2000年后流域FVC 發(fā)生增加突變,與氣候變化趨勢(shì)密切相關(guān).然而流域平均FVC 與降水、氣溫相關(guān)系數(shù)均不顯著,這可能與氣候變化對(duì)植被覆蓋影響具有明顯的空間變異相關(guān).

圖5 流域平均生長(zhǎng)季降水和氣溫變化Fig.5 Average growing season precipitation and temperature changes in the land region of QHLB

以Pearson 相關(guān)系數(shù)為指標(biāo)分析氣候和人類活動(dòng)對(duì)FVC 變化影響的空間差異(圖6).青海湖流域?qū)儆谒疅嵯拗茀^(qū)域,水量和氣溫增加均對(duì)植被生長(zhǎng)為促進(jìn)作用.1986～2000 年環(huán)青海湖東北及南部FVC顯著減小,與降水呈正相關(guān),而與氣溫呈負(fù)相關(guān),相關(guān)系數(shù)為-0.5～-0.8;流域北部FVC 增加,與降水呈負(fù)相關(guān),與氣溫呈正相關(guān),相關(guān)系數(shù)為0.3～0.8.2001～2020年環(huán)青海湖東北部及布哈河上游和中游FVC 增加,降水、氣溫與該區(qū)域FVC 均為正相關(guān);流域北部FVC 顯著減小,降水、氣溫與該區(qū)域FVC 均為負(fù)相關(guān).因此,1986～2000 年青海湖流域東北及南部FVC減小受降水減小影響,而流域北部FVC 增加與氣溫增加關(guān)系密切,2001～2020 年流域FVC 顯著增加區(qū)域受降水和氣溫共同影響顯著.

圖6 不同時(shí)期氣候因素與FVC 相關(guān)系數(shù)空間分布Fig.6 Spatial distributions of correlation coefficients between climate factors and FVC in different periods

降水集度(PCI)和降水重心(PC)是降水結(jié)構(gòu)重要指標(biāo),表征降水分配均勻度和集中發(fā)生時(shí)間.降水結(jié)構(gòu)變化對(duì)青藏高原植被生長(zhǎng)具有重要影響,分析不同時(shí)期PC、PCI 與流域FVC 變化的關(guān)系(圖6).1986～2000 年在青海湖流域FVC 減小區(qū)域,PC、PCI與FVC 均為負(fù)相關(guān),增加區(qū)域?yàn)檎嚓P(guān);2001～ 2020年環(huán)青海湖東部和南部FVC 增加區(qū)域,FVC 與PCI、PC 為負(fù)相關(guān)和正相關(guān),而在布哈河上游和沿河FVC增加區(qū)域,與PCI 和PC 為正相關(guān)和負(fù)相關(guān).以上結(jié)果表明2000 年前PCI、PI 對(duì)FVC 變化均為促進(jìn)作用,2000 年后作用方向相反,空間差異明顯.

2.3.2 人為因素 由圖7可知,1990～2017年青海湖流域平均人類活動(dòng)強(qiáng)度呈增加趨勢(shì),增加幅度不顯著,其中2000 年和2010 年是趨勢(shì)變化轉(zhuǎn)折點(diǎn).1990年和2017 年青海湖流域人類活動(dòng)強(qiáng)度空間變化格局基本一致,環(huán)湖地區(qū)、道路及流域水系周邊人類活動(dòng)強(qiáng)度增大,其它區(qū)域人類活動(dòng)強(qiáng)度基本無變化.人類足跡對(duì)青海湖流域FVC 影響空間分異明顯,主要影響區(qū)域?yàn)榄h(huán)湖地區(qū)、道路和水系周邊.

圖7 青海湖流域人類足跡時(shí)空變化Fig.7 Spatiotemporal changes of human footprint in QHLB

2.3.3 氣候和人類活動(dòng)對(duì)流域平均 FVC 的貢獻(xiàn) 以生長(zhǎng)季氣溫、生長(zhǎng)季降水和人類足跡指數(shù)為自變量,FVC 為因變量,對(duì)各變量歸一化處理后,采用多元線性回歸模型,構(gòu)建流域平均FVC 與氣候、人類活動(dòng)關(guān)系模型(FVC=0.16P+0.62T-0.09HF,P 為降水,T為氣溫,HF 為人類足跡).模擬結(jié)果表明氣候變化對(duì)流域平均FVC 為正向作用,人類足跡對(duì)流域平均FVC為抑制作用,各因素貢獻(xiàn)依次為氣溫>降水>人類足跡.考慮到人類足跡只能表征人類活動(dòng)對(duì)植被覆蓋狀況可能形成的負(fù)面效應(yīng),無法考慮生態(tài)恢復(fù)工程的正向效應(yīng).而多元線性擬合模型對(duì)流域平均FVC 變化解釋度為51%,仍有49%變化無法解釋,因此生態(tài)恢復(fù)工程在流域FVC 變化中具有重要貢獻(xiàn).

2.4 FVC 空間變化驅(qū)動(dòng)力

2.4.1 空間粒度影響分析 如圖8 所示,不同因素相對(duì)解釋力大小較為穩(wěn)定,從大到小依次為土壤、高程、植被、凍土、距青海湖距離和坡度.隨著空間粒度的增大,高程和凍土貢獻(xiàn)表現(xiàn)為先增加后減小,空間粒度為9km 時(shí)解釋力最大.其它地表因素對(duì)空間粒度的響應(yīng)不同,當(dāng)空間粒度大于臨界值時(shí),因素對(duì)流域FVC 空間變化貢獻(xiàn)消失.其中坡向和距水系距離的臨界粒度為3km,距青海湖距離的臨界粒度為5km,坡度和植被的臨界粒度為6km,土壤的臨界粒度為7km.

圖8 不同空間尺度下驅(qū)動(dòng)因素對(duì)FVC 空間變化的解釋力Fig.8 The explanatory power(q-value) of driving factors on spatial variations in VFC at different spatial scales

氣候和人為因素對(duì)FVC 貢獻(xiàn)隨粒度變化與地表因素一致,空間粒度變化不影響因素間的相對(duì)解釋力.各因素對(duì)FVC 解釋力從大到小依次為氣溫、人類足跡、降水集度、降水量和降水重心,隨著粒度增加,氣溫和人為因素貢獻(xiàn)先增加后減小,降水重心臨界粒度為4km,降水量和降水集度臨界粒度為8km.

2.4.2 空間驅(qū)動(dòng)力分析 綜合因素解釋力90%分位數(shù)和因素交互作用,以6km 空間粒度作為空間驅(qū)動(dòng)力分析的基礎(chǔ).該空間粒度下坡向和降水集度對(duì)FVC 空間變化解釋力消失,其它因素對(duì)FVC 空間變化解釋均通過顯著性檢驗(yàn)(α=0.05),解釋力從大到小依次為生長(zhǎng)季氣溫(0.41)、高程(0.34)、植被(0.25)、土壤(0.23)、降水集度(0.17)、人類足跡(0.15)、生長(zhǎng)季降水(0.13)、距水系距離(0.12)、距青海湖距離(0.11)、坡度(0.06)和凍土分區(qū)(0.05).因此該空間粒度下青海湖流域FVC 空間變化主要由氣候、地形地貌、土壤和植被共同決定.

由圖9 可見,各因素交互作用對(duì)FVC 貢獻(xiàn)均大于因素單獨(dú)作用,交互作用均為增強(qiáng)效應(yīng).生長(zhǎng)季氣溫、高程與其它因素交互作用顯著,與距水系距離、坡度、人類足跡、降水量、降水集度、土壤和植被交互作用后,貢獻(xiàn)均顯著增大,與距水系距離交互作用解釋力最大(0.61).由于氣溫是外部環(huán)境驅(qū)動(dòng)力,交互作用主要表現(xiàn)為氣溫增加對(duì)其它因素的影響,進(jìn)而改變流域植被覆蓋狀況.青海湖流域位于高寒山區(qū),高程變化較大,高程變化影響區(qū)域地理?xiàng)l件,對(duì)交互作用強(qiáng)度具有地帶控制性.與1986～ 2000 年相比,2001 年后高程、氣溫與坡度、土壤、植被和距水系距離交互作用解釋力略有提高,其它因素貢獻(xiàn)無明顯變化,表明氣溫升高對(duì)青海湖流域植被狀況改善同時(shí)具有間接促進(jìn)作用.

圖9 青海湖流域FVC 各驅(qū)動(dòng)因子交互探測(cè)解釋力(q)Fig.9 Exploratory power(q) of interactions among driving factors of FVC in QHLB

2.5 討論

本文比較了1986～2000 年和2001～2020 年2 個(gè)時(shí)段青海湖流域FVC 變化,發(fā)現(xiàn)2001 年后FVC 顯著增加區(qū)域?yàn)榄h(huán)青海湖西北部及布哈河沿河區(qū)域.

由圖10 可見,2001 年后環(huán)青海湖及布哈河降水和氣溫均增加,且氣溫增加幅度明顯大于其它區(qū)域.青海湖流域不同海拔區(qū)域植被生長(zhǎng)對(duì)氣溫和降水的敏感性不同[31].環(huán)青海湖和布哈河中下游植被生長(zhǎng)受水分和熱量限制區(qū)域,植被對(duì)降水和氣溫變化較其他區(qū)域敏感,因此該區(qū)域植被覆蓋度增加顯著.值得注意,布哈河上游為氣溫限制區(qū)域,2001 年后氣溫增加,植被生長(zhǎng)能量條件改善,FVC 顯著增加.

圖10 青海湖流域2001～2020 年氣溫和降水變化趨勢(shì)Fig.10 Trends in temperature and precipitation changes in the QHLB from 2001 to 2020

青藏高原氣溫升高引起的生長(zhǎng)季溫度仍低于植被生長(zhǎng)最適溫度[32],氣溫對(duì)青海湖流域FVC 空間變化解釋力最大,氣候變暖背景下氣溫仍是青海湖流域植被生長(zhǎng)的主控因素.降水格局變化即降水量和降水結(jié)構(gòu),也是氣候變化的重要方面,降水量對(duì)青海湖流域FVC 空間變化解釋力較低,然而考慮降水結(jié)構(gòu)(降水集度)時(shí),降水對(duì)流域FVC 空間解釋力達(dá)0.3,僅次于氣溫和高程,說明降水結(jié)構(gòu)對(duì)青海湖流域植被生長(zhǎng)具有重要影響.同時(shí),相關(guān)研究[33]表明降水量和降水集度變化對(duì)藏北高原草地生產(chǎn)力有重要影響.距河流水系距離和青海湖距離對(duì)流域FVC 變化解釋力低于氣溫和降水.距水系/青海湖距離,與區(qū)域土壤水、地表水和地下水連通性相關(guān),表明地下水對(duì)青海湖流域FVC 影響不顯著,可能由于研究區(qū)域凍土層存在減弱了地下水和土壤水間的聯(lián)系.

地理探測(cè)器因素交互作用結(jié)果表明各因子對(duì)青海湖流域植被覆蓋空間變化均存在交互效應(yīng),與其它區(qū)域研究結(jié)果一致.氣溫和高程是青海湖流域FVC 空間變化驅(qū)動(dòng)因素交互作用的主導(dǎo),這與青海湖流域地理位置密切相關(guān).氣溫∩距水系距離交互作用解釋力最大(0.61),其次為氣溫∩降水(0.58),氣溫∩降水集度(0.56).氣溫與距水系距離、降水交互后解釋力增大,表明水熱條件同時(shí)控制著青海湖流域FVC 空間變化.人類足跡與高程、氣溫交互作用后對(duì)植被覆蓋空間變化解釋力由0.15 增大到0.52和0.58,表明海拔和氣溫對(duì)人類活動(dòng)的限制.凍土∩氣溫、海拔對(duì)FVC 變化貢獻(xiàn)從0.05 增大為0.41 和0.35,一方面由于凍土分布與高程相關(guān),另一方面氣溫升高,引起永久凍土和季節(jié)凍土退化,流域土壤水減小[34],影響流域植被覆蓋狀況.

3 結(jié)論

3.1 青海湖流域植被狀況整體改善.2001 年為植被變化突變期,FVC 由減小轉(zhuǎn)為增加.空間上表現(xiàn)為環(huán)青海湖西北、布哈河上游FVC 顯著增加,流域東北部、青海湖東部沙地零星區(qū)域FVC 顯著減小.類型上表現(xiàn)為裸地、低覆蓋植被減小,中等和中高覆蓋植被增加.

3.2 青海湖流域FVC 時(shí)空變化驅(qū)動(dòng)因素不同.時(shí)間上,由氣候變化和生態(tài)恢復(fù)工程驅(qū)動(dòng).空間上,由氣溫(0.41)、高程(0.34)、降水(0.30)、地形地貌和土壤植被驅(qū)動(dòng).各因素均為顯著交互的增強(qiáng)效應(yīng),其中氣溫和高程是主控因素,與距水系距離、降水量、降水集度及人類足跡作用后,解釋力顯著增大.

3.3 空間粒度影響氣候、地表?xiàng)l件和人為因素對(duì)流域FVC 空間變化的貢獻(xiàn).空間粒度影響因素解釋力,對(duì)因素間相對(duì)貢獻(xiàn)影響不顯著.綜合考慮各因素貢獻(xiàn)和交互作用過程,青海湖流域空間驅(qū)動(dòng)力最佳分析粒度為6km.