2000—2020年喀什地區(qū)植被凈初級生產(chǎn)力時空變化及其驅(qū)動因素研究

摘要：為全面且定量地掌握喀什地區(qū)植被生態(tài)系統(tǒng)生長變化狀況，基于MODIS17A3數(shù)據(jù)集，利用趨勢分析、地統(tǒng)計分析、相關(guān)性分析等方法，研究2000—2020年喀什地區(qū)植被凈初級生產(chǎn)力（Net Primary Productivity，NPP）的時空變化特征及其驅(qū)動因素。結(jié)果表明，近20 a喀什地區(qū)植被NPP均值為

126.160 g C/（m2·a），標(biāo)準(zhǔn)差為15.801 g C/（m2·a）；20 a喀什地區(qū)NPP整體呈緩慢增長趨勢，但2010年NPP值發(fā)生突變，此次突變與降水量關(guān)系密切；喀什地區(qū)植被NPP有明顯的空間自相關(guān)性和空間異質(zhì)性，變異系數(shù)為0.050～4.344，莫蘭指數(shù)為0.910，NPP高值較為集中地分布于水源豐富的綠洲地帶。喀什地區(qū)近20 a植被生長狀況持續(xù)向好，在近20 a植被NPP時空變化中，海拔、坡度和年降水量等結(jié)構(gòu)型因素發(fā)揮決定性作用。

關(guān)鍵詞：植被；凈初級生產(chǎn)力（NPP）；時空變化；驅(qū)動因素；相關(guān)性分析；趨勢分析；地統(tǒng)計分析；喀什地區(qū)

中圖分類號：Q948.1 文獻(xiàn)標(biāo)識碼：A 文章編號：1008-9500（2024）04-0-10

DOI：10.3969/j.issn.1008-9500.2024.04.042

Research on spatiotemporal Variation and driving factors of net primary productivity of vegetation in Kashi region from 2000 to 2020

ZHU Min， Muhetaer Aimaiti， Alimu Yiming， Nuerbiya Abudukeyoumu

（College of Life and Geographic Sciences， Kashi University， Kashi 844099， China）

Abstract： In order to comprehensively and quantitatively grasp the growth and changes of vegetation ecosystem in Kashi region， based on the MODIS17A3 dataset， this paper uses trend analysis， geostatistical analysis， correlation analysis and other methods to study the spatiotemporal variation characteristics and driving factors of Net Primary Productivity （NPP） of vegetation in Kashi region from 2000 to 2020. The results show that the average NPP of vegetation in Kashi region in the past 20 years was 126.160 g C/（m2·a）， with a standard deviation of 15.801 g C/（m2·a）; in the past 20 years， the overall NPP in Kashi region showed a slow growth trend， but in 2010， there was a sudden change in NPP values， which was closely related to precipitation; the vegetation NPP in Kashi region shows significant spatial autocorrelation and heterogeneity， with a coefficient of variation of 0.050～4.344 and a Moran index of 0.910， and high NPP values are concentrated in oasis areas with abundant water sources. The vegetation growth in the Kashi region has continued to improve over the past 20 years， and structural factors such as altitude， slope， and annual precipitation have played a decisive role in the spatiotemporal variation of vegetation NPP.

Keywords： vegetation; Net Primary Productivity （NPP）; spatiotemporal variation; driving factors; correlation analysis; trend analysis; geostatistical analysis; Kashi region

凈初級生產(chǎn)力（Net Primary Productivity，NPP）是植被在單位時間、單位面積內(nèi)扣除呼吸消耗后剩余的自身固定的能量或生產(chǎn)的有機(jī)物質(zhì)，作為表征陸地生態(tài)系統(tǒng)的關(guān)鍵變量，能反映植被在自然狀態(tài)下的生產(chǎn)能力，是全球碳循環(huán)的重要組成部分[1-4]。同時，NPP也是評價生態(tài)系統(tǒng)結(jié)構(gòu)與功能協(xié)調(diào)性與生物圈人口承載力的重要指標(biāo)[5]。因此，研究區(qū)域植被NPP的時空變化及其驅(qū)動因素具有十分重要的意義。

目前，學(xué)者對NPP的研究主要集中在兩個方面。一是對NPP估算模型的研究。Potter等[6]提出CASA（Carnegie-Ames-Stanford Approach）模型，并對全球NPP進(jìn)行估算，是目前應(yīng)用最廣泛的模型之一。朱文泉等[7-9]對國外模型進(jìn)行改進(jìn)，使其更適用于我國NPP的遙感估算。二是基于地理信息技術(shù)分析NPP時空變化特征及其驅(qū)動因素。洪辛茜等[10]基于土地覆蓋和NPP產(chǎn)品分析2001—2008年西南地區(qū)NPP變化量及其驅(qū)動因素。周雷雷等[11]以胡煥庸線為界，研究中國東西部植被NPP的變化。整體來看，學(xué)者對NPP估算模型的研究較為豐富，但對NPP時空演變的研究較多集中于大尺度范圍，如西南地區(qū)、西北地區(qū)以及中國全域等，對小尺度范圍的關(guān)注甚少。但是，小尺度范圍的研究更為精細(xì)化，是區(qū)域生態(tài)環(huán)境質(zhì)量評價必不可少的內(nèi)容。此外，學(xué)者對NPP影響因素的研究大量考慮氣候因子與NPP的關(guān)系，對其他自然因子少有提及。NPP的時空分異是多個因子綜合作用的結(jié)果，并非某個因子決定，因此本研究綜合考慮多個因子對NPP的影響。

喀什地區(qū)作為絲綢之路經(jīng)濟(jì)帶的黃金地帶和中巴經(jīng)濟(jì)走廊的廊橋，是我國向西開放的重要窗口[12]。喀什地區(qū)遠(yuǎn)離海洋且三面環(huán)山，屬于典型的溫帶大陸性氣候區(qū)，降水量偏少，年均降水量介于60.1～135.6 mm[13]。喀什地區(qū)植被類型中，荒漠植被所占面積最大，針葉林面積最少，值得注意的是，喀什地區(qū)東臨塔克拉瑪干沙漠且西南部海拔過高，因此部分區(qū)域無植被覆蓋。喀什地區(qū)區(qū)位優(yōu)勢獨(dú)特，其快速、高質(zhì)量的發(fā)展對實現(xiàn)新疆維吾爾自治區(qū)（簡稱新疆）經(jīng)濟(jì)的跨越式發(fā)展有重要意義[14]，而經(jīng)濟(jì)發(fā)展與生態(tài)環(huán)境的關(guān)系尤為密切。本文基于MODIS17A3數(shù)據(jù)集，定量研究喀什地區(qū)2000—2020年的植被NPP時空變化，并選取地形、氣候等要素探討其與NPP的關(guān)系，以期為喀什地區(qū)生態(tài)環(huán)境評估與治理提供科學(xué)依據(jù)。

1 材料與方法

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

本研究使用的數(shù)據(jù)主要有喀什地區(qū)的NPP數(shù)據(jù)、高程數(shù)據(jù)、氣象數(shù)據(jù)和土地利用數(shù)據(jù)等。其中，NPP數(shù)據(jù)來自美國國家航空航天局的2000—2020年MODIS17A3數(shù)據(jù)集，時間分辨率為1 a，空間分辨率為500 m×500 m。該數(shù)據(jù)集運(yùn)用BIOME-BGC生物地球化學(xué)模型估算全球陸地植被年NPP，已在全球和區(qū)域NPP研究中廣泛應(yīng)用。利用MRT（MODIS Reprojection Tool）工具對其進(jìn)行格式轉(zhuǎn)換、投影轉(zhuǎn)換和拼接等預(yù)處理，然后在ArcGIS軟件中剔除無效值、轉(zhuǎn)換系數(shù)并變換單位，考慮研究區(qū)范圍較小和提高圖像分辨率，對圖像進(jìn)行重采樣處理，將像元大小調(diào)整為250 m×250 m。空間分辨率30 m×30 m的高程數(shù)據(jù)來源于地理空間數(shù)據(jù)云（http：//www.gscloud.cn）；氣象數(shù)據(jù)來源于國家科技資源共享服務(wù)平臺——國家地球系統(tǒng)科學(xué)數(shù)據(jù)中心（http：//www.geodata.cn）；土地利用類型數(shù)據(jù)來源于中國科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)平臺（http：//www.resdc.cn）。

1.2 NPP時空特征分析方法

1.2.1 Mann-Kendall檢驗

Mann-Kendall（MK）檢驗是一種非參數(shù)的統(tǒng)計檢驗方法，不僅可以檢驗時間序列的變化趨勢，還可以檢驗時間序列是否發(fā)生突變。檢驗結(jié)果組成UFk和UBk兩條曲線，UFk表示趨勢變化統(tǒng)計量，UBk表示序列UFk的反序列，二者在置信區(qū)間內(nèi)的交點為突變點。該方法的計算公式和參數(shù)含義見文獻(xiàn)[15-16]。

1.2.2 趨勢分析

喀什地區(qū)植被凈初級生產(chǎn)力的年際變化率反映植被凈初級生產(chǎn)力的年際變化，本研究用一元線性回歸模型對像元尺度上NPP的變化趨勢進(jìn)行模擬，如式（1）所示。若趨勢變化率大于0，則NPP總體呈上升趨勢；反之，則NPP總體呈下降趨勢。趨勢變化率的大小反映上升或下降的速率快慢。此外，采用F檢驗對植被凈初級生產(chǎn)力的趨勢變化做顯著性檢驗。

式中：l為趨勢變化率；n為研究年數(shù)，取21；i為研究年序號；yi為第i年的植被NPP。

1.2.3 空間自相關(guān)分析

空間自相關(guān)反映某個區(qū)域單元上或某個地理現(xiàn)象的屬性與相鄰區(qū)域單元或地理現(xiàn)象相同屬性的相關(guān)程度，可以度量空間域中聚集程度[17]。常用的空間自相關(guān)指數(shù)是莫蘭指數(shù)，根據(jù)其空間范圍的大小、空間自相關(guān)的測度，可分為全局莫蘭指數(shù)與局部莫蘭指數(shù)[18]。計算公式及計算方法見文獻(xiàn)[19]。

1.3 NPP影響因素分析

1.3.1 地統(tǒng)計分析

喀什地區(qū)植被NPP作為區(qū)域化變量，其空間異質(zhì)性可以通過地統(tǒng)計學(xué)來分析。通過構(gòu)建半變異函數(shù)模型對區(qū)域化變量的空間變異結(jié)構(gòu)進(jìn)行研究[20]，半變異函數(shù)模型常用來研究區(qū)域化變量隨機(jī)性和空間相關(guān)性的規(guī)律，主要參數(shù)有基臺值、塊金值和變程[21]，如式（2）所示。

式中：γ（h）為半變異函數(shù)；h為步長；N（h）為間隔距離h的樣本總配對數(shù)；Z（xi）和Z（xi+h）分別為在位置xi和xi+h處的NPP值；i為樣本配對序號，i=1，2，…，N（h）。

1.3.2 相關(guān)性分析

相關(guān)系數(shù)用于表征兩個或多個變量的相關(guān)程度。本研究在像元尺度上分析喀什地區(qū)植被NPP與影響因素的相關(guān)系數(shù)[22]，如式（3）所示。相關(guān)系數(shù)的取值范圍為[-1，1]。若相關(guān)系數(shù)大于0，則兩個變量正相關(guān)；若相關(guān)系數(shù)小于0，則兩個變量負(fù)相關(guān)；若相關(guān)系數(shù)為0，則兩個變量無線性關(guān)系。相關(guān)系數(shù)的絕對值越接近1，表明相關(guān)性越好。

式中：rxy為相關(guān)系數(shù)；xi為第i年的NPP值；x為多年NPP均值；yi為第i年的影響因素；y為多年影響因素的均值；n為研究年數(shù)，取21。

2 結(jié)果與分析

2.1 NPP時間序列分析

NPP（Y）與年份（X1）的線性擬合分析結(jié)果如圖1所示，二者的相關(guān)系數(shù)R2為0.461。喀什地區(qū)近20 a NPP年均值呈現(xiàn)上升趨勢，均值在101～153 g C/（m2·a）上下浮動。最低值出現(xiàn)在2000年和2007年，最高值出現(xiàn)在2010年，且2010年以后，NPP年均值基本都高于20 a的NPP年平均值，20 a的NPP年平均值為126.160 g C/（m2·a）。喀什地區(qū)近20 a NPP變異系數(shù)[23]為12.525%，表明近20 a NPP的整體波動性不高。

以5年為一個時間周期，分析喀什地區(qū)NPP的相關(guān)變化。由表1可知，NPP年均值持續(xù)增長，2011—2015年NPP年均值增長陡然加快，從118.200 g C/（m2·a）猛增至132.400 g C/（m2·a）。喀什地區(qū)NPP的標(biāo)準(zhǔn)差變化顯示，除2006—2010年以外，其余時間段的NPP標(biāo)準(zhǔn)差都較小，說明這些年份的NPP值與其各自時間周期內(nèi)的NPP均值相差不大。變異系數(shù)也表明2006—2010年的NPP值波動相對較大，變異系數(shù)為18.020%；2016—2020年NPP的波動最小，僅為5.253%。

借助Python軟件對喀什地區(qū)2000—2020年植被NPP變化趨勢及突變情況進(jìn)行分析。Mann-Kendall趨勢檢驗的結(jié)果顯示，相關(guān)性系數(shù)為0.002，近20 a喀什地區(qū)NPP呈增長趨勢，這與前面統(tǒng)計分析的結(jié)果一致。喀什地區(qū)NPP突變檢驗結(jié)果如圖2所示，UFk與UBk兩條曲線出現(xiàn)交點，且交點位于兩條臨界線之間，交點對應(yīng)的時間便是NPP發(fā)生突變的時間，即2010年。2010年，NPP出現(xiàn)較大增長，達(dá)到近20 a的最高值。

2.2 NPP空間分異分析

對喀什地區(qū)近20 a植被凈初級生產(chǎn)力進(jìn)行疊加分析，結(jié)果發(fā)現(xiàn)，喀什地區(qū)NPP在空間分布上具有明顯的差異性，如圖3所示。喀什地區(qū)2000—2020年

植被NPP的空間分布特征顯示，高值較為集中地分布于水源豐富的綠洲地區(qū)，這些地區(qū)海拔相對較低，土壤質(zhì)量相對較好。NPP值小于100 g C/（m2·a）的區(qū)域主要位于喀什地區(qū)北部土壤鹽漬化嚴(yán)重區(qū)和南部高海拔地區(qū)[24-25]，這些區(qū)域占總NPP面積的39.162%。NPP在100～200 g C/（m2·a）的區(qū)域面積最大，占總NPP面積的43.203%，主要分布在葉爾羌河中上游、中游兩岸以及喀什噶爾綠洲區(qū)。NPP值大于300 g C/（m2·a）的區(qū)域面積很少，集中分布于喀什地區(qū)最西邊。對NPP變化結(jié)果進(jìn)行顯著性檢驗，如表2所示，2000—2020年喀什地區(qū)NPP顯著性變化面積大于不顯著變化的面積。其中，顯著增加區(qū)域面積占全區(qū)NPP面積的57.119%，不顯著增加的面積占全區(qū)NPP面積的31.522%，顯著減少的區(qū)域所占比例最少，只有2.558%。

2000—2020年，喀什地區(qū)NPP均值變化如圖4所示。由NPP年際變化可知，喀什地區(qū)NPP持續(xù)向好的地區(qū)主要集中于疏附縣、澤普縣和葉城縣中部等，其余區(qū)域變化不顯著。基于ArcGIS軟件，對每5年一個時間周期的NPP進(jìn)行差別檢驗。結(jié)果表明，2000—2005年，NPP持續(xù)增長；2006—2010年，喀什地區(qū)南部、西南部NPP有些許減少，其余地區(qū)繼續(xù)增長；2011—2015年和2016—2020年，葉爾羌河流域周邊及水系發(fā)達(dá)區(qū)域NPP持續(xù)增長，喀什地區(qū)其余區(qū)域NPP增長都受限。

利用GeoDa軟件進(jìn)行喀什地區(qū)近20 a NPP空間自相關(guān)分析，顯著性系數(shù)P為0.001，遠(yuǎn)小于0.05，表明莫蘭指數(shù)顯著有效。全局莫蘭指數(shù)為0.910，表明喀什地區(qū)植被NPP的空間自相關(guān)性顯著。對局部莫蘭指數(shù)進(jìn)行檢驗，結(jié)果表明，高-高聚集的NPP主要集中分布于喀什地區(qū)西部、中部，這些地區(qū)都是水資源十分豐富的區(qū)域；低-低聚集的NPP主要分布于土壤鹽漬化和污染嚴(yán)重的喀什地區(qū)北部及其南部高原山區(qū)。

2.3 NPP驅(qū)動因素分析

利用GS＋軟件構(gòu)建半變異函數(shù)模型，并選用球狀模型。在理論模型中，塊金值表示由隨機(jī)因素引起的變異，基臺值表示系統(tǒng)內(nèi)的總變異，塊金值與基臺值的比值代表由地形、氣候等非人為的結(jié)構(gòu)因素所引起的變異占總空間變異的比例。比值越小，非人為因素作用越大[26]。將2000—2020年分成4個時間段進(jìn)行具體研究，半變異函數(shù)模型的相關(guān)參數(shù)如表3所示。結(jié)果表明，所有時間段選用球狀模型的擬合效果都是最佳，其決定系數(shù)都接近1，即此模型解釋因變量的能力最強(qiáng)。塊金值與基臺值的比值都接近25%，比值最小為24.43%，比值最大為27.78%。結(jié)果充分說明，無論是長期分析還是短期分析，喀什地區(qū)近20 a植被NPP的空間異質(zhì)性都主要是由結(jié)構(gòu)型因素導(dǎo)致的。基于此，在接下來的驅(qū)動因素分析中，將考慮結(jié)構(gòu)型因素的作用。

2.3.1 NPP空間分布與海拔的關(guān)系

利用自然斷點法將喀什地區(qū)海拔分為9級，分析不同海拔下植被NPP均值與所占面積的變化，如圖5所示。海拔為1 034～1 425 m時，NPP所占面積最大，占全區(qū)NPP面積的69.881%。海拔為1 425～2 016 m時，NPP均值最高，為255 g C/（m2·a）。NPP主要分布于海拔較低的區(qū)域，隨著海拔的升高，NPP均值與NPP所占面積都逐漸減少。

2.3.2 NPP空間分布與坡度的關(guān)系

采用自然斷點法將坡度分為9級，分析坡度與植被NPP的關(guān)系，如圖6所示。坡度小于4°時，NPP所占面積最大，占全區(qū)NPP面積的64.677%。之后，隨著坡度的增加，NPP所占面積急劇減小。坡度為4°～11°時，NPP均值最高，為157 g C/（m2·a），此值相比其他坡度的NPP值稍高。之后，隨著坡度增加，NPP均值逐漸減小。

2.3.3 NPP空間分布與降水量的關(guān)系

2001—2020年，喀什地區(qū)年降水量介于14.834～317.308 mm，多年平均降水量為52.166 mm。由圖7（a）可知，近20 a NPP均值與降水量均值的變化趨勢相似。NPP（Y）與年降水量（X2）的線性擬合分析結(jié)果如圖7（b）所示，二者的相關(guān)系數(shù)R02為0.356。借助SPSS軟件對其進(jìn)行相關(guān)分析，置信度小于0.05時，多年降水量與NPP的皮爾遜相關(guān)系數(shù)為0.597，表明降水量與植被NPP正相關(guān)。從逐像元尺度分析喀什地區(qū)NPP與降水量的相關(guān)性，由圖8可知，相關(guān)系數(shù)最低值為-0.669，最高值為0.854，呈正相關(guān)的區(qū)域面積遠(yuǎn)遠(yuǎn)大于負(fù)相關(guān)區(qū)域，且廣泛分布于喀什地區(qū)，呈顯著正相關(guān)的區(qū)域面積占全區(qū)NPP面積的30.622%。相關(guān)系數(shù)最高的區(qū)域是帕米爾高原山麓前沿地區(qū)；相關(guān)系數(shù)最低的區(qū)域是塔什庫爾干塔吉克自治縣北部，這里海拔較高，不適宜植被生長。

2.4 討論

植被凈初級生產(chǎn)力是陸地生態(tài)系統(tǒng)中生態(tài)功能和碳循環(huán)的重要評價指標(biāo)，對探究NPP時空變化及其驅(qū)動因素具有重要意義[27]。本研究分析喀什地區(qū)2000—2020年NPP時空變化特征，并確定地統(tǒng)計分析中降水量、海拔、坡度等因素與NPP的關(guān)系。研究發(fā)現(xiàn)，喀什地區(qū)近20 a植被NPP總體呈增長趨勢，具有空間異質(zhì)性與空間自相關(guān)性，結(jié)構(gòu)型因素在NPP時空變化中起重要作用。

本研究使用的NPP數(shù)據(jù)來源于美國國家航空航天局的MODIS17A3數(shù)據(jù)集，該數(shù)據(jù)集已廣泛應(yīng)用于區(qū)域及全球碳循環(huán)模擬、生態(tài)系統(tǒng)研究等。眾多學(xué)者也將其用于新疆植被NPP、草地NPP的研究。姜萍等[28]直接利用2000—2017年中分辨率成像光譜儀（moderate-resolution imaging spectroradiometer，MODIS）數(shù)據(jù)，研究新疆植被NPP及其對氣候變化響應(yīng)的海拔分異。高原[29]也將MODIS數(shù)據(jù)用于研究新疆NPP的時空變化。對比發(fā)現(xiàn)，本研究結(jié)果與學(xué)者研究中關(guān)于喀什地區(qū)的NPP結(jié)論是一致的。同時，本研究結(jié)果與張山清等[30]對新疆自然植被凈第一性生產(chǎn)力的研究中關(guān)于喀什地區(qū)的描述是統(tǒng)一的，與盧娜等[31]關(guān)于南疆植被總初級生產(chǎn)力與區(qū)域降水的研究結(jié)論相符。這些充分證明本研究的正確性與科學(xué)性。植被作為自然景觀最直接的反映，也是評估生態(tài)環(huán)境質(zhì)量的指標(biāo)之一，其分布特征、影響因素等都引起國內(nèi)外學(xué)者的廣泛關(guān)注。本研究的空間尺度不大，MODIS重采樣后的像元大小為250 m×250 m，可以較好反映該區(qū)域植被NPP的生長狀況與演變特征。在喀什地區(qū)NPP的時空變化中，自然要素發(fā)揮主導(dǎo)作用，但近10 a NPP的增長趨勢也表明，環(huán)境保護(hù)政策、生態(tài)建設(shè)等人為因素的作用也在不斷提高。

3 結(jié)論

本文以2000—2020年為研究時段，基于NPP數(shù)據(jù)、氣象數(shù)據(jù)和高程數(shù)據(jù)等，利用Mann-Kendall檢驗、趨勢分析、地統(tǒng)計分析和相關(guān)性分析等方法對喀什地區(qū)植被NPP時空變化及其驅(qū)動因素進(jìn)行研究。結(jié)果表明，2000—2020年喀什地區(qū)植被NPP年均值整體呈上升趨勢。2010年是一個突變時間點，2010年以前植被NPP均值基本低于近20 a的平均值，2010年以后NPP均值基本高于近20 a的年均值。它與降水量的關(guān)系十分密切。近20 a喀什地區(qū)植被NPP在空間分布上具有明顯的差異性，存在空間自相關(guān)性和空間異質(zhì)性，NPP高值較為集中地分布于喀什地區(qū)中部沿葉爾羌河的綠洲帶以及溝渠、水庫十分密集的西部。自然因素對NPP的影響遠(yuǎn)大于人為因素。自然因素中，海拔、坡度和降水量對喀什地區(qū)植被NPP的分布都有驅(qū)動影響。喀什地區(qū)NPP近20 a呈增長趨勢，與喀什地區(qū)降水量呈正相關(guān)。除此之外，生態(tài)保護(hù)政策也初顯成效，尤其是近10 a NPP增長趨勢明顯，這表明生態(tài)文明建設(shè)的效應(yīng)開始顯現(xiàn)。除自然因素外，人為因素的作用也不可忽視。研究結(jié)果可為改善喀什地區(qū)脆弱的生態(tài)環(huán)境和促進(jìn)社會經(jīng)濟(jì)的持續(xù)穩(wěn)定發(fā)展提供參考。本研究分析喀什地區(qū)近20 a植被NPP的時空變化特征及其驅(qū)動因素，取得顯著成果，但仍存在部分亟待改進(jìn)之處。總體來說，選取的因素不夠全面，后續(xù)研究還可以考慮太陽輻射、土壤類型等自然因素，兼顧人口、經(jīng)濟(jì)和環(huán)境保護(hù)政策等人為因素，利用地理探測器、神經(jīng)網(wǎng)絡(luò)等方法，將自然因素與人為因素相結(jié)合，定量分析及比較各因素對植被NPP時空變化的影響。

參考文獻(xiàn)

1 SALLABA F，LEHSTEN D，SEAQUIST J，et al.A rapid NPP meta-model for current and future climate and CO2 scenarios in Europe[J].Ecological Modelling，2015，302：29-41.

2 陳卓奇，邵全琴，劉紀(jì)遠(yuǎn)，等.基于MODIS的青藏高原植被凈初級生產(chǎn)力研究[J].中國科學(xué)（地球科學(xué)），2012（3）：402-410.

3 趙東升，吳紹洪，尹云鶴.氣候變化情景下中國自然植被凈初級生產(chǎn)力分布[J].應(yīng)用生態(tài)學(xué)報，2011（4）：897-904.

4 梁妙玲，謝正輝.我國氣候?qū)χ脖环植己蛢舫跫壣a(chǎn)力影響的數(shù)值模擬[J].氣候與環(huán)境研究，2006（5）：582-592.

5 嚴(yán)俊霞，黃 浩，高彥華，等.吉林省長時間序列陸地生態(tài)系統(tǒng)NPP時空變化特征[J].水土保持學(xué)報，2021（5）：172-180.

6 Potter C S，Randerson J T，F(xiàn)ield C B，et al.Terrestrial ecosystem production：a process model based on global satellite and surface data[J].Global Biogeochemical Cycles，1993（4）：811-841.

7 朱文泉，陳云浩，徐 丹，等.陸地植被凈初級生產(chǎn)力計算模型研究進(jìn)展[J].生態(tài)學(xué)雜志，2005（3）：296-300.

8 朱文泉，潘耀忠，陽小瓊，等.氣候變化對中國陸地植被凈初級生產(chǎn)力的影響分析[J].科學(xué)通報，2007（21）：2535-2541.

9 朱文泉，潘耀忠，張錦水.中國陸地植被凈初級生產(chǎn)力遙感估算[J].植物生態(tài)學(xué)報，2007（3）：413-424.

10 洪辛茜，黃 勇，孫 濤.我國西南喀斯特地區(qū)2001—2018年植被凈初級生產(chǎn)力時空演變研究[J].生態(tài)學(xué)報，2021（24）：1-11.

11 周雷雷，鄭詩軍，尹 捷，等.以“胡煥庸線”為界的中國東西部凈初級生產(chǎn)力變化分析[J].遙感技術(shù)與應(yīng)用，2021（4）：916-925.

12 張 偉，楊淑敏，周向玲.喀什地區(qū)沙塵暴天氣變化特征及影響因素分析[J].干旱區(qū)資源與環(huán)境，2016（6）：95-101.

13 謝 培，顧艷玲，張玉虎，等.1961—2015年新疆降水及干旱特征分析[J].干旱區(qū)地理，2017（2）：332-339.

14 李 路，孫桂麗，陸海燕，等.喀什綠洲土地利用空間格局變化特征分析[J].西南大學(xué)學(xué)報（自然科學(xué)版），2020（5）：141-150.

15 朱鵬航，于瑞宏，葛 錚，等.烏梁素海長時序水質(zhì)變化及其驅(qū)動因子[J].生態(tài)學(xué)雜志，2022（3）：546-553.

16 周妍妍.人類活動對疏勒河流域景觀格局和凈初級生產(chǎn)力的影響研究[D].蘭州：蘭州大學(xué)，2019：11-12.

17 劉立文，段永紅，徐立帥，等.黃土丘陵區(qū)耕地細(xì)碎化特征及影響因素研究：以山西省和順縣為例[J].中國農(nóng)學(xué)通報，2021（22）：90-96.

18 喬 佳，寧小李，鄧紅蒂，等.城市郊區(qū)耕地細(xì)碎化及其空間自相關(guān)性分析：以河北省保定市新市區(qū)為例[J].國土資源科技管理，2015（2）：9-15.

19 肖衛(wèi)東.中國種植業(yè)地理集聚的空間統(tǒng)計分析[J].經(jīng)濟(jì)地理，2014（9）：124-129.

20 王杰云，羅志軍，齊 松.城鎮(zhèn)空間擴(kuò)張與景觀生態(tài)風(fēng)險的耦合關(guān)聯(lián)：以江西省袁州區(qū)為例[J].水土保持研究，2021（5）：142-151.

21 木合塔爾·艾買提，迪里木拉提·玉蘇甫.烏什縣植棉田土壤養(yǎng)分空間分級特征分析[J].湖北農(nóng)業(yè)科學(xué)，2019（2）：58-64.

22 李恒凱，歐 彬，劉雨婷.基于MOD17A3的南嶺山地森林區(qū)植被NPP時空分異分析[J].西北林學(xué)院學(xué)報，2017（6）：197-202.

23 許方岳，焦鴻渤，丁雪丹，等.亞熱帶常綠闊葉林植被凈初級生產(chǎn)力時空特征[J].西北林學(xué)院學(xué)報，2019（2）：62-68.

24 李有文，莫治新，薛江鵬，等.喀什地區(qū)土壤重金屬污染評價及來源解析[J].干旱區(qū)資源與環(huán)境，2020（8）：147-153.

25 張壽雨，吳世新，賀 可.基于開墾年齡的新疆鹽漬化耕地時空特征分析[J].干旱區(qū)研究，2017（5）：972-979.

26 王雪梅，柴仲平，武紅旗.典型干旱荒漠綠洲區(qū)耕層土壤養(yǎng)分空間變異[J].水土保持通報，2016（1）：51-56.

27 高艷妮，于貴瑞，張 黎，等.中國陸地生態(tài)系統(tǒng)凈初級生產(chǎn)力變化特征：基于過程模型和遙感模型的評估結(jié)果[J].地理科學(xué)進(jìn)展，2012（1）：109-117.

28 姜 萍，丁文廣，肖 靜，等.新疆植被NPP及其對氣候變化響應(yīng)的海拔分異[J].干旱區(qū)地理，2021（3）：849-857.

29 高 原.基于遙感的新疆NPP時空變化特征及其影響因素分析[D].蘭州：西北師范大學(xué)，2015：15-16.

30 張山清，普宗朝，伏曉慧，等.氣候變化對新疆自然植被凈第一性生產(chǎn)力的影響[J].干旱區(qū)研究，2010（6）：905-914.

31 盧 娜，牛 俊.南疆植被總初級生產(chǎn)力對區(qū)域降水的響應(yīng)規(guī)律[J].中國農(nóng)業(yè)大學(xué)學(xué)報，2020（7）：103-110.