曹 會(huì),劉立文,李雨珂,溫秋園
1.中國(guó)地質(zhì)調(diào)查局牡丹江自然資源綜合調(diào)查中心,黑龍江 牡丹江 157021;2.山西農(nóng)業(yè)大學(xué) 資源環(huán)境學(xué)院,山西 太原 030031
在陸地生態(tài)系統(tǒng)中植被居于主體地位并且影響著全球環(huán)境的變化,通過分析植被局部的動(dòng)態(tài)變化可以為改善區(qū)域生態(tài)、保護(hù)物種多樣性和可持續(xù)發(fā)展提供科學(xué)指導(dǎo),因此植被覆蓋的變化分析成為當(dāng)前生態(tài)學(xué)研究的重點(diǎn)[1-2].氣候、地表植被的變化監(jiān)測(cè)研究中,植被覆蓋是非常重要的指標(biāo)及參數(shù)[3-5].而景觀格局信息,可以揭示不同類型斑塊的空間布局和結(jié)構(gòu)特征,通過研究植被的景觀格局動(dòng)態(tài)變化可以反映植被在人為因素和環(huán)境因素共同作用下的變化特征,對(duì)進(jìn)一步評(píng)估地表植被以及改善區(qū)域生態(tài)環(huán)境具有重要意義[6-9].
目前的研究主要是通過景觀組分的類型、數(shù)目以及空間的分布與配置從表面上沒有規(guī)律的景觀中發(fā)掘潛在的有意義的生態(tài)變化特征,進(jìn)而用于指導(dǎo)生態(tài)環(huán)境建設(shè)[10-12].如蘇艷琴等對(duì)福建省三明市沙縣進(jìn)行研究,利用線性光譜混合模型獲取植被覆蓋度并分析其變化及景觀格局[13];王剛等利用地表溫度、植被覆蓋度和歸一化濕度指數(shù)等信息,結(jié)合生態(tài)學(xué)方法研究了廣州市不同植被對(duì)地表溫度的調(diào)節(jié)作用[14];黃婷等利用景觀格局分析與逐步回歸的統(tǒng)計(jì)方法,研究2005—2010年錫林郭勒盟在流域尺度上的景觀格局變化對(duì)土壤保持水平的影響[15];徐凱健等分析了不同時(shí)期長(zhǎng)汀縣植被覆蓋及其空間格局的時(shí)空變化,研究人為影響與政策規(guī)劃對(duì)植被覆蓋及景觀格局的影響[16];王靜等定量研究了2000—2010年京津冀地區(qū)植被覆蓋及其景觀格局的動(dòng)態(tài)變化,評(píng)估了城市化對(duì)植被的干擾過程和生態(tài)環(huán)境的影響[17].
晉城市地處黃土丘陵區(qū)東南邊緣,眾多復(fù)雜地形地貌相互交織,導(dǎo)致該區(qū)植被覆蓋狀況差異較大,水土流失嚴(yán)重.基于此,本研究使用FRAGSTATS運(yùn)算軟件[18],通過利用2006—2015年MODIS的MOD13Q1數(shù)據(jù),得到晉城市植被覆蓋指數(shù),并將其劃分為5個(gè)等級(jí),然后通過GIS分析的方法研究植被覆蓋在不同高程、坡度和地理單元下的變化以及分析不同植被覆蓋等級(jí)下斑塊的景觀格局變化,根據(jù)該地區(qū)不同景觀指數(shù)的年際變化情況,研究晉城市植被覆蓋時(shí)空變化的特點(diǎn)和規(guī)律,從而為晉城市生態(tài)環(huán)境建設(shè)提供科學(xué)依據(jù),達(dá)到生態(tài)環(huán)境的可持續(xù)發(fā)展.
晉城市位于山西省東南部(圖1),山西與河南交界處,地理坐標(biāo)范圍:北緯35°11′—36°04′,東經(jīng)111°55′—113°07′[19].全市處于晉城盆地中,被太行、王屋、中條山所環(huán)繞,從北到南地勢(shì)逐漸降低,構(gòu)成一個(gè)簸箕狀,沁河、丹河交錯(cuò)其中,主要呈現(xiàn)出河谷平川、丘陵和中低山3種地貌,其中山地所占面積最大,河谷平川和丘陵次之.平原區(qū)地勢(shì)開闊,植被資源豐富,是工農(nóng)業(yè)生產(chǎn)的重要基地.晉城市交通便利,太焦、侯月鐵路縱穿全區(qū),公路四通八達(dá),交織縱橫[20].
圖1 研究區(qū)地貌及位置圖Fig.1 Topography and location map of the study area
晉城市的氣候特征是暖溫帶半濕潤(rùn)大陸性季風(fēng)氣候,春暖秋涼,冬冷夏熱,降雨主要集中在夏季.年平均日照時(shí)數(shù)2 563 h,年平均氣溫10.2~12℃,年平均降水量626~750 mm.較為濕潤(rùn)的氣候和復(fù)雜的地貌,導(dǎo)致晉城市植被較為豐富,其中林木覆蓋率達(dá)到33.6%,以占全市林地面積70%左右的天然次生林為主.另外,晉城市擁有一塊面積超過700 hm2的原始森林,是山西省唯一的原始森林.處于晉城市域的25×104hm2以上的天然牧草地,占土地面積近30%,自然植被豐富.
衛(wèi)重遙感數(shù)據(jù)下載于NASA所提供的共享數(shù)據(jù),數(shù)據(jù)分類為MOD13Q1產(chǎn)品,本研究所使用的是2006—2015年晉城市歸一化植被指數(shù)(NDVI)數(shù)據(jù),空間及時(shí)間分辨率分別為250 m和16 d.該數(shù)據(jù)已經(jīng)過輻射校正、大氣校正、云體掩膜等處理,可以更好地用于科學(xué)研究.用于本次研究的山西省DEM數(shù)據(jù)取自國(guó)家地理空間數(shù)據(jù)云(www.gscloud.cn),數(shù)據(jù)分辨率為30 m,依據(jù)此數(shù)據(jù)可計(jì)算出坡度等特征.
按照《生態(tài)環(huán)境狀況評(píng)價(jià)技術(shù)規(guī)范(HJ192—2015)》,獲取植被覆蓋指數(shù)(C),并結(jié)合前人研究成果[21-23],根據(jù)自然間斷分級(jí)法和晉城市地區(qū)植被覆蓋特征規(guī)律,將植被覆蓋等級(jí)分為5個(gè)級(jí)別:低植被覆蓋(C≤0.6)、較低植被覆蓋(0.6~0.7)、中度植被覆蓋(0.7~0.8)、較高植被覆蓋(0.8~0.9)及高植被覆蓋(C≥0.9)(圖2).
圖2 2006—2015年晉城市平均植被覆蓋分級(jí)圖Fig.2 Classification map of average vegetation coverage in Jincheng City during 2006-2015
景觀格局分析是一種研究景觀結(jié)構(gòu)特征和空間配置的方法,是景觀生態(tài)學(xué)研究中經(jīng)常使用的方法[24].在參照前人研究的基礎(chǔ)上[25-27],結(jié)合晉城市植被覆蓋特點(diǎn),選取斑塊類型尺度中斑塊類型面積(CA)、景觀級(jí)別中斑塊數(shù)量(NP)、斑塊密度(PD)、散布與并列指數(shù)(IJI)、景觀形狀指數(shù)(LSI)、香農(nóng)多樣性指數(shù)(SHDI)、香農(nóng)均勻性指數(shù)(SHEI)等作為評(píng)價(jià)指標(biāo).運(yùn)用FRAGSTATS3.3軟件進(jìn)行景觀指數(shù)計(jì)算并分析景觀的空間分布特征和變化特征.
景觀面積作為分析景觀格局變化的主要參數(shù),是計(jì)算景觀中其他指數(shù)的基礎(chǔ).景觀面積對(duì)反映不同類型斑塊的生物數(shù)量、生物種類、儲(chǔ)量和養(yǎng)分的差異具有重要的意義.通過分析圖3內(nèi)數(shù)據(jù)可知,2006—2015年,高植被覆蓋、較高植被覆蓋的景觀面積有所增加,低植被覆蓋、較低植被覆蓋和中度植被覆蓋景觀面積則在減少.其中,高植被覆蓋景觀面積增加最為明顯,面積增加116.07 km2,增長(zhǎng)幅度達(dá)12.24%,增速為13.24%/10 a;中度植被覆蓋景觀面積減少最多,達(dá)86.97 km2,減少幅度達(dá)9.17%,減速為8.50%/10 a;較高植被覆蓋景觀面積增加27.87 km2,增加幅度為2.93%,增速4.60%/10 a;低植被覆蓋景觀面積和較低植被覆蓋景觀面積分別減少3.6 km2和53.40 km2,減少幅度分別為0.38%和5.63%,減速分別為2.51%/10 a和6.8%/10 a.在2009年前后,低植被覆蓋、中度植被覆蓋面積變化劇烈,2009年后不同等級(jí)植被覆蓋的景觀面積出現(xiàn)反向變化,與2009年的氣候有直接關(guān)系.根據(jù)張紅英等人的研究[28],2009年的溫度和降水量為2006—2015年最低值,不利于植被生長(zhǎng),從而導(dǎo)致2009年低植被覆蓋和較低植被覆蓋面積增多而其他等級(jí)植被覆蓋面積減少.
圖3 2006—2015年晉城市不同植被覆蓋等級(jí)面積及比例變化Fig.3 Changes of area and proportion by vegetation coverage levels in Jincheng City during 2006-2015
3.2.1 不同高程下植被覆蓋的變化
根據(jù)自然間斷分級(jí)法結(jié)合研究區(qū)地形特征,將研究區(qū)高程分為5類:<600 m、600~1000 m、1000~1400 m、1400~1800 m、>1800 m,分析研究區(qū)不同高程下植被覆蓋指數(shù)的年際變化.
從表1中可以看出,晉城市2006—2015年不同高程下的植被覆蓋指數(shù)均呈現(xiàn)出上升的趨勢(shì),當(dāng)高程增加時(shí),植被覆蓋指數(shù)也隨之增加而且高程越高穩(wěn)定性越大.到2015年,高程小于600 m的植被覆蓋指數(shù)達(dá)到0.82,相比于2006年增長(zhǎng)8.31%;高程600~1000 m的植被覆蓋指數(shù)達(dá)到0.81,相比于2006年增長(zhǎng)6.17%;高程1000~1400 m的植被覆蓋指數(shù)達(dá)到0.94,相比于2006年增長(zhǎng)6.07%;高程1400~1800 m的植被覆蓋指數(shù)達(dá)到1,相比于2006年增長(zhǎng)2.15%.在2009年,各種高程下的植被覆蓋指數(shù)都是最低點(diǎn),主要是由于2009年晉城市降雨量偏低[29],不利于作物生長(zhǎng).總體上來說,該地區(qū)的植被覆蓋在這10年間是增加的.
表1 2006—2015年晉城市不同高程植被覆蓋指數(shù)變化Table 1 Changes of vegetation coverage index by elevations in Jincheng City during 2006-2015
2.2.2 不同坡度下植被覆蓋指數(shù)的變化
根據(jù)第二次全國(guó)土地調(diào)查耕地坡度分級(jí)體系將坡度劃分為5類:0~2°、2~6°、6~15°、15~25°、>25°,研究不同坡度下的植被覆蓋指數(shù)的年際變化.
通過分析表2可知,晉城市2006-2015年不同坡度下的植被覆蓋指數(shù)均呈現(xiàn)上升趨勢(shì),隨著坡度的增加,植被覆蓋指數(shù)也在升高,而且坡度越大穩(wěn)定性越高[30].到2015年,坡度0~2°的植被覆蓋指數(shù)達(dá)到0.75,相比于2006年增長(zhǎng)5.32%;坡度2~6°的植被覆蓋指數(shù)達(dá)到0.77,相比于2006年增長(zhǎng)5.49%;坡度6~15°的植被覆蓋指數(shù)達(dá)到0.84,相比于2006年增長(zhǎng)6.22%;坡度15~25°時(shí)植的植被覆蓋指數(shù)達(dá)到0.92,相比于2006年增長(zhǎng)6.34%;坡度大于25°的植被覆蓋指數(shù)達(dá)到0.99,相比于2006年增長(zhǎng)5.77%.在2009年時(shí),不同坡度下的植被覆蓋指數(shù)均為最小值,在2012年時(shí),不同坡度下的植被覆蓋指數(shù)均為最大值,這與每年的降水量有著密切的關(guān)系[29].
表2 2006—2015年晉城市不同坡度植被覆蓋指數(shù)變化Table 2 Changes of vegetation coverage index by slope gradients in Jincheng City during 2006-2015
3.2.3 不同地理單元下植被覆蓋指數(shù)的變化
分析植被覆蓋指數(shù)變化的另一個(gè)重要影響因子就是地理單元.通過對(duì)不同縣域2006—2015年的植被覆蓋指數(shù)進(jìn)行統(tǒng)計(jì),分析不同地理單元下晉城市的植被覆蓋指數(shù)的年際變化(表3).其下屬縣域包含沁水縣、高平市、陵川縣、陽城縣、澤州縣和城區(qū).
表3 2006—2015年晉城市不同地理單元下植被覆蓋指數(shù)變化Table 3 Changes of vegetation coverage index by geographical units in Jincheng City during 2006-2015
由表3可知,各縣市的植被覆蓋指數(shù)在2006—2015年均呈現(xiàn)增長(zhǎng)趨勢(shì),晉城市整體植被覆蓋指數(shù)增長(zhǎng)6.07%.沁水縣、高平市、陵川縣、陽城縣、城區(qū)和澤州縣植被覆蓋指數(shù)隨時(shí)間呈現(xiàn)增長(zhǎng)趨勢(shì),分別增長(zhǎng)5.77%、4.44%、5.30%、7.34%、3.86%、6.61%.除陵川縣和城區(qū)外,其他縣域在2009年的植被覆蓋指數(shù)均為最小值.在晉城市所有縣域中陵川縣整體植被覆蓋指數(shù)最好,每年植被覆蓋指數(shù)均達(dá)到0.85以上,高平市整體植被覆蓋指數(shù)最差,每年的植被覆蓋指數(shù)均在0.8以下.
斑塊數(shù)量和斑塊密度可以反映出景觀的破碎程度,其值越大,破碎度越高,值越小,破碎度越低.從圖4a、b中可以看出,2006—2015年晉城市區(qū)域的斑塊數(shù)量和斑塊密度呈現(xiàn)波浪式變化,總體為上升趨勢(shì),說明晉城市景觀異質(zhì)性增強(qiáng),破碎度提高.在2013年斑塊數(shù)量和斑塊密度達(dá)到最大值,破碎化程度最為劇烈;在2010年斑塊數(shù)量和斑塊密度均為最小值,破碎化程度相對(duì)較小.
景觀形狀指數(shù)是衡量斑塊形狀復(fù)雜程度的重要指數(shù),其取值范圍是LSI≥0,LSI值越大,表明景觀復(fù)雜度越高.由圖4c可知,2006—2015年晉城市景觀形狀指數(shù)總體上出現(xiàn)降低趨勢(shì),相比于2006年,2015年LSI值降低了0.88,說明在研究期間晉城市景觀復(fù)雜程度減弱;在2013年時(shí)LSI達(dá)到最大值,研究區(qū)景觀復(fù)雜程度最高,2015年LSI為最小值,研究區(qū)景觀復(fù)雜程度最低;從折線圖中可以看出,晉城市LSI值一直處于不斷波動(dòng)變化中,表明晉城市植被覆蓋受外界因素影響較大.
散布與并列指數(shù)用來反映不同覆蓋等級(jí)植被的混合程度,其取值在0~100之間,當(dāng)IJI值較小時(shí)說明某一斑塊類型幾乎不與其他類型相鄰,混雜程度低,反之,混雜程度高.從圖4d可以看出,除2006年和2015年外,晉城市IJL值都在50以上,植被混雜程度大,異質(zhì)性較高.2006年后IJL值開始上升并趨于穩(wěn)定,到2014年急劇下降并在2015年成為最小值點(diǎn),說明研究區(qū)植被混雜程度先升高后降低,總體上呈現(xiàn)降低趨勢(shì).
SHEI即香農(nóng)均勻性指數(shù),取值在0~1之間.當(dāng)SHEI值較小時(shí),表明景觀僅由少數(shù)斑塊類型組成,斑塊優(yōu)勢(shì)度高;SHEI接近1時(shí)表明各斑塊類型均勻分布,斑塊優(yōu)勢(shì)度低.如圖4e所示,除2006年和2015年外,晉城市香農(nóng)均勻性指數(shù)均在0.8以上,植被景觀沒有特別明顯的優(yōu)勢(shì)類型,分布較為均勻.晉城市SHEI波動(dòng)較頻繁,在2015年SHEI值降低到最小值0.51,植被景觀均勻度減少,說明高植被覆蓋類型斑塊逐漸處于優(yōu)勢(shì)地位,晉城市景觀格局越來越由少數(shù)斑塊類型所控制.
SHDI即香農(nóng)多樣性指數(shù),其取值大于零,SHDI值越小,全部景觀中所含的斑塊類型越少,SHDI值越大,說明斑塊類型數(shù)量越多,各種斑塊類型在景觀中呈現(xiàn)均衡式分布.從圖4f中可以看出,SHDI總體為下降趨勢(shì),到2015年SHDI值降低到1.34,比2006年減少0.11,表明研究期間高植被覆蓋類型斑塊對(duì)景觀的支配程度加強(qiáng),景觀多樣性下降.其中2006—2009年SHDI值增加,2009—2012年SHDI值降低,說明2006—2009年景觀多樣性在升高,斑塊分布趨于穩(wěn)定;2009—2012年景觀多樣性急劇減少,景觀多樣性降低,2009年是晉城市景觀多樣性變化的重要拐點(diǎn).
圖4 2006—2015年晉城市景觀指數(shù)變化圖Fig.4 Changes of landscape index in Jincheng City during 2006-2015
(1)低植被覆蓋
2006—2015年低植被覆蓋指數(shù)斑塊數(shù)量較少,從2006年的292上升到2015年的310,斑塊破碎度較低但有略微提升,植被發(fā)生輕微退化(圖5).低植被覆蓋景觀形狀指標(biāo)均較小并且變化波動(dòng)較小,說明低植被覆蓋斑塊邊界簡(jiǎn)單,形狀較為規(guī)則.最大斑塊所占面積比例(LPI)較小,表明低植被覆蓋在整個(gè)景觀中的優(yōu)勢(shì)度低,其中2009年出現(xiàn)異常高值,這主要是因?yàn)楫?dāng)年晉城市降水量偏低,低植被覆蓋優(yōu)勢(shì)度提高.從平均最近距離(MNN)和聚合度(AI)來看,低植被覆蓋指數(shù)斑塊之間相隔距離遠(yuǎn),分布較離散且主要由連通度較好的大斑塊組成.
圖5 2006—2015年晉城市不同植被覆蓋等級(jí)斑塊類型指數(shù)變化圖Fig.5 Variation diagrams of patch types by vegetation coverage levels in Jincheng City during 2006-2015
(2)較低植被覆蓋
2006—2015年較低植被覆蓋斑塊數(shù)量多,呈現(xiàn)上升趨勢(shì),在2011—2013年斑塊數(shù)量起伏波動(dòng)較大,斑塊破碎度高且仍在提高.從最大斑塊所占面積比例和聚合度方面觀察,LPI和AI均較小,說明較低植被覆蓋在整個(gè)景觀中并不處于優(yōu)勢(shì)地位且連通度差,并且在2010—2012年LPI和AI明顯降低,2012年后又迅速上升,可能是因?yàn)?012年外界因素不利于較低植被覆蓋作物生長(zhǎng),從而使較低植被覆蓋優(yōu)勢(shì)度降低,斑塊連通性減弱.
(3)中度植被覆蓋
2006—2015年中度植被覆蓋斑塊數(shù)量較多,呈現(xiàn)上升趨勢(shì),從2006年的442增加到2015年的601,斑塊數(shù)量增多159,斑塊破碎度較高且仍在提高.景觀形狀指數(shù)較高并且在2011年后逐漸趨于穩(wěn)定,說明被植被覆蓋斑塊邊界復(fù)雜,形狀不太規(guī)則,直到2011年景觀形狀逐漸趨于穩(wěn)定.最大斑塊所占面積比例每年波動(dòng)變化大,說明中度植被覆蓋受外界干擾較大.
(4)較高植被覆蓋
2006—2015年較高植被覆蓋斑塊數(shù)量多,出現(xiàn)略微下降,斑塊破碎度高但有輕微降低;景觀形狀指數(shù)每年都是5個(gè)植被覆蓋等級(jí)中的最大值,說明較高植被覆蓋復(fù)雜度最高,形狀最不規(guī)則.最大斑塊所占面積比例基本維持在最低點(diǎn),表明較高植被覆蓋并不占優(yōu)勢(shì).另外較高植被覆蓋平均最近距離小,兩兩斑塊之間距離較近,分布集中.
(5)高植被覆蓋
2006—2015年高植被覆蓋指數(shù)斑塊數(shù)量較少,呈現(xiàn)上升趨勢(shì),斑塊破碎度較低但正在逐漸升高.景觀形狀指數(shù)低并且總體處于穩(wěn)定狀態(tài),說明高植被覆蓋斑塊邊界簡(jiǎn)單,形狀較為規(guī)則.最大斑塊所占面積比例大,但每年的波動(dòng)起伏大,說明高植被覆蓋在整個(gè)景觀中處于優(yōu)勢(shì)地位但受外界干擾較大.從聚合角度觀察,高植被覆蓋指數(shù)AI最大并且較為穩(wěn)定,表明高覆蓋植被斑塊連通性好,面積集中.
本文以晉城市為研究對(duì)象,提取了2006—2015年的植被覆蓋指數(shù),結(jié)合研究區(qū)具體情況和前人研究成果,將其劃分為低植被覆蓋區(qū)、較低植被覆蓋區(qū)、中度植被覆蓋區(qū)、較高植被覆蓋區(qū)和高植被覆蓋區(qū).運(yùn)用空間分析和景觀生態(tài)學(xué)方法,對(duì)晉城市2006—2015年植被覆蓋指數(shù)的時(shí)空變化進(jìn)行研究,分析其植被覆蓋指數(shù)結(jié)構(gòu)、景觀格局的年際變化,得到以下結(jié)論:
(1)研究區(qū)以高植被覆蓋為主,其中高植被覆蓋指數(shù)和較高植被覆蓋指數(shù)面積在增多,而其中以高植被覆蓋指數(shù)斑塊面積增加最多;中度、較低和低植被覆蓋指數(shù)面積在減少,其中中度植被覆蓋指數(shù)面積減少最多.另外每年的高植被覆蓋指數(shù)和中度植被覆蓋指數(shù)均占了很大比重,但中度植被覆蓋面積不斷減少,向高植被覆蓋轉(zhuǎn)變.
(2)研究區(qū)10年間的植被覆蓋指數(shù)均有不同程度的增加,并且隨著高程越高,坡度越大,植被覆蓋指數(shù)也在逐漸增加.在2009年,由于降水量偏少,該地區(qū)不同高程和坡度下的植被覆蓋指數(shù)均成為了這10年間的最小值.在各個(gè)縣域,植被覆蓋指數(shù)都是增加的,其中沁水縣、陵川縣和陽城縣增加明顯,高平市、城區(qū)和澤州縣僅有輕微增加.總體來說,晉城市在10年間植被覆蓋指數(shù)是增加的,貫徹了國(guó)家的退耕還林,退耕還草政策.
(3)景觀水平上,10年間晉城市景觀破碎度增強(qiáng),景觀多樣性降低,高植被覆蓋斑塊對(duì)景觀的支配程度加強(qiáng)并逐漸處于優(yōu)勢(shì)地位,植被景觀均勻度減少,越來越由少數(shù)景觀植被類型所控制.類型水平上,低植被覆蓋指數(shù)和高植被覆蓋指數(shù)景觀破碎度程度小,形狀較為規(guī)則,高植被覆蓋在整個(gè)景觀中處于優(yōu)勢(shì)地位.
本研究在選取景觀類型水平和景觀格局水平兩個(gè)尺度上的指數(shù)進(jìn)行計(jì)算時(shí),主要是通過前人研究成果確定景觀格局分析指數(shù),因此還需進(jìn)一步結(jié)合研究區(qū)的特征選取不同的指數(shù)計(jì)算分析.此外晉城市地區(qū)的植被景觀既受人類活動(dòng)影響,又受自然因素制約.本文只考慮了地形因子,沒有考慮降雨量、濕度、溫度、日照等的影響.隨著降雨量、濕度、溫度、日照等條件的改變,植被景觀會(huì)產(chǎn)生什么樣的變化也值得深入研究.