曾 迅

（南京農(nóng)業(yè)大學(xué)，南京 210095）

糧食安全是關(guān)乎國(guó)計(jì)民生的重大戰(zhàn)略問題。對(duì)區(qū)域耕地生產(chǎn)力進(jìn)行科學(xué)評(píng)價(jià)，研究耕地生產(chǎn)力的時(shí)空演變規(guī)律對(duì)于挖掘耕地生產(chǎn)潛力、保障糧食安全等具有重大價(jià)值。傳統(tǒng)耕地生產(chǎn)力評(píng)價(jià)多依據(jù)采樣數(shù)據(jù)和社會(huì)統(tǒng)計(jì)數(shù)據(jù)[1-3]。隨著遙感技術(shù)的發(fā)展，遙感數(shù)據(jù)逐步應(yīng)用生產(chǎn)力評(píng)價(jià)領(lǐng)域。已有學(xué)者利用歸一化植被指數(shù)（NDVI）、增強(qiáng)型植被指數(shù)（EVI）、凈初級(jí)生產(chǎn)力（NPP）和總初級(jí)生產(chǎn)力（GPP）[4-6]等指標(biāo)進(jìn)行生產(chǎn)力評(píng)估，或通過構(gòu)建CASA模型、GLO-PEM模型、PSN模型及VPM模型等進(jìn)行生產(chǎn)力計(jì)算[5]。其中NPP數(shù)據(jù)在實(shí)時(shí)、連續(xù)監(jiān)測(cè)和年際波動(dòng)和長(zhǎng)期變化趨勢(shì)的探測(cè)方面優(yōu)勢(shì)顯著[7]，能較好地表征區(qū)域耕地生產(chǎn)力。

本文基于NASA平臺(tái)提供的MODIS-NPP數(shù)據(jù)及土地覆蓋數(shù)據(jù)，通過空間自相關(guān)分析法、變異系數(shù)法及Sen-Man Kendall趨勢(shì)檢驗(yàn)法，對(duì)北京市2001—2005年耕地生產(chǎn)力進(jìn)行時(shí)空特征分析，揭示該時(shí)間段北京市耕地生產(chǎn)力的空間分布情況、穩(wěn)定性及變化趨勢(shì)。

1 研究區(qū)域、數(shù)據(jù)與方法

1.1 研究區(qū)域

北京市位于華北平原北部，北緯39°56′，東經(jīng)116°20′，總面積16 410.54 km2，下轄16個(gè)區(qū)縣。城市地形呈現(xiàn)“西北高，東南低”的態(tài)勢(shì)，平均海拔為43.5 m。氣候?yàn)闇貛О敫珊蛋霛駶?rùn)大陸性季風(fēng)氣候，土壤類型以棕壤、褐土和潮土等為主[8]。

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

（1）NPP數(shù)據(jù)集。采用NASA官網(wǎng)提供的MOD17A 3HGF數(shù)據(jù)產(chǎn)品，空間分辨率為500 m，從中提取北京市2001—2005年NPP時(shí)間序列數(shù)據(jù)。

（2）土地覆蓋數(shù)據(jù)。采用NASA官網(wǎng)提供的MCD12 Q1數(shù)據(jù)產(chǎn)品，空間分辨率為500 m，從中提取北京市2001—2005年土地覆蓋數(shù)據(jù)。

（3）糧食產(chǎn)量數(shù)據(jù)。來源于《北京市統(tǒng)計(jì)年鑒》。

1.3 研究方法

首先，從MOD17A3HGF數(shù)據(jù)集中提取2001—2005年北京市耕地年均MODIS-NPP數(shù)據(jù)，通過線性回歸模型擬合耕地NPP數(shù)據(jù)與耕地作物產(chǎn)量之間的關(guān)系，驗(yàn)證耕地NPP數(shù)據(jù)對(duì)耕地生產(chǎn)力的表征能力。而后通過空間自相關(guān)分析法、變異系數(shù)法、Sen-Mann Kendall法分析北京市耕地生產(chǎn)力的空間集聚特征、穩(wěn)定性及變化趨勢(shì)。技術(shù)路線如圖1所示。

圖1 技術(shù)路線

1.3.1 耕地生產(chǎn)力空間集聚特征研究

對(duì)于耕地生產(chǎn)力的空間集聚特征研究，本文先采用全局Moran's I系數(shù)判斷耕地NPP是否在空間上存在聚集性，而后利用ArcGIS中的“聚類和異常值分析”進(jìn)行對(duì)北京市耕地NPP進(jìn)行局部自相關(guān)分析，探究耕地生產(chǎn)力的局部空間集聚。

1.3.2 耕地生產(chǎn)力穩(wěn)定性研究

對(duì)于耕地生產(chǎn)力穩(wěn)定性研究，本文采用CV表示耕地NPP在一段時(shí)間內(nèi)的波動(dòng)程度，基于CV數(shù)值利用自然斷點(diǎn)法將耕地生產(chǎn)力穩(wěn)定性分為低波動(dòng)、較低波動(dòng)、較高波動(dòng)及高波動(dòng)4個(gè)等級(jí)。

變異系數(shù)Cv的計(jì)算公式為：

式中：CV為變異系數(shù)；NPPi為時(shí)間序列i的NPP數(shù)值：NPP為整段時(shí)間序列中NPP的均值；n為NPP時(shí)間序列中數(shù)據(jù)的個(gè)數(shù)。

1.3.3 耕地生產(chǎn)力變化趨勢(shì)研究

對(duì)于耕地生產(chǎn)力的變化趨勢(shì)研究，本文采用Sen-Mann Kendall趨勢(shì)檢驗(yàn)法。先進(jìn)行Sen斜率估計(jì)，再進(jìn)行Mann-Kendall（MK）檢驗(yàn)時(shí)間序列趨勢(shì)的顯著性，其中Sen斜率估計(jì)的斜率β計(jì)算公式為

式中：NPPi和NPPj分別為時(shí)間序列i和j中的數(shù)值。

對(duì)于Mann-Kendall（MK）檢驗(yàn)，統(tǒng)計(jì)量S定義為

S的方差為

式中：n為NPP時(shí)間序列中數(shù)據(jù)的個(gè)數(shù)。

標(biāo)準(zhǔn)化后的檢驗(yàn)統(tǒng)計(jì)量Z計(jì)算式為

在α=0.05的置信水平上判斷NPP時(shí)間序列變化趨勢(shì)的顯著性。

2 結(jié)果與分析

將北京市2001—2005年的北京市各區(qū)縣糧食產(chǎn)量與對(duì)應(yīng)耕地NPP年均數(shù)據(jù)進(jìn)行回歸分析，實(shí)現(xiàn)耕地生產(chǎn)力到NPP的轉(zhuǎn)換。結(jié)果表明，R2在0.568 7與0.770 1之間，北京市NPP數(shù)據(jù)能較好地表征耕地生產(chǎn)力。

從時(shí)間變化上看，2001—2005年北京市耕地NPP年際變化整體呈現(xiàn)上升趨勢(shì)，5年內(nèi)的均值為403.64gC/m2，如圖2所示。從空間分布上看，北京市耕地多年NPP空間分布總體上呈現(xiàn)北部高，南部低，由北向南遞減的特征。

圖2 北京市2001—2005年耕地年際NPP變化情況

2.1 北京市耕地生產(chǎn)力空間集聚特征分析

利用ArcGIS10.6對(duì)北京市2001—2005年耕地平均NPP進(jìn)行全局空間自相關(guān)分析。結(jié)果表明，NPP的全局Moran's I值大于0，Z得分169.86，表明北京市耕地NPP存在正向的全局空間自相關(guān)性。再對(duì)北京市耕地NPP進(jìn)行局部自相關(guān)分析，顯示為“高—高集聚”的耕地占全市總耕地的23.8%，主要分布于西北部、西南部及中部地區(qū)，說明該區(qū)域耕地NPP較高的地方，周圍耕地NPP也較高。顯示為“低—低集聚”的耕地占全市總耕地的28.6%，主要分布于北京市的南部及東南部地區(qū)，說明該區(qū)域耕地NPP較低的地方，周圍耕地NPP也較低。

2.2 北京市耕地生產(chǎn)力穩(wěn)定性分析

計(jì)算2001—2005年北京市耕地生產(chǎn)力的變異系數(shù)CV，結(jié)果顯示，北京市生產(chǎn)力CV分布在0～1.51之間。耕地生產(chǎn)力穩(wěn)定性基于自然斷點(diǎn)法可分為4個(gè)等級(jí):低波動(dòng)區(qū)（CV≤0.16）、較低波動(dòng)區(qū)（0.16＜CV≤0.32）、較高波動(dòng)區(qū)（0.32＜CV≤0.89）、高波動(dòng)區(qū)（CV＞0.89）。北京市耕地生產(chǎn)力穩(wěn)定性存在顯著的空間差異，變異系數(shù)總體上呈現(xiàn)出東南高，西北低的態(tài)勢(shì)。其中，低波動(dòng)區(qū)占全市耕地面積的43.9%，主要分布于延慶縣、密云縣和順義區(qū)等北京市北部和中部地區(qū);較低波動(dòng)區(qū)占全市耕地面積的52.9%，主要分布于通州區(qū)、大興區(qū)、順義區(qū)等北京市南部;較高波動(dòng)區(qū)和高波動(dòng)區(qū)分別占全市耕地面積的1.8%、1.4%，零星分布于北京市各區(qū)縣。總體而言，北京市2001—2005年耕地生產(chǎn)力較為穩(wěn)定。

2.3 北京市耕地生產(chǎn)力變化趨勢(shì)分析

基于MATLAB，采用Sen-Man Kendall趨勢(shì)檢驗(yàn)法對(duì)2001—2005年北京市耕地生產(chǎn)力變化趨勢(shì)進(jìn)行分析。結(jié)果顯示，耕地生產(chǎn)力整體呈改善趨勢(shì)。呈明顯改善趨勢(shì)的耕地占全市耕地的7.3%，主要分布在延慶區(qū)西部，在懷柔區(qū)、密云區(qū)、昌平區(qū)、順義區(qū)及平谷區(qū)也有零散分布。呈輕微改善的耕地占全市耕地的92.5%，遍布北京市各區(qū)耕地。呈輕微退化的耕地僅占全市耕地0.2%，零星分布于順義區(qū)西南部、平谷區(qū)中部及密云區(qū)南部。無呈現(xiàn)嚴(yán)重退化的耕地。

3 結(jié)論

（1）本文利用MODIS-NPP數(shù)據(jù)表征耕地生產(chǎn)力，構(gòu)建了包含空間自相關(guān)分析法、變異系數(shù)法及Sen-Man Kendall趨勢(shì)檢驗(yàn)法的時(shí)空特征分析框架，得到2001—2005年間北京市耕地生產(chǎn)力空間集聚特征、穩(wěn)定性及變化趨勢(shì)。

（2）北京市耕地生產(chǎn)力空間分異顯著，北部和中部耕地生產(chǎn)力較高，且呈現(xiàn)“高—高”集聚狀態(tài)，南部耕地生產(chǎn)力較低，且呈現(xiàn)“低—低”集聚狀態(tài)。北京市耕地生產(chǎn)力整體而言較為穩(wěn)定，低波動(dòng)區(qū)與較低波動(dòng)區(qū)分別占全市耕地面積的43.9%和52.9%。

（3）2001—2005年間，北京市耕地生產(chǎn)力整體呈現(xiàn)改善趨勢(shì)。呈明顯改善趨勢(shì)的耕地占全市耕地的7.3%，主要分布在延慶區(qū)西部，在懷柔區(qū)、密云區(qū)、昌平區(qū)、順義區(qū)及平谷區(qū)也有零散分布。呈輕微改善的耕地占全市耕地的92.5%，遍布北京市各區(qū)耕地。