土地利用變化對三峽庫區(qū)重慶段植被凈初級生產(chǎn)力的影響

趙 曉,周文佐,田 羅,何萬華,章金城,劉東紅,楊 帆

西南大學(xué)地理科學(xué)學(xué)院, 重慶 400715

植被凈初級生產(chǎn)力(Net Primary Productivity,NPP)指植物在單位時(shí)間單位面積上由光合作用所固定的有機(jī)物總量減去自養(yǎng)呼吸后的剩余部分[1]。NPP反映植被在生態(tài)系統(tǒng)中的生產(chǎn)能力,是判定生態(tài)系統(tǒng)碳積累和調(diào)節(jié)生態(tài)過程的主要因子,在碳平衡和全球變化中具有重要作用[2]。土地利用變化將直接影響植被NPP,而NPP的變化反映了植被對土地利用變化的響應(yīng)。Defries等[3]研究發(fā)現(xiàn)全球每年因土地利用變化減少的NPP約占5%；Imhoff等[4]研究表明美國由城市化引起的土地轉(zhuǎn)換每年造成0.04 PgC的碳損失；Milesi等[5]研究發(fā)現(xiàn)美國東南部由城市化引起的土地利用變化導(dǎo)致NPP每年下降0.4%；孫政國等[6]研究表明城市化使土地利用和覆蓋方式改變,導(dǎo)致NPP大量損失；Xu等[7]研究表明三峽工程使移民安置區(qū)的NPP在2000—2010年減少了8%,主要驅(qū)動(dòng)機(jī)制為土地利用變化。三峽庫區(qū)地理位置特殊,自然災(zāi)害頻發(fā),加之三峽水利工程建設(shè)的影響,導(dǎo)致庫區(qū)土地利用格局和生態(tài)環(huán)境發(fā)生劇烈變化,影響了庫區(qū)生態(tài)系統(tǒng)的結(jié)構(gòu)和功能[8]。研究庫區(qū)土地利用變化對植被NPP的影響有利于合理優(yōu)化庫區(qū)土地利用格局,對探求更多的碳增長空間、調(diào)節(jié)碳循環(huán)具有重要意義[9]。

植被NPP的估算主要有兩種方式：一是實(shí)地調(diào)查法；二是模型估算法。其中估算模型主要分為4類：氣候模型、遙感模型、過程模型和遙感-過程耦合模型[10]。本文使用的MOD17A3數(shù)據(jù)利用改進(jìn)的BIOME-BGC模型[11]與光能利用率模型進(jìn)行耦合估算全球陸地生態(tài)系統(tǒng)年NPP,屬于遙感—過程耦合模型。該數(shù)據(jù)已在全球和區(qū)域的植被生長狀況、生物量估算、全球變化等研究中得到廣泛應(yīng)用和驗(yàn)證[12- 13]。國內(nèi)學(xué)者也在東北地區(qū)、陜西省、黃河三角洲、南嶺山地等不同地形、不同大小的區(qū)域驗(yàn)證了在我國使用MOD17A3數(shù)據(jù)的適用性[1,14- 16]。

已有研究主要從土地利用的面積變化、類型轉(zhuǎn)變、時(shí)空特征、驅(qū)動(dòng)因素等方面來分析三峽庫區(qū)土地利用變化[17- 19],而對庫區(qū)土地利用變化對NPP的影響研究較少。研究土地利用變化對NPP的影響多從土地利用面積變化、土地轉(zhuǎn)移、影響因素等角度分析[9,20],較少從景觀格局的角度去探究。因此,本文以三峽庫區(qū)重慶段為研究區(qū),從景觀生態(tài)學(xué)的角度分析土地利用變化與區(qū)域NPP變化的關(guān)系,量化土地利用變化對區(qū)域NPP變化的貢獻(xiàn),并從景觀格局的角度分析影響區(qū)域NPP變化的因素,以期為優(yōu)化區(qū)域土地利用結(jié)構(gòu),提高區(qū)域固碳能力提供有益參考。

1 研究區(qū)概況

三峽庫區(qū)重慶段(28°28′—31°44′N,105°49′—110°12′E)位于我國西南部,長江上游的下段,海拔高度在148—2793 m之間,地形以山地、丘陵為主。研究區(qū)西起江津區(qū),東抵巫山縣,南起武隆縣,北至開縣,包括了重慶主城區(qū)、涪陵區(qū)、長壽區(qū)、豐都縣等共22個(gè)區(qū)縣,面積約4.6×104km2,占三峽庫區(qū)總面積的86%。研究區(qū)處于亞熱帶季風(fēng)濕潤氣候區(qū),年均溫15—18℃,具有春季較短,夏季炎熱,秋季多雨,冬季溫暖,全年降水豐富,濕度大等氣候特征[21]。主要地帶性植被為亞熱帶常綠闊葉林,主要植被類型包括農(nóng)作物、林地和草地。地勢較為平坦的地區(qū)主要為草地和農(nóng)作物,其中農(nóng)作物類型主要為水稻、玉米、小麥；海拔較高的地區(qū)如巫溪、巫山縣則以林地為主,森林植被主要為馬尾松和櫟類[22]。

2 數(shù)據(jù)來源與研究方法

2.1 數(shù)據(jù)來源及預(yù)處理

本文數(shù)據(jù)主要包括遙感數(shù)據(jù)和土地利用數(shù)據(jù)。遙感數(shù)據(jù)為EOS/MODIS的MOD17A3數(shù)據(jù),來源于美國國家地質(zhì)調(diào)查局MODIS數(shù)據(jù)分發(fā)中心(https://modis.gsfc.nasa.gov)。該數(shù)據(jù)是基于BIOME-BGC模型計(jì)算出全球陸地植被凈初級生產(chǎn)力年際變化的資料,像元值即為年NPP值,表示全年每日光合作用產(chǎn)生的有機(jī)物總量的加和減去全年自養(yǎng)呼吸后剩余的部分[15],空間分辨率為1 km。研究區(qū)涉及的數(shù)據(jù)行列號為h27v05和h27v06,2000—2015年共32景。采用MRT(MODIS Reprojection Tool)工具對數(shù)據(jù)進(jìn)行投影轉(zhuǎn)換、拼接等預(yù)處理,最后將無效值剔除[23]。土地利用數(shù)據(jù)來源于中國科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心(http://www.resdc.cn),包括2000、2005、2010、2015年共4期數(shù)據(jù),空間分辨率為1 km。該土地利用數(shù)據(jù)的分辨率較低,研究區(qū)內(nèi)較小的綠地?zé)o法識別,因此本文主要考慮研究區(qū)內(nèi)的主要植被類型,將非建成區(qū)的土地利用類型劃分為耕地、林地和草地,不考慮未利用地。

2.2 研究方法

2.2.1 NPP時(shí)空特征分析

簡單差值法：利用同一地區(qū)不同時(shí)相的圖像間的差值來表示NPP變化,公式[14]為：

一元線性回歸法：分別對2000—2005年、2005—2010年、2010—2015年、2000—2015年進(jìn)行NPP年際變化趨勢分析,其變化趨勢用多年柵格數(shù)值回歸方程的斜率(slope)[24]來反映,斜率使用最小二乘法獲得,公式為：

式中,n為研究時(shí)段年數(shù),i為年序號,VNPP,i為第i年的NPP值。slope>0表示NPP年際變化趨勢增加；反之則降低。目前對slope的劃分沒有統(tǒng)一標(biāo)準(zhǔn)[25],本文采用等間距劃分法,將NPP變化分為5個(gè)等級：明顯增長(slope≥15)、輕度增長(5≤slope<15)、穩(wěn)定不變(-5≤slope<5)、輕度降低(-15≤slope<-5)和明顯降低(slope≤-15)。

2.2.2 景觀格局變化對NPP的影響

景觀格局變化采用土地利用轉(zhuǎn)移矩陣[26]來描述。本文通過構(gòu)建2000—2015年三峽庫區(qū)重慶段的土地利用轉(zhuǎn)移矩陣,分析16年來研究區(qū)景觀格局的變化。

景觀格局變化對區(qū)域NPP的影響可以通過計(jì)算土地利用變化的面積對NPP總量變化的貢獻(xiàn)率(R)來表示[27]：

式中,ΔS和ΔNPP分別表示土地利用面積變化量和NPP變化量,S1和NPP1分別表示某一時(shí)段初始時(shí)間的土地利用面積和單位面積NPP值；ΔS×NPP1表示土地利用面積變化對NPP總量的影響量；S1×ΔNPP表示氣候變化對NPP總量的影響量；ΔNPP×ΔS表示氣候變化和土地利用共同作用對NPP總量的影響。因本文只探討土地利用對NPP變化的影響,故只計(jì)算土地利用變化對NPP的貢獻(xiàn)量。

本文基于以往的研究[18,28- 29],從景觀水平上選取了3個(gè)景觀格局指數(shù)來表征景觀格局的變化：表征景觀異質(zhì)性的香農(nóng)多樣性指數(shù)(SHannon′s Diversity Index, SHDI)、表征景觀破碎度的斑塊密度指數(shù)(Patch Density, PD)以及表征景觀聚集度的聚合度指數(shù)(Aggregation Index, AI)。本文采用移動(dòng)窗口法(20 km)通過Fragstats 4.2軟件計(jì)算2000、2005、2010、2015年的景觀格局指數(shù),再通過SPSS計(jì)算各年份各景觀格局指數(shù)與NPP的相關(guān)系數(shù),并利用主成分分析法提取3個(gè)景觀格局指數(shù)的主成分,分析主成分與NPP間的關(guān)系。

3 結(jié)果

3.1 MOD17A3數(shù)據(jù)驗(yàn)證

因?qū)崪y數(shù)據(jù)缺乏,本文采用與文獻(xiàn)值對比的方式對MOD17A3數(shù)據(jù)進(jìn)行驗(yàn)證。本文數(shù)據(jù)模擬值為590.94 gC m-2a-1,略高于其他模型模擬值,但仍處于文獻(xiàn)值324—626 gC m-2a-1范圍內(nèi)(表1)。該數(shù)據(jù)模擬值偏高可能與不同模型參數(shù)不同以及數(shù)據(jù)分辨率高低有關(guān)[35]。陳雅如等[36]使用BIOME-BGC模型模擬三峽庫區(qū)1981—2014年植被NPP,其平均值范圍：270.83—552.16 gC m-2a-1,通過3種方式驗(yàn)證表明該模型能夠真實(shí)反映三峽庫區(qū)的情況。考慮到地形因素,劉玉安等[37]對太陽總輻射和氣溫進(jìn)行地形校正,構(gòu)建改進(jìn)的CASABTC(CASA based on topographical correction)模型模擬了大別山區(qū)植被NPP,并與MOD17A3數(shù)據(jù)、實(shí)測值和文獻(xiàn)值比較,驗(yàn)證了模型在山區(qū)的適用性。CASABTC模擬值、MOD17A3數(shù)據(jù)和實(shí)測值三者較為一致,從側(cè)面說明了MOD17A3數(shù)據(jù)與實(shí)測值誤差較小,表明MOD17A3數(shù)據(jù)在山區(qū)應(yīng)用的可行性。基于以上幾點(diǎn)認(rèn)為MOD17A3數(shù)據(jù)在本研究區(qū)應(yīng)用具有可行性。

表1 MOD17A3數(shù)據(jù)與其他模型模擬值對比

CASA：光能利用率模型Carnegie-Ames-Stanford Approach；GLOPEM-CEVSA：遙感-過程耦合模型AVHRR Global Production Efficiency Model-Carbon Exchange between Vegetation, Soil and Atmosphere；*表示修正后的模型,**為本文數(shù)據(jù)結(jié)果

3.2 NPP時(shí)空變化特征

利用ArcGIS 10.3計(jì)算得到2000—2015年NPP均值的空間分布(圖1)。研究區(qū)多年NPP均值主要集中分布在400—600、600—800(單位：gC m-2a-1)兩個(gè)范圍區(qū)間。400—600范圍主要分布在長江以北的忠縣、長壽等耕地地區(qū)；600—800范圍主要分布在東北部的巫溪、巫山、奉節(jié)以及長江以南的石柱、武隆,江津區(qū)的四面山等海拔較高的林地地區(qū)；低值區(qū)(NPP<200和200—400)主要分布在重慶主城區(qū)附近以及長江沿岸兩側(cè)的生態(tài)消落地帶；高值區(qū)(800—1000、1000—1200以及NPP>1200)分布較為零散。16年來NPP均值總體變化范圍為76—1282 gC m-2a-1,空間特征表現(xiàn)為從東到西逐漸減少。

采用簡單差值法利用ArcGIS 10.3計(jì)算得到16年來研究區(qū)NPP的空間變化(圖2)。2015年與2000年相比,大多數(shù)地區(qū)的NPP值增加,只有少數(shù)地區(qū)的NPP值減少。其中,NPP值增加的區(qū)域約占研究區(qū)面積的86.85%。增加幅度主要在0—100間,約占研究區(qū)面積的57.44%,包括武隆、豐都、石柱等區(qū)縣；增加幅度在100—200的區(qū)域約占29.16%,主要有開縣、云陽、奉節(jié)、巫山；NPP值減少的區(qū)域主要集中在主城區(qū),其中NPP減少差值大于200的區(qū)域僅占總面積的0.38%。

2000—2015年,研究區(qū)NPP多年均值為590.94(單位：gC m-2a-1)。2000、2005、2010、2015年NPP均值分別為580.51、590.96、588.76、619.58。16年來NPP均值變化趨勢表現(xiàn)為：2001年發(fā)生伏旱[22],NPP值開始減少,到2002年增加到峰值668.02,2003年庫區(qū)開始蓄水,研究區(qū)生態(tài)環(huán)境得到改善[35],植被生長良好,NPP值趨于穩(wěn)定水平(NPP值在590.94上下波動(dòng)),到2011年發(fā)生極端高溫事件,導(dǎo)致NPP達(dá)到最低值464.64,2012年后又開始增加,至2014年后恢復(fù)至穩(wěn)定水平(圖3)。從整體來看,16年間NPP均值變化有增有減,其總體趨勢基本保持不變。2000—2005年,NPP變化呈輕度增長,面積約占45%,穩(wěn)定不變的區(qū)域約占29%,明顯增長的地區(qū)(9%)主要分布在研究區(qū)的中部(忠縣、豐都)；2005—2010年,NPP變化表現(xiàn)為穩(wěn)定不變的區(qū)域面積占比最大,約為40%,其次為輕度降低(26%),明顯降低的區(qū)域(5%)主要集中在研究區(qū)的中部沿河流地段和主城區(qū)；2010—2015年,NPP變化表現(xiàn)為明顯增長,約占65%,輕度增長面積占28%,主要分布在東北部地區(qū)；2000—2015年,絕大部分地區(qū)NPP值保持穩(wěn)定,增長區(qū)主要集中在研究區(qū)東北部,減少區(qū)主要集中在主城區(qū)附近(圖3、圖4)。

圖1 2000—2015年三峽庫區(qū)重慶段NPP均值分布圖 Fig.1 Distribution of the average annual NPP from 2000 to 2015 in the Chongqing TGRATGRA：三峽庫區(qū) Three Gorges Reservoir Area

圖2 三峽庫區(qū)重慶段2000—2015年NPP差值分布圖Fig.2 Difference of NPP from 2000 to 2015 in the Chongqing TGRA

圖3 2000—2015年三峽庫區(qū)重慶段NPP變化趨勢Fig.3 NPP variation trend in the Chongqing TGRA from 2000 to 2015

圖4 2000—2015年三峽庫區(qū)重慶段NPP變化斜率空間分布圖Fig.4 Spatial distribution of NPP variation slope in Chongqing TGRA from 2000 to 2015

利用分辨率為1 km的4期土地利用數(shù)據(jù),使用ArcGIS 10.3中的分區(qū)統(tǒng)計(jì)模塊(Zonal Statistics as Table)分別統(tǒng)計(jì)2000、2005、2010、2015年研究區(qū)NPP總量,得到NPP總量分別約為25.6 TgC(1 TgC=1012gC)、26.2 TgC、26.6 TgC、28.5 TgC,呈總體增長趨勢。其中耕地對NPP總量的貢獻(xiàn)最大,約為11.8 TgC,占NPP總量的44%；其次是林地,NPP總量約為10.7 TgC,占NPP總量的40%；草地約占NPP總量的14%；其他用地約占NPP總量的2%(圖5)。

圖5 2000、2005、2010、2015年三峽庫區(qū)重慶段NPP總量變化Fig.5 Total NPP variation of the Chongqing TGRA in 2000,2005,2010,2015

3.3 區(qū)域景觀格局變化

研究區(qū)景觀類型以耕地和林地為主,耕地主要分布在低海拔的西南部,林地主要分布在高海拔的東北部。2000—2015年,其他用地面積增加826 km2,3個(gè)時(shí)段變化率分別為：12.6%、30.1%、45.6%,增速不斷加快；耕地不斷減小,減小面積約970 km2,變化率分別約為：-1.2%、-1.5%、-2.4%,減速逐漸加快；林地表現(xiàn)為先增長后減少,總體呈增長趨勢,增長面積約681 km2,變化率分別為3.8%、0.4%、-0.2%,增勢逐漸減緩；草地面積也在不斷減少,16年減少了156 km2,但其變化率(-2.0%、-0.2%、-0.3%)在逐漸降低。

16年來耕地面積不斷減小,主要轉(zhuǎn)換為其他用地(777 km2)、林地(289 km2)和草地(118 km2)；林地面積增加,主要來源于耕地和草地(207 km2)的轉(zhuǎn)換；草地面積減少,除了轉(zhuǎn)換為林地外,還轉(zhuǎn)換為耕地(83 km2)和其他用地(53 km2)(表2)。

表2 2000—2015年三峽庫區(qū)重慶段土地利用轉(zhuǎn)移矩陣/km2

TGRA：三峽庫區(qū) Three Gorges Reservoir Area

2000—2015年,研究區(qū)景觀格局指數(shù)總體變化不大。SHDI各年份平均值為2.28,總體呈上升趨勢,表明斑塊類型增加；PD均值為0.32,基本保持穩(wěn)定不變,表明景觀破碎化程度變化不大；AI平均值約為66.66,AI減小,表明同一土地利用類型集中分布度降低(表3)。3個(gè)指數(shù)反映了研究區(qū)景觀格局的變化：景觀多樣性增加,土地利用類型集中度降低,景觀破碎化程度較低且維持穩(wěn)定狀態(tài)。

表3 2000、2005、2010、2015年三峽庫區(qū)重慶段景觀格局指數(shù)

SHDI：香農(nóng)多樣性指數(shù) Shannon′s Diversity Index；PD：斑塊密度 Patch Density；AI：聚合度 Aggregation Index

3.4 景觀格局變化對NPP變化的影響分析

3.4.1 土地利用變化對NPP的影響

由于其他用地(主要為城市建成區(qū))中的植被較少,盡管該類型面積變化較大,但對NPP的影響較小(圖5),因此本文只探討耕地、林地和草地對NPP變化的影響。2000—2005年,林地變化最大,增長面積約650 km2,NPP相應(yīng)增長了約0.02 TgC,增長幅度(0.02%)不大；草地面積減小,NPP相應(yīng)也減少；2005—2010年,林地持續(xù)增長,NPP相應(yīng)增加,增幅(5.77%)較大；草地面積繼續(xù)減少,但NPP總量增加；2010—2015年,林地面積開始減少,減少幅度較小,NPP仍然保持增加(5.32%)趨勢；草地面積持續(xù)減少,其NPP總量繼續(xù)增加。2000—2015年,耕地面積持續(xù)減小,其NPP表現(xiàn)為先增加后減少再增加(圖6)。

圖6 土地利用面積變化對NPP的影響Fig.6 Effect of land-use changes on NPP

根據(jù)貢獻(xiàn)率(R)計(jì)算出3個(gè)時(shí)段土地利用面積變化對NPP的影響量和貢獻(xiàn)率(表4)。2000—2005年,土地利用變化對NPP的增長表現(xiàn)為促進(jìn)作用,面積變化使NPP增加了0.17 TgC,貢獻(xiàn)率為26.49%；2005—2010年,土地利用變化抑制了NPP增長,使NPP減少了0.15 TgC,貢獻(xiàn)率為59.76%；2010—2015年,土地利用變化對NPP也表現(xiàn)為抑制作用,NPP總量減少了0.29 TgC,貢獻(xiàn)率為17.27%。

表4土利用變化對NPP總量變化的貢獻(xiàn)率

Table4Contributionrateofland-usechangestototalNPPvariation

時(shí)間段Period影響量Influence amount/(TgC)貢獻(xiàn)率(R)Contribution rate2000—20050.1726.492005—2010-0.1559.762010—2015-0.2917.27

圖7 景觀格局指數(shù)與NPP相關(guān)性分析Fig.7 Correlation between Landscape indexes and NPP SHDI：香農(nóng)多樣性指數(shù) SHannon′s Diversity Index；PD：斑塊密度 Patch Density；AI：聚合度 Aggregation Index

3.4.2 景觀格局指數(shù)變化與NPP的相關(guān)性分析

將景觀格局指數(shù)分布圖和NPP分布圖按20 km×20 km的大小劃分成連續(xù)的格網(wǎng)作為樣方[38],以該尺度分析景觀格局指數(shù)與NPP的相關(guān)性。剔除無效值和不對應(yīng)的值后共計(jì)102個(gè)格網(wǎng)。2000—2015年,SHDI、PD與NPP呈正相關(guān),SHDI各年份相關(guān)系數(shù)為0.228、0.216、0.205、0.198,PD各年份相關(guān)系數(shù)分別為0.367、0.210、0.231、0.198。而AI與NPP呈負(fù)相關(guān),各年份相關(guān)系數(shù)分別為-0.386、-0.343、-0.370、-0.321(圖7)。

分別對4個(gè)年份的3個(gè)景觀格局指數(shù)進(jìn)行主成分分析,方差大于1的成分只有一個(gè),提取該主成分。第一主成分中各年份的特征值分別約為2.368、2.332、2.375、2.282,累積載荷分別約為76%、79%、78%、79%。4個(gè)年份中第一主成分與NPP均呈顯著正相關(guān)(顯著性水平為0.01)。將各年份的第一主成分與NPP值進(jìn)行線性擬合(圖8),可以看到四個(gè)年份中的NPP值均隨第一主成分的增大而增大。

圖8 景觀指數(shù)第一主成分與NPP擬合關(guān)系Fig.8 Relationship between the first principal component of landscape indexes and NPP**表示顯著相關(guān)；樣本數(shù)n=102

4 討論

4.1 不同景觀類型與NPP變化的關(guān)系

景觀格局變化是土地利用變化的直接反映。不同時(shí)期、不同景觀類型,土地利用變化對NPP的影響不同(圖6)。2000—2005年,城市擴(kuò)張、退耕還林等因素導(dǎo)致耕地減少,林地大量增加。耕地雖減少,但其NPP總量增長,這可能與耕地的產(chǎn)量提高有關(guān)。造林帶的形成增加了林地面積,但該時(shí)段處于退耕還林初期,人造林多為幼苗,林地雖增加,但其NPP產(chǎn)量不高。2005—2010年,耕地持續(xù)減少,其NPP總量減少。查閱該時(shí)段重慶市統(tǒng)計(jì)年鑒發(fā)現(xiàn)耕地產(chǎn)量下降。該時(shí)期發(fā)生的旱災(zāi)(2006年)、雪災(zāi)(2008年)等導(dǎo)致耕地產(chǎn)量下降,間接導(dǎo)致NPP的減少。該時(shí)段林地增長減緩,但幼齡林隨著林齡的增加,NPP產(chǎn)量增加,故林地NPP總量增加。2010—2015年,耕地繼續(xù)減少,林地也開始減少,但兩者的NPP總量繼續(xù)增加,這與林齡增長和耕地產(chǎn)量提高有關(guān)。成方妍等[39]雖指出幼齡林成長為中齡林,NPP產(chǎn)量會有所下降,但黃從德等[40]研究表明重慶地區(qū)幼齡林和中齡林碳儲量呈增長趨勢,成熟林的碳儲量呈下降趨勢。草地由于水體淹沒、城鎮(zhèn)重建等因素導(dǎo)致16年來數(shù)量減少,但受到近年來生態(tài)系統(tǒng)保護(hù)政策的影響,草地NPP總量也表現(xiàn)為增長。土地利用面積變化影響NPP的增長,但不同景觀類型受不同因素的影響,導(dǎo)致NPP變化存在差異。因此,合理規(guī)劃景觀類型比例對于提高區(qū)域NPP增長具有重要意義。

土地利用變化對NPP的貢獻(xiàn)率在3個(gè)時(shí)段上存在差異。2000—2005年土地利用變化對NPP變化表現(xiàn)為促進(jìn)作用,能夠解釋NPP增長區(qū)域,而對NPP減小區(qū)域影響??；后兩個(gè)時(shí)段土地利用變化對NPP的影響表現(xiàn)為抑制作用,但NPP總量仍持續(xù)增加,可能與這兩個(gè)時(shí)段氣候變化對NPP的影響占主導(dǎo)地位有關(guān),且影響NPP變化的還有如林齡、地形等因素,也可能使NPP總量增加[41]。

4.2 景觀格局變化對NPP變化的影響

不同空間分辨率的數(shù)據(jù)對景觀格局將產(chǎn)生粒度效應(yīng)[42]。已有研究表明在景觀水平上斑塊密度PD隨空間粒度變粗,其波動(dòng)幅度趨于平緩；聚集度AI隨時(shí)間推移呈下降趨勢,對粒度變化敏感度低,即降幅不隨粒度變粗而變化；多樣性指數(shù)SHDI隨時(shí)間推移單調(diào)上升且增幅不隨粒度變粗而變化,對于粒度變化不敏感[36,43]。本文研究中,PD基本維持不變,AI由67.23下降到66.08,SHDI由2.28增加到2.31,其變化趨勢與陳雅如、朱明研究結(jié)果相符。

16年來研究區(qū)景觀異質(zhì)性增加,可能與造林帶的形成以及質(zhì)量效應(yīng)有關(guān)。造林帶的形成豐富了植被類型,使景觀多樣性增加。質(zhì)量效應(yīng)指一個(gè)物種的某些個(gè)體從能夠成功建群的區(qū)域轉(zhuǎn)移到不適其生存的區(qū)域[44]。庫區(qū)生態(tài)消落帶因質(zhì)量效應(yīng)的存在,增加了物種的多樣性,使景觀異質(zhì)性增加。研究區(qū)景觀破碎化程度基本保持不變,與張煜星等[45]研究表明生態(tài)工程的實(shí)施發(fā)揮了顯著成效,景觀破碎化明顯緩解的結(jié)果不符。這可能與粒度效應(yīng)所產(chǎn)生的指數(shù)對斑塊變化不敏感以及造林帶的形成并不意味著景觀破碎化趨勢的變好有關(guān)。研究區(qū)景觀聚集度下降,可能是林地主要分布在高海拔的山地,而耕地分布在海拔較低的地區(qū),新形成的造林帶與天然的林地并不連續(xù),造成研究區(qū)景觀聚集度下降。

SHDI、PD與NPP呈正相關(guān),表明隨著景觀類型和景觀破碎度的增加,NPP總量增加。研究區(qū)景觀破碎度雖變化不大,但因新增加的造林帶以及質(zhì)量效應(yīng)、林齡效應(yīng)的存在,植被類型增加,幼齡林和部分中齡林使NPP總量增加；AI與NPP呈負(fù)相關(guān),而較為集中分布的天然植被主要為固碳能力逐漸下降的成熟林,導(dǎo)致NPP減少。景觀指數(shù)主成分分析中,4個(gè)年份的NPP均與第一主成分呈顯著正相關(guān),表明要提高區(qū)域固碳能力需要適當(dāng)增加區(qū)域植被類型、對森林進(jìn)行撫育更新,并控制成熟林?jǐn)?shù)量。需要指出的是本文中景觀格局指數(shù)與NPP的相關(guān)性不高,這可能與所使用的土地利用數(shù)據(jù)空間分辨率不高有關(guān)。另外,第一主成分的累積載荷值最高只達(dá)到79%,雖未達(dá)到大于85%的提取要求,但未影響景觀格局指數(shù)與NPP之間的線性關(guān)系。

5 結(jié)論

三峽庫區(qū)重慶段2000—2015年NPP均值空間分布特征表現(xiàn)為從東北部山區(qū)到西部丘陵區(qū)逐漸減少；NPP空間變化有增有減,增長區(qū)主要為研究區(qū)東北部山區(qū),減少區(qū)主要集中在重慶主城區(qū)內(nèi)；16年來NPP總量呈總體上升趨勢,共增加了2.9 TgC。土地利用的類型轉(zhuǎn)變、面積變化、耕地產(chǎn)量和林地的林齡效應(yīng)等對區(qū)域植被固碳能力、NPP空間變化產(chǎn)生影響；從景觀格局的角度來看,不同景觀類型受不同因素的影響,NPP變化存在差異。景觀格局的優(yōu)化,即斑塊類型和景觀破碎度的增加可促使NPP增長,而斑塊聚集度的增長將導(dǎo)致NPP的下降。因此,在進(jìn)行區(qū)域規(guī)劃時(shí),要考慮到不同的景觀類型植被固碳能力不同,合理的景觀類型比例才能促使區(qū)域NPP增長；另外,豐富區(qū)域植被類型,增加景觀多樣性和破碎度,重視對幼齡林的培育并控制成熟林?jǐn)?shù)量,改善區(qū)域景觀格局,也能提高區(qū)域植被固碳能力,增加NPP的增長空間。