成渝雙城經(jīng)濟(jì)圈NPP 時(shí)空變化及氣候因子驅(qū)動(dòng)分析

梁帥, 辛宇, 辜寄蓉

四川師范大學(xué)地理與資源科學(xué)學(xué)院, 四川 成都 610068

政府間氣候?qū)ｉT委員會(huì)(IPCC)第六次評(píng)估報(bào)告提出：自第一次工業(yè)革命以來，人類活動(dòng)所排放的溫室氣體是導(dǎo)致全球氣候變暖加劇的重要原因，如何實(shí)現(xiàn)碳中和和碳達(dá)峰以應(yīng)對(duì)氣候變化挑戰(zhàn)成為近年來研究熱點(diǎn)[1]。植被凈初級(jí)生產(chǎn)力(NPP)是指植被在單位面積、單位時(shí)間內(nèi)固定有機(jī)干物質(zhì)的總量[2]，其作為陸地碳循環(huán)的重要組成部分，是判定生態(tài)系統(tǒng)碳循環(huán)的重要指標(biāo)[3-4]。探究植被NPP 動(dòng)態(tài)變化及其驅(qū)動(dòng)因素，可以反映生態(tài)系統(tǒng)固碳能力的時(shí)空分異規(guī)律，揭示生態(tài)環(huán)境與氣候變化之間的關(guān)系，為生態(tài)系統(tǒng)保護(hù)及“產(chǎn)碳端”和“固碳端”的動(dòng)態(tài)平衡提供參考意見[5]。早在19 世紀(jì)80 年代就已經(jīng)有學(xué)者開始研究植被NPP，但由于技術(shù)局限性，早期大多研究都停留在實(shí)地調(diào)查再構(gòu)建回歸模型計(jì)算NPP[6-8]，這種獲取植被NPP 的傳統(tǒng)方法不僅耗時(shí)、耗力且難以開展。隨著遙感技術(shù)的發(fā)展，利用遙感數(shù)據(jù)與數(shù)學(xué)模型計(jì)算植被NPP 已經(jīng)成為大尺度區(qū)域研究NPP的重要方法[9]。這種方法能獲取高時(shí)密NPP 數(shù)據(jù)，避免了地面實(shí)測(cè)數(shù)據(jù)時(shí)間不連續(xù)、區(qū)域尺度小等缺點(diǎn)，并能夠利用數(shù)據(jù)對(duì)其進(jìn)行長(zhǎng)時(shí)間序列監(jiān)測(cè)及探索其驅(qū)動(dòng)因子，為決策者制定氣候政策達(dá)到碳達(dá)峰、碳中和等目標(biāo)提供重要參考。近年來，許多學(xué)者利用MODIS 數(shù)據(jù)分析了不同時(shí)間、區(qū)域下的植被NPP 動(dòng)態(tài)變化規(guī)律，研究發(fā)現(xiàn)近20 年間在我國(guó)植被NPP 在空間上呈現(xiàn)由西北向東南依次遞增的格局，在時(shí)間上總體呈上升趨勢(shì)[10-14]，此外研究表明NPP 的時(shí)空變化受氣候、高程以及人類活動(dòng)組合影響[15]，其中氣候因子直接影響植被覆蓋度變化對(duì)植被固碳趨勢(shì)的貢獻(xiàn)度即NPP 和GPP[16-18]。

成渝地區(qū)雙城經(jīng)濟(jì)圈作為長(zhǎng)江上游的重要城市群，其生態(tài)環(huán)境保護(hù)及可持續(xù)發(fā)展對(duì)整個(gè)長(zhǎng)江流域的社會(huì)經(jīng)濟(jì)高質(zhì)量發(fā)展有著舉足輕重的作用。在“十四五”時(shí)期，成渝地區(qū)生態(tài)文明建設(shè)進(jìn)入?yún)f(xié)同推進(jìn)減污降碳、實(shí)現(xiàn)生態(tài)環(huán)境質(zhì)量改善由量變到質(zhì)變的重要階段，監(jiān)測(cè)成渝地區(qū)雙城經(jīng)濟(jì)圈的植被NPP 時(shí)空變化及其氣候因子影響，為長(zhǎng)江上游生態(tài)文明建設(shè)及生態(tài)屏障構(gòu)建提供決策性意見?，F(xiàn)已有許多學(xué)者對(duì)川渝、西南地區(qū)的植被動(dòng)態(tài)變化做了眾多研究，例如茆楊、何宏昌等學(xué)者基于MODIS-NPP 數(shù)據(jù)對(duì)西南地區(qū)植被NPP 的時(shí)空變化及影響因素進(jìn)行了研究[10,19]，郭斌、王耠熠等學(xué)者利用MODIS 及氣象數(shù)據(jù)分析了近20 年來川西高原若爾蓋地區(qū)的NPP 時(shí)空變化特征及氣候響應(yīng)[20-21]?；贛ODIS-NPP 數(shù)據(jù)對(duì)2001—2020 年成渝地區(qū)雙城經(jīng)濟(jì)圈植被NPP 時(shí)空變化特征進(jìn)行研究，并利用ANUSPLIN 軟件插值氣象站點(diǎn)數(shù)據(jù)分析了植被NPP 對(duì)氣候變化的響應(yīng)，為長(zhǎng)江上游城市群的經(jīng)濟(jì)發(fā)展與生態(tài)環(huán)境保護(hù)的動(dòng)態(tài)平衡提供參考價(jià)值。

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

1.1 研究區(qū)概況

成渝地區(qū)雙城經(jīng)濟(jì)圈，位于“一帶一路”和長(zhǎng)江經(jīng)濟(jì)帶交匯處，以成都和重慶為中心，包含四川15 個(gè)市和重慶29 個(gè)區(qū)縣（見圖1），是西南地區(qū)經(jīng)濟(jì)水平最高，發(fā)展?jié)摿ψ畲蟮某擎?zhèn)化區(qū)域[22]。地形與四川盆地基本一致，中部多為平原，周圍由低山、丘陵等環(huán)繞，氣候適宜，自然資源豐富，生態(tài)稟賦優(yōu)良[23]。

圖1 研究區(qū)地理位置示意圖Fig. 1 Schematic diagram of geographical location of the study area

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

植被NPP 數(shù)據(jù)來源于美國(guó)NASA EOS/MODIS的2002—2021 年（實(shí)際計(jì)算植被NPP 年份為2001—2020 年）的MOD17A3HGF6.0 版本的全球NPP 數(shù)據(jù)，空間分辨率為500 m?；贕EE 平臺(tái)對(duì)MODIS數(shù)據(jù)進(jìn)行預(yù)處理，獲取統(tǒng)一坐標(biāo)的成渝地區(qū)雙城經(jīng)濟(jì)圈20 年的植被NPP 數(shù)據(jù)。對(duì)預(yù)處理后的各年份NPP 數(shù)據(jù)乘以轉(zhuǎn)換因子并剔除多年份NPP 平均值為0 的像元。

數(shù)字高程數(shù)據(jù)來源于地理空間數(shù)據(jù)云平臺(tái)，分辨率為30 m。氣象數(shù)據(jù)來自中國(guó)氣象局?jǐn)?shù)據(jù)信息中心的“中國(guó)地面氣候資料日值數(shù)據(jù)集(V3.0)”，其數(shù)據(jù)時(shí)間尺度為2001 年1 月—2020 年12 月，包括日照時(shí)數(shù)、風(fēng)速、降水、濕度、水面蒸發(fā)、日平均溫度等6 個(gè)氣象要素逐日數(shù)據(jù)?；诔捎宓貐^(qū)雙城經(jīng)濟(jì)圈內(nèi)的氣象站點(diǎn)日觀測(cè)數(shù)據(jù)，采用ANUSPLIN 軟件，以DME 為協(xié)變量插值生成空間分辨率為1 km的氣象柵格數(shù)據(jù)，再經(jīng)過重采樣生成所用的氣象數(shù)據(jù)，空間分辨率為500 m。對(duì)比其他插值方法，ANUSPLIN 插值精確度得到顯著提高[15]。

1.3 研究方法

1.3.1 植被NPP 時(shí)空演變分析方法

采用一元線性回歸模型來分析2001—2020 年成渝地區(qū)雙城經(jīng)濟(jì)圈植被NPP 整體變化趨勢(shì)及空間異質(zhì)性。

式中， θ是NPP 年際變化趨勢(shì)的斜率，n 為統(tǒng)計(jì)時(shí)段的累計(jì)年數(shù)，NPPi為第i 年植被npp 值，若 θ＞0 則表明植被NPP 呈上升趨勢(shì)，若θ＜0 則表明植被NPP呈下降趨勢(shì)，若θ=0 則表明植被NPP 幾乎無(wú)變化。

采用F 檢驗(yàn)法對(duì)植被NPP 的變化趨勢(shì)進(jìn)行顯著性分析，其公式如下：

變異系數(shù)(Cv)表示為標(biāo)準(zhǔn)差和平均值的比值，是用來描述數(shù)據(jù)離散程度的指標(biāo)[24]，采用Cv來反映近20 年植被NPP 的波動(dòng)特征，其公式如下：

式中，NPPi為第i 年植被npp 值，是20 年內(nèi)NPP 的平均值；Cv值越小，植被NPP 變化越小，反之亦然[25]。參考ArcGIS 軟件中的自然分類法[26]，將Cv分為以下五個(gè)等級(jí)：穩(wěn)定（Cv＜=0.10）,較穩(wěn)定（0.10＜Cv＜=0.20），不穩(wěn)定（0.20 ＜Cv＜=0.30），極不穩(wěn)定（Cv＞0.30）。

1.3.2 相關(guān)性分析

基于像元偏相關(guān)分析方法計(jì)算植被NPP 與氣溫、降水影響因子的偏相關(guān)系數(shù)，其公式如下：

式中，Rxy,z表示為控制變量z 之后x 與y 的偏相關(guān)系數(shù)，x 表示NPP 值，y 表示降水量，z 表示氣溫，Rxy、Rxz、Ryz分別表示NPP 與降水量、NPP 與氣溫以及氣溫與降水量的線性相關(guān)系數(shù)。采用t 檢驗(yàn)法對(duì)偏相關(guān)系數(shù)進(jìn)行顯著性檢驗(yàn)。

復(fù)相關(guān)分析法是研究幾個(gè)要素與某一要素的相關(guān)關(guān)系[27]。采用復(fù)相關(guān)分析法探究氣溫與降水對(duì)植被NPP 的共同影響作用，其公式如下：

式中，Rxy表示NPP 與降水量的線性相關(guān)系數(shù)，Rxz,y表示NPP 與氣溫的偏相關(guān)系數(shù)。采用F 檢驗(yàn)法對(duì)復(fù)相關(guān)系數(shù)進(jìn)行顯著性檢驗(yàn)。

1.3.3 驅(qū)動(dòng)力分區(qū)準(zhǔn)則

根據(jù)陳云浩等學(xué)者研究的植被變化驅(qū)動(dòng)因子分區(qū)準(zhǔn)則，將引起成渝地區(qū)雙城經(jīng)濟(jì)圈NPP 變化的驅(qū)動(dòng)力進(jìn)行分區(qū)，根據(jù)偏相關(guān)系數(shù)與復(fù)相關(guān)系數(shù)的顯著性檢驗(yàn)結(jié)果，將驅(qū)動(dòng)力分為下表五類（見表1）[27]。當(dāng)復(fù)相關(guān)顯著時(shí)：氣溫與NPP 偏相關(guān)的像元定義為氣溫驅(qū)動(dòng)；降水與NPP 偏相關(guān)的像元定義為降水驅(qū)動(dòng)，氣溫、降水分別與NPP 偏相關(guān)的像元定義為氣溫降水強(qiáng)驅(qū)動(dòng)；氣溫、降水與NPP 偏相關(guān)不顯著的像元定義為氣溫降水共同驅(qū)動(dòng)。不滿足復(fù)相關(guān)顯著的像元定義為非氣候因子驅(qū)動(dòng)。

表1 NPP 變化驅(qū)動(dòng)因素分區(qū)準(zhǔn)則Tab. 1 Zoning criteria for NPP change drivers

2 結(jié)果與分析

2.1 植被NPP 時(shí)空演化分析結(jié)果

2.1.1 植被NPP 空間分布特征

研究區(qū)內(nèi)2001—2020 年植被NPP 均值分布如圖2所示，成渝地區(qū)雙城經(jīng)濟(jì)圈植被NPP 平均值為563.14 g C·m-2·a-1，年均值分布在＜400 g C·m-2·a-1、400—500 g C·m-2·a-1、500—600 g C·m-2·a-1、600—700 g C·m-2·a-1、 700 —800 g C·m-2·a-1、 ＞800 g C·m-2·a-1的區(qū)域分別占總面積的5.49%、28.79%、34.31%、14.33%、8.95%、8.13%。從空間分布來看，研究區(qū)內(nèi)植被NPP 呈現(xiàn)四周高中間低的分布格局，其與四川盆地地形特征有著密切聯(lián)系。結(jié)合不同海拔下植被NPP 均值分布圖來看，隨著海拔的升高，植被NPP 均值也在增加，在海拔＞3 000 m 的情況下，均值開始顯著降低，但通過整體趨勢(shì)線可以看出植被NPP 呈增加趨勢(shì)。在德陽(yáng)、成都和雅安一帶海拔較高的山脈植被NPP 均值高，在海拔較低的盆地內(nèi)部的各城市中心植被NPP 均值最低，在重慶黔江區(qū)、彭水縣等低山丘陵處植被NPP 均值較高。

圖2 成渝地區(qū)雙城經(jīng)濟(jì)圈2001—2020 年植被NPP 均值分布圖及不同海拔下植被NPP 均值分布圖Fig. 2 Mean NPP distribution of vegetation in Chengdu-Chongqing economic circle in 2001-2020 and mean NPP distribution of vegetation at different altitudes

2.1.2 植被NPP 時(shí)空變化分析

2.1.2.1 植被NPP 年際變化趨勢(shì)分析

利用一元線性回歸法及變異系數(shù)分析研究區(qū)植被NPP 的時(shí)空變化特征（見圖3a 和3b），并對(duì)變化趨勢(shì)做顯著性檢驗(yàn)，將結(jié)果分為：顯著減少、減少、穩(wěn)定、增加、顯著增加五類（見表2 和圖3c）。

表2 植被NPP 變化趨勢(shì)顯著性結(jié)果統(tǒng)計(jì)Tab. 2 Statistical results of significance of vegetation NPP change trend

圖3 成渝地區(qū)雙城經(jīng)濟(jì)圈2001—2020 植被NPP 變化趨勢(shì)Fig. 3 Variation trend of vegetation NPP in Chengdu-Chongqing economic circle from 2001 to 2020

2001—2020 年植被NPP 變異系數(shù)如圖3a 所示，可以看出成渝地區(qū)雙城經(jīng)濟(jì)圈近20 年來植被NPP 較為穩(wěn)定，大部分區(qū)域變異較弱，變異較強(qiáng)的地方存在于成都市和重慶市主城區(qū)周邊以及沿長(zhǎng)江一帶。從斜率分布圖(見圖3b)及近20 年植被NPP 變化分類圖（見圖3c）來看，其變異較強(qiáng)的區(qū)域都呈現(xiàn)顯著減少趨勢(shì)；在植被NPP 均值較高的德陽(yáng)、成都和雅安一帶的山脈以及重慶市的黔江區(qū)、彭水縣等低山丘陵處變化甚微，趨于穩(wěn)定；在城市群內(nèi)部盆地及丘陵地帶植被NPP 呈現(xiàn)增加趨勢(shì)；結(jié)合表2 統(tǒng)計(jì)結(jié)果來看，其植被NPP 呈現(xiàn)顯著增加趨勢(shì)的面積最大，占研究區(qū)總面積的70.30%；呈現(xiàn)減少趨勢(shì)的面積最小，其面積占比為0.43%。從2001—2020 年植被NPP 年均值變化（見圖3b 和3d）來看，整體呈現(xiàn)增加趨勢(shì)，其中植被NPP 呈上升趨勢(shì)的面積占比達(dá)95.08%，呈下降趨勢(shì)的面積占比為4.92%；2005 年研究區(qū)植被NPP 年均值最小為471.93 g C·m-2·a-1，植被NPP 年均值最大出現(xiàn)在2020 年為644.99 g C·m-2·a-1。

2.1.2.2 不同海拔下植被NPP 動(dòng)態(tài)變化趨勢(shì)分析

將成渝地區(qū)雙城經(jīng)濟(jì)圈海拔劃分成不同梯度等級(jí)，統(tǒng)計(jì)各海拔梯度下植被NPP 的變化情況，從圖4可以看出，在研究區(qū)內(nèi)，隨著海拔的增加，植被NPP的年變化趨勢(shì)率不斷下降，在海拔小于500 m 的區(qū)域植被NPP 增長(zhǎng)趨勢(shì)最顯著，年變化趨勢(shì)率達(dá)到8.72 g C·m-2·a-1；在海拔大于3 000 m 的區(qū)域植被NPP 最穩(wěn)定，年變化趨勢(shì)率為0.95 g C·m-2·a-1；在海拔500—1 000 m、1 000—2 000 m、2 000—2 500 m及2 500—3 000 m 高程區(qū)間內(nèi)，年變化趨勢(shì)率分別為7.24 g C·m-2·a-1、4.62 g C·m-2·a-1、2.99 g C·m-2·a-1、1.94 g C·m-2·a-1。從整體來看，成渝地區(qū)雙城經(jīng)濟(jì)圈植被NPP 整體處于增長(zhǎng)趨勢(shì)，在海拔大于3 000 m的區(qū)域人類活動(dòng)少、植被覆蓋度高其植被NPP 均值較高且變化小，在海拔小于500 m 的區(qū)域，由于城市擴(kuò)張、人類活動(dòng)干擾其植被NPP 均值相對(duì)較低，近年來天然林保護(hù)、長(zhǎng)江流域生態(tài)屏障構(gòu)建等措施致使植被NPP 均值得到顯著增加。

圖4 成渝地區(qū)雙城經(jīng)濟(jì)圈各海拔梯度植被NPP 年際變化圖Fig. 4 Interannual variation of vegetation NPP at different elevations in Chengdu-Chongqing economic circle

2.2 植被NPP 對(duì)氣候因子的響應(yīng)

2.2.1 植被NPP 與氣溫、降水的相關(guān)性分析

通過一元線性趨勢(shì)分析法，可知研究區(qū)20 年里氣溫和降水都呈現(xiàn)增長(zhǎng)趨勢(shì)，年變化趨勢(shì)率分別為0.0067℃a-1、8.14 mma-1。

研究發(fā)現(xiàn)成渝地區(qū)雙城經(jīng)濟(jì)圈植被NPP 與氣溫、降水的偏相關(guān)系數(shù)范圍都在-0.83～0.87 之間，且相關(guān)性空間格局存在異質(zhì)性（見圖5）。分析表明研究區(qū)植被NPP 與氣溫呈顯著負(fù)相關(guān)、負(fù)相關(guān)的區(qū)域面積占比分別為0.39%、25.96%，呈顯著正相關(guān)、正相關(guān)的區(qū)域面積占比分別為6.22%、67.43%；與降水呈顯著負(fù)相關(guān)、負(fù)相關(guān)的區(qū)域面積占比分別為0.21%、7.33%，呈顯著正相關(guān)、正相關(guān)的區(qū)域面積占比分別為23.40%、69.06%。

圖5 成渝地區(qū)雙城經(jīng)濟(jì)圈植被NPP 與氣溫、降水的偏相關(guān)系數(shù)（a、b）及顯著性（c、d）Fig. 5 Partial correlation coefficients (a, b) and significance (c, d) of vegetation NPP with air temperature and precipitation in Chengdu-Chongqing economic circle

從空間上來看，研究區(qū)植被NPP 與氣溫正負(fù)關(guān)系并存，在其西南部即眉山市、樂山市和雅安市等地植被NPP 與氣溫呈負(fù)相關(guān)關(guān)系，在重慶市黔江區(qū)及彭水縣部分地區(qū)存在負(fù)相關(guān)關(guān)系且零星出現(xiàn)顯著負(fù)相關(guān)關(guān)系；在研究區(qū)中部及東北部植被NPP 與氣溫正相關(guān)關(guān)系明顯，尤其是在重慶城市群及東部城市群的部分地區(qū)存在顯著正相關(guān)關(guān)系。研究區(qū)大部分區(qū)域植被NPP 與降水之間都呈現(xiàn)正相關(guān)關(guān)系，并且在成都城市群與東部城市群的區(qū)域正相關(guān)關(guān)系顯著；只有在東部城市群零星存在負(fù)相關(guān)關(guān)系。其結(jié)果表明了氣溫在不同區(qū)域?qū)χ脖籒PP 變化的影響不同，而降水在研究區(qū)內(nèi)幾乎都起著積極影響。

2.2.2 植被NPP 變化氣候因子驅(qū)動(dòng)分析

植被NPP 的變化除了受氣溫、降水等單方面因子的作用，同時(shí)受氣溫與降水等因子的共同驅(qū)動(dòng)作用，通過植被NPP 與氣候因子的相關(guān)性及其顯著性分析，將驅(qū)動(dòng)因子進(jìn)行分區(qū)，從圖6 可知，植被NPP與氣溫、降水驅(qū)動(dòng)因子的復(fù)相關(guān)系數(shù)范圍為0.002～0.867，其高值區(qū)主要分布在東部城市群及成都市與德陽(yáng)市的西部，在東部城市群的東北部、雅安市與宜賓市的屏山區(qū)系數(shù)值較低。研究發(fā)現(xiàn)成渝地區(qū)雙城經(jīng)濟(jì)圈植被NPP 變化主要受非氣候因子驅(qū)動(dòng)，其面積占比達(dá)86.87%；受氣候因子影響的區(qū)域面積占研究區(qū)總面積的13.13%，其中受氣溫驅(qū)動(dòng)、降水驅(qū)動(dòng)、氣溫降水共同驅(qū)動(dòng)、氣溫降水強(qiáng)驅(qū)動(dòng)的面積占比分別為1.58%、8.40%、0.16%、2.99%。從空間上來看，受氣溫驅(qū)動(dòng)的區(qū)域主要分布在東部城市群的平行嶺谷地帶及重慶城市群的江津區(qū)、綦江區(qū)；受降水驅(qū)動(dòng)的區(qū)域分布最廣，包括了成都城市群的西北部及南部部分區(qū)域、南部城市群的瀘州市部分區(qū)域、重慶市黔江區(qū)部分區(qū)域；受氣溫、降水共同驅(qū)動(dòng)的區(qū)域十分有限，分布稀少；在南充儀隴、蓬安縣一帶受氣溫、降水強(qiáng)驅(qū)動(dòng)的區(qū)域最為密集。

圖6 復(fù)相關(guān)系數(shù)及驅(qū)動(dòng)力分區(qū)圖Fig. 6 Partition diagram of complex correlation coefficient and driving forces

3 結(jié)論與討論

基于MOD17A 產(chǎn)品NPP 數(shù)據(jù)，利用一元線性趨勢(shì)法分析2001—2020 年成渝地區(qū)雙城經(jīng)濟(jì)圈植被NPP 變化特征，并結(jié)合氣象數(shù)據(jù)，運(yùn)用相關(guān)分析法，研究植被NPP 變化的氣候因子驅(qū)動(dòng)力。

研究結(jié)果表明成渝地區(qū)雙城經(jīng)濟(jì)圈植被近20 年間NPP 在471.93—644.99 g C·m-2·a-1范圍之間，均值為563.14g C·m-2·a-1，比全國(guó)平均值514.48 g C·m-2·a-1高48.66g C·m-2·a-1，其總體呈現(xiàn)增長(zhǎng)趨勢(shì)，與馮婉[12]、張鳳英[28]等學(xué)者的關(guān)于成渝城市群地區(qū)的研究結(jié)果一致。從年際變化上來看，研究區(qū)內(nèi)植被NPP 整體呈現(xiàn)緩慢增長(zhǎng)的趨勢(shì)，增長(zhǎng)率為7.53 g C·m-2·a-1。在植被NPP 均值較高的德陽(yáng)、成都和雅安一帶的山脈以及重慶市的黔江區(qū)、彭水縣等低山丘陵處趨于穩(wěn)定；在城市群內(nèi)部盆地及丘陵地帶植被NPP 呈現(xiàn)增加趨勢(shì)，這一現(xiàn)象得力于“退耕還林、還草政策”、“長(zhǎng)防林工程”、“長(zhǎng)江保護(hù)修復(fù)攻堅(jiān)戰(zhàn)行動(dòng)計(jì)劃”[29]。從空間上來看，研究區(qū)內(nèi)植被NPP 呈現(xiàn)四周高中間低的分布格局，隨著海拔的升高，植被NPP 值不斷增加。海拔在2 500—3 000 m 的高程區(qū)間內(nèi)植被NPP均值為805.21 g C·m-2·a-1，海拔小于500 m 的高程區(qū)間內(nèi)植被NPP 均值為491.06 g C·m-2·a-1。這與四川盆地地形息息相關(guān)，這可能是因?yàn)檠芯繀^(qū)靠近川西高原的部分高海拔區(qū)域，由于氣候影響，生態(tài)環(huán)境脆弱，植被覆蓋率低，NPP 低且較為穩(wěn)定，在盆地城市內(nèi)部由于地勢(shì)平坦，水熱條件好，其城市不斷擴(kuò)張，NPP 值較低，在低山丘陵地區(qū)，因氣候適宜而少受人類活動(dòng)影響，植被覆蓋廣，NPP 值最高，該結(jié)果與胥曉[30]、張銳[31]等學(xué)者的研究結(jié)論相符合。

此外，基于氣象站點(diǎn)數(shù)據(jù)，利用ANUSPLIN 軟件對(duì)氣象數(shù)據(jù)進(jìn)行插值，分析氣候因子對(duì)植被NPP變化的影響。結(jié)果表明研究區(qū)植被NPP 與氣溫、降水的相關(guān)性存在空間異質(zhì)性，氣溫在眉山市、樂山市和雅安市等地以及重慶市黔江區(qū)、彭水縣部分地區(qū)對(duì)植被NPP 起著抑制作用，在重慶城市群及東部城市群的部分地區(qū)促進(jìn)了植被NPP 的增長(zhǎng)，而降水與植被NPP 在整個(gè)研究區(qū)內(nèi)都呈正相關(guān)的關(guān)系，這可能是因?yàn)樗疅釛l件決定了植被的生長(zhǎng)，而降水的升高使植被獲得更多的水分，促進(jìn)植被生長(zhǎng)，使得NPP 不斷增加，但研究表明氣溫對(duì)植被NPP 的影響是個(gè)復(fù)雜的過程，氣溫升高能夠?yàn)橹脖惶峁└嗟臒崃看龠M(jìn)植被光合作用從而導(dǎo)致NPP 增加，但過高的溫帶增長(zhǎng)會(huì)對(duì)植被產(chǎn)生高溫脅迫[10,32]。此外，研究表明研究區(qū)內(nèi)受到非氣候因子影響的區(qū)域面積占比達(dá)到86.87%，植被NPP 分異較大程度受到了人為因素干擾，但僅針對(duì)氣候因子來說，氣溫、降水仍然是影響植被NPP 變化的主要因素。

以往大多數(shù)研究多選擇長(zhǎng)江流域、川西高原或橫斷山區(qū)這樣的生態(tài)區(qū)內(nèi)[4,12,33]，但對(duì)成渝這樣的城市群體研究甚少，成渝城市群作為長(zhǎng)江上游最大的城市群體，其植被NPP 變化趨勢(shì)對(duì)長(zhǎng)江上游生態(tài)屏障構(gòu)建及生態(tài)修復(fù)有著舉足輕重的影響。